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Abstract 

This thesis concerns the observation of a new particle and the measurements of its prop¬ 
erties, from the search of the Higgs boson through its decay into two photons at the CMS 
experiment at CERN’s Large Hadron Collider (LHC), on the full LHC “Run I” data col¬ 
lected by the CMS detector during 2011 and 2012, consisting of proton-proton collision 
events at = 7 TeV with L = 5.1 fb“^ and at ^/s = 8 TeV with L = 19.7 fb“^, with 
the hnal calibration. In particular, an excess of events above the background expectation 
is observed, with a local signihcance of 5.7 standard deviations at a mass of 124.7 GeV, 
which constitutes the observation of a new particle through the two photon decay channel. 
A further measurement provides the precise mass of this new particle as 124.72lQ3g GeV = 
124.72]'^Q;32(stat)]'^Q;Jg(syst) GeV. Its total production cross section times two photon decay 
branching ratio relative to that of the Standard Model Higgs boson is determined as 1.12^g'23 
= 1.12]'^Q;2i(stat)]'^Q;g9(syst), compatible with the Higgs boson expectation. Further extrac¬ 
tions of its properties relative to the Higgs boson, including the production cross section 
times decay branching ratios for separate Higgs production processes, couplings to bosons 
and to fermions, and effective couplings to the photon and to the gluon, are all compatible 
with the expectations for the Standard Model Higgs boson. 
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Title: Professor 
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Chapter 1 


Introduction 


The night before Friday, June 15, 2012. I started writing the “unblinding” slides. My 
colleagues from two different teams continued progressing independently towards the final 
plots. This was going to be a sleepless night for these searchers in a working group of the 
Compact Muon Solenoid (CMS) experiment at the Large Hadron Collider (LHC) of the 
European Organization for Nuclear Research (CERN), searching for the Higgs boson [H) 
through its decay into two photons. 

At this point, the Higgs boson remained as the last undetected elementary particle pre¬ 
dicted by the Standard Model (SM) of particle physics The Standard Model endeavors 
to describe the fundamental components of matter and interactions except for gravity—the 
strong, electromagnetic, and weak interactions—in terms of the corresponding elementary 
fields and field quanta, the elementary particles, of which spin-1/2 fermions are matter com¬ 
ponents and spin-1 vector bosons are interaction mediators. The fermions consist of leptons 
and quarks, while the vector bosons consist of gluons for the strong interaction, photons 
for the electromagnetic interaction, and W and Z bosons for the weak interaction. The 
main component of the Standard Model is the electroweak theory, which provides a unified 
description of the electromagnetic interaction and the weak interaction as the electroweak 
interaction. The fundamental assumption underlying this theory is the symmetry between 
these two interactions—the electroweak symmetry. However, the manifest symmetry does 
not allow the particles associated with the interactions to possess mass. This works for the 
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massless photon, but not for the massive W and Z bosons. To resolve this inconsistency and 
to formulate the current form of the electroweak theory in the 1960s, a mechanism invented 
independently by Robert Brout and Frangois Englert, Peter Higgs, and Gerald Guralnik, 
Garl Hagen, and Tom Kibble [6 11 was applied to break the electroweak symmetry. A 
scalar held permeating over space is introduced, which generates masses of W and Z bosons 
by interacting with them. This scalar held could also generate masses of quarks and charged 
leptons through the additional Yukawa interaction. 

The observation of weak neutral currents—mediated by the Z boson—by the Gargamelle 


experiment at GERN in 1973 12,13 , and then the direct observations of the W and Z bosons 


by the UAl and UA2 experiments at GERN’s proton-antiproton collider in 1983 14 16 


conhrmed the prediction by the electroweak theory of the existences of the W and Z bosons. 
And these experimental conhrmations established this theory as the theoretical cornerstone 
of particle physics. Still, the crucial point of the theory lacked experimental evidence—the 
mechanism for the electroweak symmetry breaking and mass generation. Observation of the 
quantum of the scalar held, the scalar boson with spin-0, conventionally called the Higgs 
boson, is the key. 

The search for the Higgs boson had been one of the central tasks of experimental particle 
physicists since the observations of W and Z bosons. The Standard Model predicts the 
couplings of the Higgs boson to the other elementary particles—proportional to the mass 
squared of bosons and just to the mass for fermions—so that its production and decay rates 
at a given mass could be calculated and compared with observations at collider experiments. 
But the mass of the Higgs boson, rriH, is not predicted, which adds to the complications 
of the search. There were indirect constraints on the Higgs mass from the probability con¬ 


servation of WW scattering, niff < ~1 TeV 17 20 , and from the precision electroweak 


measurements, niff < 158 GeV at 95% conhdence level (GL) [^, but still a wide range of 
Higgs mass hypotheses had to be explored. Before the search at the LHG, direct searches at 
GERN’s Large Electron-Positron Gollider (LEP) and Fermilab’s proton-antiproton collider 


Tevatron excluded the mass range mjj < 114.4 GeV 22 and 162 GeV < niff < 166 GeV 23 
respectively, with no evidence of the Higgs boson. 

The LHG was designed to collide two beams of protons composed of quarks bound by 
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gluons, at center-of-mass energies up to a/s = 14 TeV and instantaneous luminosities up to 
Linst = 10^^ cm“^s“^—about 7 times the collision energy and 0(10) times the intensity of the 


Tevatron 24 , the previous most powerful hadron collider—which allows the LHC to create 
particles with masses up to the TeV scale and to relatively quickly accumulate proton-proton 
(pp) collision events for physics analyses—with one potential Higgs boson with a mass of 


125 GeV produced from about 4 billions of inelastic collisions at 7 TeV 25,26 . The design 


and construction of the LHC 27 , along with its two largest experiments CMS 28 and 


ATLAS (A Toroidal LHC Apparatus) each having a comprehensive particle detector 
and a collaboration of thousands of physicists and engineers from all over the world, were 
centered on proving or excluding the existence of the Higgs boson. 


Higgs bosons are produced at the LHC through four major processes from pp collisions: 
gluon fusion (ggH), vector boson {W/Z boson) fusion (VBF), associated production with a 
W or Z boson ( VB), and associated production with a pair of top quarks {itH). Gluon fusion 
is the dominant process. The other three processes occur much less frequently than gluon 
fusion, but with additional particles present along with the Higgs boson, whose features 
are used to identify the Higgs event. The Higgs boson decays immediately—with a lifetime 
about 10“^^ s for a Higgs mass of 125 GeV. The Higgs search is therefore conducted through 
its decay channels as well as its production processes. There are hve main decay channels 
in terms of the sensitivity to the Higgs search: Higgs decaying into two photons {H — 77 ), 
two Z bosons to four charged leptons (ZZ —)■ 4£) (the charged lepton here refers to electron 
or muon), two W bosons to two charged leptons and two neutral leptons—two neutrinos 
{W^W~ —)■ 2i2i'), and either two tau leptons or two bottom quarks or bb). The 

77 —)■ 77 channel—having a hnal state of two energetic and isolated photons which are 
clearly identihed and whose energies are measured with excellent resolution—is one of the 
most promising channels in the search range 110 GeV < mu < 150 GeV. The diphoton 
signature allows the reconstruction of a narrow signal peak in the diphoton mass {m^^) 
spectrum—corresponding to the Higgs resonance with a natural width negligible relative to 
the detection resolution—on top of a smoothly falling background from known SM physics 
processes, which yields eloquent evidence if the Higgs boson exists. 


About one year ago before this June night, when the LHG had just ramped up the 


21 







luminosity of pp collisions at ^/s = 7 TeV to produce significant amount of data for physics 
analyses, I came to CERN to work on this thesis—searching for the Higgs boson through the 
two photon decay channel by analyzing the data collected by the CMS detector, together 
with my colleagues from MIT and the 77 —)■ 77 working group of the CMS collaboration. 
Despite the appealing feature of a signal peak in the diphoton mass spectrum, the two 
photon decays are very rare—about one out of hve hundred decays from the already rarely 
produced Higgs boson, assuming a mass of 125 GeV. The great challenge facing this channel 
is to identify the small peak from a background that is several orders of magnitude larger. 

To fully unfold the power of the diphoton mass spectrum, we made the analysis strat¬ 
egy to classify diphoton events according to the expected signal-to-background ratio {S/B) 
under the peak assuming the existence of the SM Higgs boson. Specihcally, we developed 
the analysis using advanced multivariate analysis (MVA) techniques to fold in all the rele¬ 
vant diphoton information of an event—variables related to photon identihcation, diphoton 
kinematics and mass resolution—into a single event classiher, and used it to optimize the 
classihcation of the events. The MVA analysis significantly improved the expected Higgs 
sensitivity, equivalent to adding about more than 40% of the data, with respect to our initial 
analysis using the traditional cut-based techniques, which selected and classihed diphoton 
events by applying simple cuts on a subset of MVA input variables. We therefore used the 
MVA analysis as our main analysis, with the cut-based analysis as a cross-check later. For 
both analyses, the additional feature of the VBF Higgs production process—a pair of jets 
fragmented from a pair of quarks present in the final state along with the two photons—was 
utilized to further select events into high S/B classes. Finally, any possible Higgs signal 
of a mass in 110 GeV < rriH < 150 GeV was extracted by a simultaneous likelihood £t to 
the reconstructed diphoton mass spectra over 100 GeV < < 180 GeV of all the event 

classes, using parametric signal and background models. For each event class, the signal 
model for any Higgs mass was derived from simulation, 

and the data/simulation discrepancies are corrected and validated through control sam¬ 
ples. The background model was derived directly from the data utilizing the smoothly-falling 
nature of the background shape. The expected background under the emerging signal peak 
for any Higgs mass was constrained by the large number of events in the diphoton mass 
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sidebands of the signal region. To cross-check the background modeling, we used an alter¬ 
native MVA analysis, which extracted the signal by counting events in the signal region—in 
classes dehned by both the diphoton event classiher used in the main MVA analysis and the 
diphoton mass—and was less sensitive to the exact shape of the background. 

By early 2012, we observed an excess of events above the background expectation with 
a local signihcance of about 3 standard deviations at a mass around 125 GeV in if —)■ 77 
channel, from both the cut-based analysis and the newly developed MVA analyses, on the 
2011 dataset collected by the CMS detector from pp collision events at ^/s = 7 TeV with 
L = 4.8 fb“^ (1 barn (b) equals 10“^^ cm^). Taking into account the probability that 
the background fluctuated at any mass point within our entire search range, the global 


signihcance was below 2 standard deviations 30,31 . Among other search channels, the 


if —> ZZ —>■ 4i channel observed an excess near 120 GeV but not as signihcant 32 . And 


the CMS combined result of the hve main decay channels was driven by if —)■ 77 , showing 


an excess with a local signihcance just above 3 standard deviations 33 . At the same time. 


the ATLAS experiment also observed an excess with a local signihcance above 3 standard 
deviations near 125 GeV from the combined search, driven by its two photon channel as 


well, and with a global signihcance of about 2 standard deviations 34 


The excess in if —)• 77 channel remained when we rerun our analyses on the full 2011 
dataset with the integrated luminosity increased to L = 5.1 fb“^ and with improved detector 
calibration. To determine the source of the excess—an upward huctuation of background vs. 
a real signal—the analysis of the 2012 data was critical. We improved and re-optimized our 
analyses to accommodate the enhanced collision energy to i/s = 8 TeV and the increased 
overlapping pp collisions. To avoid the possibility of biasing the results, we developed the 
analysis in a strict “blind” manner, i.e. we did not look at the diphoton mass spectrum or 
extract the observed results in the potential signal region 110 GeV < < 150 GeV until 

our analysis procedure was hxed and fully verihed. All the other Higgs search channels also 
progressed with a “blind” procedure as well. 


Finally, we reached this night before June 15. We had gotten our analyses pre-approved by 
the collaboration one week ago, by showing the performance of the various components and 
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the expected results of the analyses—both the expected exclusion limit of the signal strength 
under the background-only hypothesis, and the signihcance of the excess over background 
assuming the existence of the SM Higgs boson—on the first 2012 data with L = 1.5 fb“^ 
combined with the 2011 data. And we just obtained the “green light” this afternoon after a 
further review to unblind the cut-based analysis on the 2012 data with L = 3.1 fb“^ certihed 
right after the pre-approval. The unblinding of our MVA analyses—the ones into which we 
were putting the most effort—was postponed until the next week, as we decided to wait for 
an updated simulation sample important for “training” the MVA event classiher earlier this 
morning around 3:00 am. 

I stayed in a room of the CERN hostel this night, within three minutes walk from my 
office at CERN Building 32. The past couple of weeks were so intense—producing the results 
for the pre-approval and then working for the “unblinding”—that I almost lived at CERN, 
with few hours of sleep at the desk plus ~10 cups of coffee many of the days. What was I 
feeling when I wrote down the title “Search For A Standard Model Higgs Boson Decaying 
Into Two Photons—Unblinding”? Was my heart beating as fast as I am feeling now? Or 
even faster? Representing the if —)■ 77 gronp, I was going to nnblind onr Higgs search results 
to the entire CMS collaboration in the coming afternoon—the resnlts still sitting in the dark 
waiting to be uncovered. 

To reach this night, we had gone through many sleepless nights working on the analyses 
over the past year—from the development of different analysis ingredients, the various cor¬ 
rections and validations for the signal modeling, the large amount of tests and justihcations 
of the backgronnd model, to the multiple rounds of producing results on the dataset to keep 
updated with the increased luminosity and improved calibrations—accompanied by countless 
meetings, presentations, emails, messages, discussions, and also multiple rounds of docnmen- 
tation preparations. Onr different teams, running independent analysis frameworks, had also 
gone throngh constrnctive while sometimes herce competitions on the analysis methods, but 
ultimately to striving together for the hnal resnlts. 

The nights were mixed with mornings for the past few days. To re-optimize and hnalize 
the analysis ingredients for the “unblinding” data condition, validate the inputs and ontputs, 
process the datasets and simulation samples to select events and compute variables for the 
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final event classification and diphoton mass spectrum reconstruction, and synchronize the 
event selection and variable computation among different teams—within a week—equals to 
a huge effort of the whole group. In the end, two teams synchronized to a satisfactory level. 
Each team was going to produce a set of “unblinded” results to cross-check. 

The night deepened, folding all the sleepless nights into the dark sky, towards the un¬ 
known. The slides grew, with analysis descriptions and validation plots, reaching the blank 
region for the hnal plots. The last round of cross-checks started between the two teams, with 
information flowing through an email thread. Finally, the expected results and the event 
yields agreed. Time to look into the signal region. Around 3:00 am, both teams unblinded 
the diphoton mass spectra of the 2012 dataset- 


Clear excesses jumping from the falling spectra of multiple event classes 

around 125 GeV! 

About the same place of the excess that we observed from the 2011 dataset! 
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“A real signal is there!!!” 


It was not spoken out. 

But I heard 
The yells 
Bursting out, 

From the hearts 
At different ends. 

Enormous excitement 
Flowed out, 

Permeating silently 
In the air. 

I would 

Jump through the window 
Into the sky. 

With the speed of light; 

But in the end. 

Stood together 

With the grand Jura Mountain, 
Quietly upon the ground. 
Watching the pairs of photons 
Passing through layers of nights 
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Later in the morning, we had a quick gathering together with more colleagues from the 
iL —)■ 77 group, in a small meeting room at CERN. Everybody looked extremely excited 
despite not sleeping much. We tried not to speak loudly, since we had to keep the “secret” 
until the “unblinding” event in the afternoon—the event that the unblinded Higgs search 
results from all the main decay channels were presented to the entire collaboration for the 
Erst time. Still more plots to make. We soon went back to the work, with more colleagues 
from different teams joining to help. At this time, all our hearts were bound together. I kept 
modifying the slides with suggestions from my colleagues—except for a short break around 
11:00 am to meet with the CMS management—while new plots continually came, of various 
statistical results or diphoton mass spectra, being updated with hner granularity or refined 
style. As time approached the “unblinding” event, I started putting the final version of the 
plots from my colleagues onto the pages, one after another, each plot a great trust falling 
silently upon my heart. 

Time passed 3:00 pm. I hnally hnished the slides with the last plot just sent from one of 
my teammates. Our other colleagues, after this “super quick and collaborative effort” (quote 
from said teammate), had left earlier to the “unblinding” event, which had already started 
at the CERN “Filtration Plant”. The H —)■ W~^W~ —)■ 2i2u channel was the first to unblind. 
The if —)■ 77 channel was the second, starting at 3:30 pm. We saved the slides onto a flash 
disk and walked quickly to the conference room. Soon, we arrived. My teammate opened 
the door. 

A hot current flowed out. 

The room was packed with CMS colleagues. 

All the seats were taken. 

Many colleagues sat on the floor or stood against the wall. 

There were probably also many colleagues connecting through the video link. 

“Good Luck!” said he. 

One of the “Good Luck” s I having received from my colleagues since the morning. 
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I carried the slides and moved slowly, through the held of colleagues. My body got heavier 


and heavier, as the lives of more and more people connected to my own 

-My teammates always giving the strongest support; My hf —)■ 77 colleagues striving 

together for the unblinding over last night and through this day; And the entire group having 
been working extremely hard together since last year especially during the past couple of 
weeks; Representing whom I was going to unblind our Higgs search results; 

-The CMS colleagues having worked on different stages and aspects of the experiment 

during the past 20 years, leading to the hnal data analysis for the Higgs search—from design, 
construction, commission, to operation; from hardware, software, to computing; from data 
taking, calibration, reconstruction, to validation; The colleagues working on the different 
Higgs decay channels trying to answer the same question; And the colleagues working on the 
different physics topics from the precision measurements within the Standard Model to the 
searches beyond the Standard Model, all trying to deepen and enlarge the same drawing of 
fundamental particle physics; Many of whom were in this room or through the video link, 
waiting to see and to listen to the results; 

-The ATLAS colleagues working towards the same goal; 

-The LHC colleagues providing the most powerful and intense proton beams; 

-The generations of experimental physicists searching for the Higgs boson; 

- The theoretical physicists whose work led to the prediction of this scalar boson 

about half a century ago; 

-And all the physicists from the experimental and theoretical sides working together 

to reach the current understanding of the fundamental components of matter and interactions 
during the past century; 

-And all the human beings, craving to understand the nature and ourselves, asking 

and searching, across the vast space and time ... 

Some of us got together, in this space, at this time. 

The iL —j- 77 presentation was starting. 

My heart was beating violently. My mind was calm. 

“Please everybody, get ready for the next 15 minutes.” 


These 15 minutes would become a part of our common memories 35 












This unblinded if —)■ 77 cut-based results on the combined 2011 (-^/s = 7 TeV, L = 
5.1 fb“^) and 2012 (\/s = 8 TeV, L = 3.1 fb“^) datasets provided the first convincing 
evidence of the existence of a new particle —the local signihcance of the observed 
excess above the expected SM background was about 4 standard deviations at a mass near 
125 GeV. 


Night of the day, CERN Building 32, fourth floor. 
“We just experienced a historic moment.” 

Yes, we did. Not just in the history of science ... 

What else from the day I still remember? 

The smile from the bottom of everyone’s heart. 


Our hnal —)■ 77 results from the main MVA analysis on the combined 2011 and 2012 

datasets, with the 2012 luminosity increased to L = 5.3 fb“^, kept showing an excess of 


events near 125 GeV, with a local signihcance of 4.1 standard deviations 36 . This excess 
was the most signihcant among all the main decay channels, followed by the excess observed 
from the H —)■ ZZ —)■ 4£ channel with a local signihcance of 3.2 standard deviations also 


near 125 GeV 36 . The local signihcance of the observed excess combining the 77 —77 and 
H —)■ ZZ —>■ 4£ channels reached 5.0 standard deviations. The combined signihcance of the 
observed excess of all the hve main decay channels was 4.9 standard deviations (updated 
to 5.0 standard deviations later) near 125 GeV [^. Meanwhile, the ATLAS experiment 
also observed an excess with a local signihcance of 5.0 standard deviations (updated to 5.9 
standard deviations later) near 125 GeV, again with the 77 —)■ 77 providing the largest excess 


with a local signihcance of 4.5 standard deviations 38 


These results led to the announcements of the discovery of a new particle from both 
experiments at GERN in a joint seminar with the 36th International Gonference on High 
Energy Physics (IGHEP) on July 4, 2012. 

This new particle was identihed as a boson with integer spin other than 1 because of its 
decay into two photons. 
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Since the discovery of the new particle, we have continued to verify its observation, and 
to further measure its properties and check its compatibility with the SM Higgs boson, with 
improved inputs and analyses. In particular, we had about three times more 2012 data 
collected by the CMS detector before the end of the LHC “Run I”, with better detector 
calibration and more accurate simulation. We rehned all the major components of the main 
MVA analysis, from the diphoton event classifier, the event classihcation procedure, to the 
modeling of signal and background diphoton mass spectrum. Moreover, we extended the 
analysis to employ the additional features of all the Higgs production processes to select 
events in high S/B classes and to separate signal events from different production processes 
sensitive to Higgs couplings to bosons and to fermions, respectively. These improvements 
signihcantly enhance the Higgs search sensitivity—almost doubling the expected signihcance. 
They allow precise measurement of the mass of the new particle and extraction of its total 
production rate relative to that of the Higgs boson (signal strength). They also allow to 
extract the signal strengths of different Higgs production processes, and to further extract 
the couplings of the new particle to bosons and to fermions relative to those of the Higgs 
boson (coupling strengths). 

This thesis concludes the Odyssey of searching for the Higgs boson through its decay 
into two photons that I have experienced together with my colleagues since 2011, with a 
standalone observation of a new particle and the measurements of its mass, signal strengths, 
and coupling strengths, using the rehned and extended main MVA analysis, on the full LHC 
“Run I” data collected by the CMS detector, consisting of pp collision events at ^/s = 7 TeV 
with L = 5.1 fb“^ in 2011 and at y/s = 8 TeV with L = 19.7 fb“^ in 2012, with the hnal 
calibration. 
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More introduction to the Standard Model and the Higgs boson, the Higgs searches at 
the LHC, and the MVA techniques used in this analysis can be found in the rest of this 
chapter. The introduction of the CMS detector and the event reconstruction is in Chapter 
An overview of this analysis is in Chapter The further descriptions of the analysis 
components are in Chapter |4]j^ The hnal iL —)■ 77 results are in Chapter followed by a 
review of other Higgs results from the CMS and ATLAS experiments in Chapter and the 
conclusion in Chapter [T^ The natural units i.e. h = c = 1 are used throughout this thesis. 


Again, the hnal results from the main MVA analysis are produced and cross-checked by 
two highly synchronized analysis frameworks in the CMS if —?■ 77 group, and cross-checked 
by alternative cut-based and MVA analyses. More details of the main MVA analysis and the 


descriptions of the alternative analyses, are in our analysis note 39 and paper 40 , where 
the results presented are randomly chosen from one of the frameworks. Additional results 
including hypothesis tests between spin -0 and spin -2 models are also in the note/paper, 
which are all consistent with the SM Higgs boson. 


Again, there have been many sleepless nights, which are now only in our memories. 
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1.1 The Standard Model and the Higgs Boson 


The Standard Model Si. based on the relativistic qnantnm gange held theory, describes 
the elementary particles and their interactions except for gravity. Elementary particles are 
depicted as the qnanta of excitation of their corresponding helds, which inclnde spin-1/2 
fermions as fnndamental components of matter and spin-1 vector bosons as mediators of 
interactions. The spin-1/2 fermions consisting of leptons and qnarks are gronped into three 
generations with the higher generation a heavier copy of the lower one, as snmmarized in 


Table The vector bosons mediating three kinds of interactions, weak, electromagnetic, 
and strong—listed in increasing strength—are shown in Table 1.2[ All qnarks and leptons 
participate in weak interactions. The electrically charged particles inclnding charged leptons, 
qnarks, and participate in electromagnetic interactions. Qnarks and glnons, which carry 
color charge, participate in strong interactions. 


Table 1.1: Spin-1/2 fermions: leptons and qnarks (and corresponding anti-particles) in three 
generations. 


Generation 

I 

II 

III 

Leptons 

Electron 

Nentrino 

Ve {Ve) 

Mnon 

Nentrino 


Tan 

Nentrino 



Electron 

e- (e+) 

Mnon 

p- (p+) 

Tan 

T- (r+) 

Qnarks 

Up 

U iu) 

Charm 

c (c) 

Top 

t {t) 

Down 

d {d) 

Strange 

s (s) 

Bottom 

b (b) 


Table 1.2: Spin-1 vector bosons and their corresponding interactions. 


Vector Boson 

Interaction 

W boson 

iu± 

Weak 

Z boson 

z 

Weak 

Photon 

7 

Electromagnetic 

Glnon 

9 

Strong 


The fnndamental mechanism nnderlying the Standard Model is to generate interactions 
by reqniring local gange symmetries. In particnlar, its symmetry gronp is SU{3)c <8) SU{2)l 
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® in which SU{3)c determines the strong interaction while SU(2)l ® U{1)y deter¬ 

mines the electroweak interaction. Though the insight of symmetry enables the derivation 
of the interactions in a systematic way, it forces the bosons mediating the interactions to be 
massless, which is consistent with the massless photon and gluon but apparently not with 
the massive W and Z bosons. A direct breaking of the symmetry would allow for massive 
W and Z bosons but make the theory no longer renormalizable, i.e. the infinities in the 
calculation of observables are not removed. To solve this inconsistency, the Higgs mecha¬ 
nism [6[jII] is employed instead to preserve the renormalizability and at the same time allow 
for massive W and Z bosons. It introduces a doublet of complex scalar fields, which has 
a symmetric potential under SU{2)l ® U{1)y and degenerate vacuum states with non-zero 
expectation values. The SU{2)l ® H(l)y (electroweak) symmetry is spontaneously broken 


by choosing a particular vacuum state, while the renormalizability of the theory is kept 41 


Only one of the four real scalar fields in the doublet remains, which is the Higgs field. W 
and Z bosons then get mass through the interaction with the Higgs field, and the degrees of 
freedom of the three disappearing scalar fields turn into the longitudinal polarizations of W 
and Z. This spontaneous symmetry breaking would also provide mass for fermions, except 
for neutrinos whose mass generation mechanism is unknown, by adding Yukawa interaction 
between fermions and the Higgs field. The particle corresponding to the excitation of the 
Higgs field is the Higgs boson (H). It is neutral, colorless, and has spin (J), parity (P) and 
charge conjugation (C) = O"'"'''. The Standard Model does not predict the mass of the 

Higgs boson but its couplings to bosons and fermions, which are proportional to the boson 
mass squared and to the fermion mass, respectively. With the couplings provided, the Higgs 
cross section for any production process, and its width and corresponding branching ratio 
for any decay mode are predicted for any Higgs mass hypothesis, mu- For more detailed 


introduction on the Standard Model and Higgs boson see References 42 -44 


In case a signal is observed, the Standard Model Higgs production cross sections and 
decay branching ratios at a given mass hypothesis, and its couplings are compared to the 
experimental observations to quantify the compatibility between the signal and the Higgs 
boson. For example, in the search of the Higgs boson through one of its decay modes, 
the compatibility between the observed signal and the Higgs boson is first quantified by 
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extracting its relative total cross section for all production processes times branching ratio 
with respect to the Standard Model Higgs prediction, namely the signal strength, ^h- Given 
sufficient data, the signal strength for each production process is extracted to make a more 
detailed comparison. Depending on the available production processes and the decay mode, 
the compatibility is further quantihed by measuring the relative coupling (coupling strength) 


to bosons, Ky, the relative coupling to fermions, Kf, or both 45 


The last elementary particle of the Standard Model missing experimental conhrmation 
has been the Higgs boson. The search of the Higgs boson is one of the central tasks for 
experimental particle physics, as its experimental observation is crucial to verify the current 
understanding of the electroweak symmetry breaking. Prior to the observation of the Higgs 
boson-like excess in 2012 at the LHC, searches at the Large Electron-Positron Collider (LEP) 
excluded the Standard Model Higgs boson below a mass of 114.4 GeV (95% conhdence 
level) [^. These exclusions were extended in 2012 by searches at the Tevatron, which 
excluded 100 GeV < tuh < 103 GeV and 147 GeV < tuh < 180 GeV (95% conhdence level); 
but also reported a small (3.0 standard deviations) excess at uih = 120 GeV shortly before 
the LHC observation (4^. 


1.2 Search for the Higgs Boson at LHC 


The Large Hadron Collider (LHC) [^, constructed by the European Organization for Nu¬ 
clear Research (CERN), is the highest energy collider of protons (or heavy ions) and allows 
the study of the physics at the TeV scale. Four major experiments are conducted at LHC, 


ALICE 47 , ATLAS 29 , CMS 28 and LHCb 48 . ATLAS and CMS use multi-purpose 


detectors and explore a broad range of particle physics topics, with the search for Higgs 
boson as one of the main goals. 

The LHC is the last element of the CERN accelerator complex as shown in Figure 
49 . It is installed in a circular tunnel with 27 km in circumference, which ranges from 
45 m to 170 m in depth beneath the surface at the outskirts of Geneva. It mainly consists 
of 8 radio frequency cavities for acceleration of each particle beam, 1232 superconducting 
dipole magnets for beam bending, and 392 superconducting quadruple magnets for beam 


34 








focusing. The magnets, cooled by superfluid helium to 1.9 K, are designed to provide a 
magnetic held of 8.33 T. Protons, extracted from hydrogen gas, are hrst accelerated by a 
successive set of accelerators and then injected separately into the two beam pipes of the 
LHC. The two proton beams are designed to run oppositely with 2808 proton bunches per 
beam and about 10^^ protons per bunch, which collide (bunch crossing) every 25 ns at 
center-of-mass energy of up to y/s = 14 TeV and with a peak instantaneous luminosity Ljnst 
= 10^"^ cm“^s“^. The actual bunch crossing rate is every 50 ns, and the collision energy 
and peak instantaneous luminosities are 7 TeV and about 4 x 10^^ cm“^s“^ in 2011, and 
8 TeV and about 8 x 10^^ cm“^s“^ in 2012. The high instantaneous luminosity leads to 
the presence of inelastic pp interactions with low momentum transfer (pileup interactions) in 
the same bunch crossing with the interesting inelastic pp interaction with large momentum 
transfer (hard interaction). The interactions are distributed in space approximately following 
three dimensional Gaussian distribution. The corresponding standard deviation in the beam 
direction and in its perpendicular directions, is about 6 cm (5 cm) for 7 TeV (8 TeV), and 
0(10 /im), respectively. 


Higgs bosons are produced at the LHC through the interactions of the partons from the 
incoming protons. The main production processes, in decreasing order of cross sections, 
are gluon fusion {ggH), vector boson fusion (VBF), associated production with a W or Z 
(WB or ZB, VB for the combined WB and ZB) and associated production with tt (iiB). 
The corresponding leading order Feynman diagrams are shown in Figure and the cross 
section for each process as a function of Higgs mass mn at 7 TeV (8 TeV) is shown on the 


left (right) in Figure 1-3 25,45 


Gluon fusion is the dominant process, whose cross section for tuh = 125 GeV at 8 TeV 
is 19.27 pb, about 7 times the sum of the cross sections of all the other processes. In this 
process, two gluons produces a Higgs boson through a loop of quarks, mainly the heavy top 
quark. This indirect production is due to the fact that the gluon is massless and the Higgs 
boson couples to a boson proportional to its mass squared. The production rate of ggB 
process relative to the Standard Model expectation is proportional to Ky. It is also sensitive 
to the existence of any new colored particles in the loop too heavy to be produced directly. 


which manifests as a deviation of the effective Higgs coupling strength to gluon, Kg 45 
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Figure 1-1: The accelerator complex of CERN. 
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(c) VH 



(d) ttH 


Figure 1-2: Leading order Feynman diagrams for Higgs production processes: (a) ggH (gluon 
fusion) (b) VBF (vector boson fusion) (c) VH (associated production with a IF or Z) (d) 
tiH (associated production with tf). 



Figure 1-3: The Higgs production cross section for each process as a function of Higgs mass 
niH at 7 TeV (left) and 8 TeV (right), along with the theoretical uncertainty bands. From 
top to bottom: ggH, VBF, WH, ZH, iiH. 
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from 1 . 


The cross sections for VBF, VH and iiH processes are mnch smaller than that for ggH 
process, which for mn = 125 GeV at 8 TeV are 1.578 pb, 1.1199 pb and 0.1293 pb, respec¬ 
tively. Despite the low production rates of VBF, VH and itH processes, they are interesting 
processes to be deployed for two reasons. First, in these processes, the Higgs boson is pro¬ 
duced along with other particles whose signature is used to identify the events, and thus 
improves the signal to background ratio. Furthermore, these processes provide additional 
information in Higgs coupling to bosons and fermions. In the VBF process, two quarks 
radiate W bosons or Z bosons, which annihilate to produce the Higgs boson. A pair of 
quarks are present in the hnal state moving oppositely close to the beam direction, which 
fragment into two jets with large opening angle. In the VH process, a quark and an anti¬ 
quark produces a IF or Z boson which in turn radiates a Higgs boson. The W ot Z further 
decays leptonically or hadronically. For the leptonic decay, a lepton (muon or electron) plus 
a neutrino are produced from the W while a pair of leptons are produced from the Z. For 
the hadronic decay, a pair of quarks are produced which fragment to a pair of jets. Both 
the production rates of VBF and VH processes relative to the Standard Model expectation 
are proportional to Ky. In the iiH process, two gluons produce a pair of top and anti-top 
quarks, and a Higgs boson in association. Each top quark decays to a bottom quark plus a 
IF boson. The bottom quark fragments to a so called b-jet, while the IF boson decays in 
the way as mentioned above. The production rate of iiH relative to the Standard Model 
expectation is proportional nj, as that of ggH. 

The Higgs boson—whose lifetime is about 10“^^ s at mn = 125 GeV — decays immediately 
after its production. The Higgs search is therefore conducted through its decay channels as 
well as production processes as explained. The main channels in terms of the sensitivity 
include H —)■ 1F+H^“, H —>■ ZZ, if —)■ 77 where the Higgs boson decays to a pair of bosons, 
and H ^ bb and H —)■ r+r“ where the Higgs boson decays to a pair of fermions. The Higgs 
decay branching ratios in the mn range between 80 GeV and 200 GeV are shown on the 
left in Figure [l-^ j^|^. The H ^ bb channel dominates in the niH range well below the 
WW production threshold. The H —)• H^+IF” channel and H ^ ZZ channel are dominant 
in the rriH range just below and beyond this threshold, because W and Z have much larger 
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mass than the other decay particles and so larger couplings to the Higgs boson. Comparing 
to the other four channels, the if —>■ 77 channel has a much smaller branching ratio across 
the mass range, which reaches the maximum between 120 GeV and 130 GeV and has the 
value 0.228% at rriH = 125 GeV. Despite its small branching ratio, the if —)■ 77 channel has 
a clear signature with two energetic and isolated photons, and allows the reconstruction of 
the narrow Higgs resonance in the diphoton mass spectrum. This makes it one of the most 
sensitive channels for the Higgs discovery in the low mass range, and also one of the only two 
channels—the other is H ^ ZZ with four leptons in the hnal state (ff —?■ ZZ —)■ 4i) —for 
the precision measurement of the Higgs mass. In addition, its production rate is sensitive 
to both the Higgs coupling to bosons and fermions as well as the existence of new charged 
heavy particles. We therefore choose to search for the Higgs boson through the ff —)■ 77 
channel. More details of this channel and our search strategy are given below. 



Figure 1-4: Higgs decay branching ratios for the various channels (left), and total decay 
width (right). 
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1.2.1 Higgs Boson to Two Photons Decay Channel 

Signal and Backgronnd 


The Higgs boson decays to two photons through a loop of massive charged particles, mainly 
W boson and top quark, since the photon is massless while the Higgs boson only couples to 


massive particles. The leading order Feynman diagrams are shown in Figure where the 
W loop and top quark loop interfere destructively. The loop makes the decay rate of the 
two photon channel smaller than those of the other four main channels, for which the Higgs 
boson couples directly to the vector boson or the fermion at the leading order. 



Figure 1-5: Leading order Feynman diagrams for Higgs boson decaying to two photons. 


Even though the Higgs boson decays to two photons at a very small rate, this channel is 
one of the most sensitive channels for the Higgs search in the low rriff region thanks to the 
two isolated high energy photons in the hnal state. The photons—each carrying an energy of 
62.5 GeV for a Higgs boson with mj/ = 125 GeV decaying at the rest—are clearly detected 
and identihed. Their energies and momenta are well measured, from which the diphoton 
mass, rrijy, is reconstructed using the kinematic formula: 

— cos{0^^))^ ( 1 . 1 ) 


where E'^^ and E'^'^ are the measured single photon energies and 9^^ is the measured angle 
between the momenta of the two photons. 

The total decay width of the Higgs boson in the interested mu range is very narrow— 


about 4 MeV at mu = 125 GeV—as shown on the right in Figure 1-4 25,45 . A good 
resolution of both the measured photon energy and the measured open angle therefore leads 
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to a narrow peak of diphoton mass spectrum associated with the Higgs resonance. Because 
the distribution of the background events is expected to be continuously falling, this narrow 
peak provides an eloquent evidence to the existence of the Higgs boson. 

From the amplitude and the location of the peak, the relative total Higgs production 
cross section times the branching ratio to two photons with respect to the Standard Model 
Higgs expectation—the signal strength, and the Higgs mass are measured precisely. The 
rate of the decay, mediated through a loop of particles involving W boson and top quark, is 
sensitive to the magnitudes of both Ky and Kf as well as their relative sign—same sign for 
destructive interference between W loop and top quark loop as expected by the Standard 
Model while opposite sign for constructive interference. It is also sensitive to any possible 
new heavy charged particles in the loop, whose existence is quantihed through measuring 
the effective Higgs coupling strength to photon, Kj 

The dominant background consists of “irreducible” and “reducible” components. The 
“irreducible” component is real (prompt) diphoton events. The “reducible” component in¬ 
cludes dijet and 7 -|- jet events, in which jets are misidentihed as photons. A jet typically 
fakes a photon when it results in a narrow concentration of photonic energy in the detector 
due to the decays of high energy neutral mesons, especially tt^’s. The 7r° decays into two 
photons with small opening angle, which may appear as a single photon. 

Factors for Sensitivity 

The main challenge for the Higgs search through the two photon channel is that the signal 
is much smaller than the background. The expected inclusive signal (S) to background (B) 
ratio S/B under the signal peak, at a Higgs mass of 125 GeV, is about 2% for events at 
8 TeV preselected for the hnal analysis, as evaluated from the numbers in Table 17.11 In order 
to achieve optimal sensitivity of the Higgs search and properties measurements, we need to 
separate the signal and background as much as possible, and further we need to understand 
the background under the signal peak, well. 

The good separation between signal and background depends on the following factors 
related to photons: 

• Good diphoton mass resolution for a narrow diphoton mass peak—requiring good 
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resolution of both single photon energy and diphoton opening angle. 

• Effective separation between prompt photon and a jet faking a photon. 

• Utilization of differences in diphoton kinematics between signal and background. 

In addition, the selections of VBF, VH and itH events according to the features of other 
physics objects produced along with the diphoton, the so called Higgs production tags, are 
another important factors for signal/background separation. Furthermore, these production 
process tags also separate the different signal production processes sensitive to different Higgs 
couplings, which allows the measurement of the signal strengths for individual processes and 
improves the sensitivity of the Higgs coupling strengths. 


Search Strategy 


We use the Compact Muon Solenoid (CMS) detector to detect diphoton events from the 
pp collisions. The CMS detector, with a homogeneous and hne-grained electromagnetic 
calorimeter (ECAL), allows us to identify the photons and to measure their energies with high 
resolution. The photon momentum is obtained using the direction from the reconstructed 
vertex of the associated diphoton production to the photon location in the ECAL, since 
the photon trajectory is not directly measured. The diphoton vertex is selected from all 
the vertices in a bunch crossing, and the efficiency of selecting the correct vertex drives the 
resolution of the diphoton opening angle. The multiple sub-detectors of CMS further allow 
the reconstruction of other physics objects used for the Higgs production tags, including 
electrons, muons, jets and the signature of neutrinos—the imbalance to the total momentum 
projection in the transverse plane with respect to the beam direction (the transverse missing 
energy). 

We design our analysis to maximally separate the signal and background by optimizing 
the diphoton mass resolution for a given phase space and classifying diphoton events accord¬ 
ing to expected S/B under the signal diphoton mass peak. We use Multivariate Analysis 


(MVA) techniques, especially Boosted Decision Trees (BDT) [^ 52 as introduced in Section 
1.3[ to address the key photon factors as follows: 
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Correct the single photon energy and select the diphoton vertex with high efficiency 
to narrow the expected diphoton mass peak for a given phase space. 


• Estimate the energy resolntion of each photon and the probability of selecting the right 
diphoton vertex to bnild a diphoton mass resolntion estimator. 

• Combine all the single photon level information into a photon identihcation classiher 
between prompt photon and fake photon. 

• Combine all the diphoton event level information, including the diphoton mass res¬ 
olution estimator, the photon identihcation classiher for each photon, and diphoton 
kinematics, into a diphoton event classiher which provides a measure of S/B. 

We then use the features of other physics objects produced along with the diphoton to 
select the events into high S/B Higgs production tagged classes, and use the diphoton event 
classiher to select the untagged events into classes with boundaries optimized for the Higgs 
sensitivity. 

We hnally extract the Higgs signal by simultaneous likelihood ht to the diphoton mass 
spectra of all event classes. The expected background under the emerging signal mass peak 
for each event class is constrained directly by the large number of events from data in the 
sidebands of signal region, utilizing the smoothly falling nature of the background shape. 


1.3 Boosted Decision Trees 


Boosted Decision Trees 50 52 is one of the popular MVA techniques, which are employed 


in experimental particle physics to estimate a function mapping a set of input variables 
of an event to its identity as signal or background (classihcation), or to the value of its 
certain property (regression). We choose to use BDT in this analysis for its ability to handle 
large number of input variables and their correlations, as well as its simple mechanism. We 
use BDT to combine all the relevant information in an event into a single variable, which 
maximally separates signal from background for classihcation, or precisely and accurately 
estimates the target property for regression. 
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To construct a classification BDT 52 , or to train a BDT, we provide a signal sample and 


a background sample from Monte Carlo simulated events with known identity, and a selected 
set of input variables with distinguishing power of = {xi,X 2 , ...,Xn}- A single decision tree 
is hrst trained, which is to cut the variable phase space into several signal dominated or 
background dominated hypercube regions, following a certain rule to optimize the separation 
between signal and background, and to label the events in the regions accordingly as “signal” 


or “background”. A demonstration plot of a single decision tree is shown in Figure 1^ It 
has a tree structure, with a root node in magenta representing the entire variable phase 
space, intermediate nodes in yellow representing the split phase spaces, and terminal nodes 
in blue for “signal” regions (SIG) while in red for “background” regions (BKG). The nodes are 
connected by arrows labeled with a variable Xj under consideration and the corresponding cut 
value, which specihes how a parent node is split into two daughter nodes. The tree building 
starts from the root node, with the number of signal and background events reweighted such 
that both the total weights for signal and for background equal to the number of the signal 
events. The node is then split by selecting a single variable and a cut value on it. There are 
several possible splitting criteria. We use the Gini Index dehned as: 


Gini Index = • (1 — ps), 


( 1 . 2 ) 


where ps represents the fraction of the signal weights of the total signal plus background 
weights in a node. The Gini Index is maximal at the root node with ps equal to 0.5. The 
splitting variable and the cut value are chosen to maximize the decrease of the Gini Index 
from the parent node to the two daughter nodes, for which the relative fraction weighted sum 
of the Gini Indices of the two daughter nodes is used. The splitting continues iteratively until 
the predetermined limit is reached, such as the maximum depth of the tree or the minimum 
number of events in a node. The limit is set to decrease the bias due to statistical fluctuation 
of the training samples, the overtraining. The terminal nodes with ps greater (less) than 0.5 
are labeled as “signal” (“background”), and the events in the nodes are assigned a score +1 
(-!)• 


For a single decision tree, some of the events in the terminal nodes are easily misclassihed. 


44 




ROOT 



Figure 1-6: A demonstration plot of a single decision tree. 

and the classification result is susceptible to overtraining. To decrease the misclassihcation 
rate and the effect of overtraining as much as possible, a procedure called “boosting” is 
used, which is basically to train a set of trees and assign a score to an event as the weighted 
average of the scores of all the trees. In our analysis, we employ two boosting procedures. 
Gradient Boost and Adaptive Boost (AdaBoost). The expression for the Gradient Boost is 
as following: 

M 

P e (1.3) 

m=l 

where F{~^;P) represents the function with the set of parameters P corresponding to the 
BDT made up of M trees, represents the function corresponding to the rrith tree, 

am represents the parameters of the rrith tree including the splitting variables and cut values 
at each node, and /3m is the weight on the rrith free. The parameter set P is determined by 
minimizing the deviation between the estimates provided by F{lt ; P) and the true identities 
of the training events, measured by the loss function: 

N 

L(F{^;P),y) = ^ ln(l + F(V;P) e F)}".„ y e (1-4) 

n=l 


45 


where Fn{lt]P) represents the estimated value for the nth event, and i/n is the true value 
+1 or —1 of the nth event, and N is the total number of events. The AdaBoost is obtained 
by minimizing a different type of loss function: 

N 

P),y) = J2 (1.5) 

n=l 


The trained BDT function is then used to assign a score to any event, given its values of 
input variables. The score is a quasi-continuous variable, varying from —1 to +1. The more 
signal-like an event is, the higher value it gets. 

To train a regression BDT |^, we typically provide a sample of Monte Carlo simulated 
events, a target variable corresponding to the desired event property—whose value is known 
for a training event, and a set of other input variables related to the property. The trained 
BDT function provides an evaluation of the property for any event based on its input vari¬ 
ables, which is the weighted average of values estimated by individual decision trees. In the 
case of this analysis, we use the regression in a more generalized way, which regress a prob¬ 
ability density function of the reconstructed energy over the true energy for a photon. We 
provide a known functional form, and set the parameters of the function as target variables. 

For our analysis, we use Toolkit for Multivariate Data Analysis (TMVA) within 


CERN’s ROOT framework 53 to train the classihcation BDTs, while the approach described 


in Reference 39 to train the regression BDT. 
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Chapter 2 


The CMS Experiment at the LHC 


2.1 The Compact Muon Solenoid Detector 


The Compact Muon Solenoid (CMS) detector 28,54 was built to shed light on the mech¬ 


anism of electroweak symmetry breaking by searching for the Higgs boson, to look for de¬ 
viations from the Standard Model by making precise measurements of the Standard Model 
processes, and to search for direct evidence of new physics such as supersymmetry, dark 
matter and extra dimensions. 


An overview of the CMS detector is shown in Figure 2-1 55 . It is a cylindrical detector 


28.7 m long and 15 m in diameter, which is centered at the collision point (LHC point 5) 
with its longitudinal axis along the beam pipe. It is composed of a superconducting solenoid 
magnet and multiple sub-detectors inside and surrounding the magnet. The solenoid provides 
a 3.8 T magnetic held along the longitudinal direction of the detector to bend charged 
particles in the transverse direction. Going from the beam pipe to the solenoid, there is 
a tracker measuring the momenta of charged particles, an electromagnetic calorimeter to 
primarily measure the energies of photons and electrons, and a hadronic calorimeter for 
measuring the energies of charged and neutral hadrons. Outside the solenoid, there are 
muon chambers measuring momenta of muons, which are interleaved with the steel return 
yoke of magnetic hux return. 

To describe the CMS detector we use both right-handed Cartesian coordinates and polar 
coordinates, with the nominal collision point as the origin in both cases. For the Cartesian 
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coordinates, the a:-axis and y-axis are in the transverse plane pointing along the inward radial 
direction of the LHC ring and along the upward vertical direction, respectively, while the 
2 ;-axis is parallel to the beam. For the polar coordinates, 0, r and 6 represent the azimuthal 
angle from the x-axis in the transverse plane, the radial distance in the plane, and the polar 
angle from the z-axis in the y-z plane, respectively. 

The fractions of proton momenta carried by the two colliding partons are generally un¬ 
equal in pp collisions, which leads to the non-zero total collision momentum along the z-axis. 
Under the boost in the z direction, to approximately make an Lorentz-invariant description 
of the hard collision events, with highly relativistic incoming and outgoing particles, the 
pseudorapidity y, dehned as —ln[tan(6'/2)], and the transverse momentum '^t, dehned as 
the projection of the momentum in the transverse plane, are used. The magnitude of the 
transverse momentum is denoted as pt- The presence of high px particles is a signature of 
hard collision events, which is used later for the event selection. 


2.1.1 Tracker 


The tracker [28||56|[57| measures the hit positions of charged particles along their trajectories 
passing through it, which are used to reconstruct the trajectories (tracks), momenta and 
production vertices of the particles. The r-z plane cross section of the tracker is shown in 


Figure 2-2 28 . It is a full-silicon based detector, consisting of an inner silicon pixel detector 


and a silicon strip detector with acceptance \r]\ < 2.5. Silicon is used due to its fast response, 
desired for making measurements from the high luminosity LHC pp collisions, and its good 
spatial resolution. The granularity of the detector decreases with an increase of distance 
from the collision point, which corresponds to a decrease of particle flux. 

The pixel detector has three cylindrical layers in the barrel region at effective radii of r = 
4.4 cm, 7.3 cm and 10.2 cm within \z\ < 26.5 cm, and two disks at \z\ = 34.5 cm, 46.5 cm 
in the endcap on each side within about 6 cm < |r| < 15 cm. It has 66 million pixels each 
with dimensions of 100 pm x 150 pm, which results in an occupancy of about 0.1 permill 
per pixel per bunch crossing. It measures the positions of charge particles hitting its silicon 
wafers with a single point resolution from 15 pm to 20 pm. 

The silicon strip detector consists of inner and outer parts within \z\ < 282 cm and 20 cm 
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MUON CHAMBERS 

Barrel: 250 Drift Tube, 480 Resistive Plate Chambers 
Endcaps: 468 Cathode Strip, 432 Resistive Plate Chambers 


PRESHOWER 

Silicon strips —16m^ —137,000 channels 


FORWARD CALORIMETER 

Steel + Quartz fibres -2,000 Channels 


CRYSTAL 

ELECTROMAGNETIC 
CALORIMETER (ECAL) 
-76,000 scintillating PbW 04 crystals 


HADRON CALORIMETER 

Brass + Plastic scintillator —7,000 channels 


CMS DETECTOR 


Total weight ; 14,000 tonnes 

Overall diameter : 15.0 m 

Overall length : 28.7 m 

Magnetic field : 3.8 T 


SUPERCONDUCTING SOLENOID 

Niobium titanium coil carrying —18,000A 


STEEL RETURN YOKE 

12,500 tonnes 


SILICON TRACKERS 

Pixel (100x150 ^un) -16m^ -66M channels 
Microstrips (80x180 pm) —200m^ —9.6M channels 


Figure 2-1: An overview of Compact Muon Solenoid detector. 


< |r| < 116 cm. The inner part has 4 layers in the Tracker Inner Barrel (TIB), and 3 disks 
in the Tracker Inner Endcap (TIE) on each side. The outer part has 6 layers in the Tracker 
Outer Barrel (TOB), and 9 disks in the Tracker EndCap (TEC) on each side. The whole 
silicon strip detector has 9.3 million strips with thickness of 320 pm or 500 pm and pitches 
from 80 pm to 184 pm. It measures the r-0 or z-(j) positions of charged particles hitting the 
strip detector, with resolutions of 23 pm to 53 pm in the 0 direction. 

The thickness of the tracker material t measured in number of radiation length Xq as a 


function of rj from simulation is shown in Figure 2-3 58 , and has a maximum of about 2. 


The amount of material of the CMS full-silicon based tracker, is much larger than that of 
a tracker utilizing gas detector, e.g. the tracking system of CDF detector at Tevatron has 
a thickness of 0{l% Xq) (^. As a result, the measurement of electron momentum from 
the tracker, and the measurement of electron or photon energy from the electromagnetic 
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Figure 2-2: The cross section of the tracker in r-z plane. 

calorimeter, suffer more from the effects degrading the measurement resolution, including 
multiple scattering, electron bremsstrahlung or photon conversion. The silicon is chosen as 
the tracker material despite of this disadvantage because its fast response and good spatial 
resolution are must for the high luminosity LHC environment. 

For tracks with pt = 100 GeV, the momentum resolution is about 1-2% in the barrel. 
The resolutions for transverse impact parameter dxy and for longitudinal impact parameter 
dz are 0(10 pm). 


2.1.2 The Electromagnetic Calorimeter 


The Electromagnetic Calorimeter (ECAL) 1^,60,61 measures the energies of photons and 


electrons through the electromagnetic (EM) shower they produce traversing the calorimeter. 
An electromagnetic shower for a photon or an electron starts as an electron-positron pair 
production by the impacting photon or the Bremsstrahlung by the impacting electron, and 
develops to a cascade of electrons, positrons and photons through repeating processes of pair 
productions and Bremsstrahlung. The CMS ECAL is designed to measure the photon and 
electron energies with high resolution, which is essential for if —)■ 77 sensitivity. It is homo- 
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Figure 2-3: The thickness of the tracker along with the beam pipe and the support tube, t, 
measured in number of radiation length, Xq, as a function of rj from simulation. 


geneous, fine-grained, and almost hermetic. It is also compact enough to be put inside the 
solenoid to reduce the number of radiation lengths in front, and thus reduce the probability 
of photon conversion and electron Bremsstrahlung before a photon or electron entering the 
ECAL, which improves the resolution of the photon and electron energy measurements. 

It includes a barrel component 


An overview of the ECAL is shown in Figure 2-4 


covering \r]\ < 1.479 and an endcap component on each side covering 1.479 < \t]\ < 3. Both 
barrel and each endcap consist of one layer of lead tungstate (PbW 04 ) crystals, which have 
short radiation length, small Moliere radius, good transparency and fast response as desired. 
Each crystal is coupled to a photodetector: an avalanche photodiode (APD) in the barrel 
and a vacuum phototriode (VPT) in the endcap—which is subject to higher radiation. An 
impacting photon or electron generates an electromagnetic shower through the interaction 
with the crystal and transfers its energy into the shower. The energy of the developed shower 
is then deposited into the crystals, and the crystals emit scintillation light in proportion to the 
deposited energy. The scintillation light in each crystal is converted into photoelectrons and 
amplified by its coupled photodetector, which are further converted into voltages and ADC 
(Analogue-to-Digital Converter) counts. The ADC counts are hnally converted to energy as 
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the measurement of the energy deposit in the crystal, which is later used to reconstruct the 
total energy of the photon or electron. 



Preshower 


Dee 


Modules 


Crystals in a 
supermodule 


Preshower 


Supercrystals 


Endcap crystals 


Figure 2-4: An overview of the electromagnetic calorimeter. 

The ECAL barrel component starts from r = 129 cm. It has 61200 crystals grouped into 
36 supermodules, and each supermodule covers a half barrel in and 20° in 0. Four modules 
are inside each supermodule, and each module has 400 or 500 crystals. Each crystal is in a 
truncated-pyramid shape covering 0.0174 (Ap) x 0.0174 (A0), with a cross section of 22 mm 
X 22 mm in the front end and 26 mm x 26 mm in the back end comparable to the square 
of Moliere radius. The crystal is 23 cm deep, corresponding to 25.8 radiation lengths, which 
well contains the total energy of the electromagnetic shower. The crystal is oriented with 
its axis 3° off from pointing to the nominal collision vertex to avoid particles from passing 
through inter-crystal gaps. 

The ECAL endcaps start from \z\ = 315.4 cm. Each endcap has 7324 crystals grouped 
into two semi-circular parts (“Dees”). Each part consists of 138 units of 5 x 5 crystals 


52 




(supercrystals) and 18 partial supercrystals. Each crystal has a cross section of 28.62 mm 
X 28.62 mm in the front end and 30 mm x 30 mm in the back end. The crystals are 22 cm 
deep, corresponding to 24.7 radiation lengths. The crystals are oriented with their axes 2° to 
8° off from pointing to the nominal collision vertex. To help resolving the two photons from a 
neutral meson, a preshower detector is added in front of each ECAL endcap covering 1.653 < 
{rjl < 2.6. It has two lead disks with thickness of 2 Xq and 1 Xq to generate electromagnetic 
showers. A silicon strip detector with a pitch of 1.9 mm is behind each lead disk to measure 
the shower shape. 


The variation of a crystal response with time and the variation of responses among 


crystals are calibrated and corrected as described References 28,62 . The time variation of 
the response of each crystal is mainly due to changes in crystal transparency from irradiation, 
and subsequent recovery. These variations are monitored by a laser system consisting of lasers 
with wavelength A = 440 nm (near the wavelength of scintillation peak) and A = 796 nm. 
For each crystal, the laser pulse is injected during the beam gap, and a time dependent 
correction factor is computed from the change of the crystal response. The variation of the 
relative responses among crystals are calibrated by a series of methods (intercalibration), 
which use the energy deposition symmetry in 0, the mass of diphotons from 7r° [r]^) decays, 
and the ratio between ECAL energy and tracker momentum of electrons from W and Z 
decays, respectively. A correction factor (intercalibration constant) is obtained for each 
crystal, which is the weighted average of the correction factors from all methods. 

The relative energy resolution cte/E measured in 2006 from the test beam of electrons 
with energy E reconstructed by summing energy deposits in 3 x 3 crystals is 


{aE/EY = {2.8%/VEf + {0.12/Ef + (0.30%)^, (2.1) 


where the hrst term is the stochastic term mainly associated with the shower fluctuation, the 
second term is due to the noise from the readout electronics, the third term is determined 
by the accuracy of calibration. 

The impact positions of photons and electrons in the ECAL are also measured. The 
position resolution in the barrel is 3 mrad in 0 and 0.001 in rj, and the resolution in the 
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endcap is 5 mrad in 0 and 0.002 in r; 62 . 


2.1.3 The Hadronic Calorimeter 

The Hadronic Calorimeter (HCAL) |28p^ measures energies of hadrons through the hadronic 
showers they produce passing through the calorimeter. The HCAL is a sampling calorimeter 
which includes four parts: HCAL Barrel (HB), HCAL Endcaps (HE), HCAL Outer (HO) 
and HCAL Forward (HE). The various sub-components of the HCAL are shown in Figure 


2-5 28 , where a quarter of CMS is portrayed in the r-z cross section. The HB and the HE 


cover respectively |? 7 | < 1.3 and 1.3 < \r]\ < 3, which use brass plates as absorbers to gener¬ 
ate showers and plastic scintillators as the active material to measure the shower energies. 
The HO is added in the central barrel outside the solenoid to increase the HCAL thickness, 
which consists of scintillators and uses the solenoid as absorber. The HE, which uses steel 
as absorber and quartz hbers as active material, covers 3 < < 5.2. The total thickness 

measured in number of nuclear interaction lengths, including the ECAL, ranges from 10 to 
15, depending on rj. 

2.1.4 Muon Detector 


The muon detector 28,64 measures the hit positions of muons along their trajectories 


passing through the detector, which are used to reconstruct their trajectories and momenta. 

, with the various 


A quarter view of CMS in the r-z cross section is shown in Figure 2-6 


sub-components of the muon system labeled. There is a barrel system covering |? 7 | < 1.2 
and endcaps covering 0.9 < \ri\ < 2.4. The barrel is composed of Drift Tube (DT) Chambers 
and Resistive Plate Chambers (RPC). The endcap includes Cathode Strip Chambers (CSC) 
and RPCs as well. The DTs and CSCs are used for precision position measurements, and 
the RPCs are used for fast triggering. 


2.1.5 Trigger 


The trigger system 28,65,66 is used to select potentially interesting physics events to be 


read out and recorded for offline use out of the design rate of 40 MHz pp bunch crossings 
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Figure 2-5: A quarter of CMS r-z cross section. The components of the hadronic calorimeter 
are labeled, including HCAL Barrel (HB), HCAL Endcap (HE), HCAL Outer (HO) and 
HCAL Forward (HF). 

(the actual bunch crossing rate is 20 MHz). The accept rate is subjected to the constraints of 
the detector readout speed, event processing power, and storage space. The trigger consists 
of two levels: the Level-1 (LI) trigger based on hardware and the High-Level Trigger (HLT) 
based on software. The LI uses coarse information from ECAL, HCAL and muon detectors 
to determine approximate candidates of physics objects such as electrons, photons, muons, 
jets and transverse missing energy, and it selects events at a rate of up to 100 kHz. The HLT 
uses detailed information from the entire detector and reconstructs the physics objects in a 
similar way to the offline reconstruction used in the hnal analysis. It further selects events 
with good quality and high pt objects at a rate of 0(100 Hz). 
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Figure 2-6: A quarter view of CMS in the r-z cross section. The components of the muon 
detector are labeled, including Drift Tube (DT) Chambers, Cathode Strip Chambers (CSC) 
and Resistive Plate Chambers (RPC). 

2.2 Event Reconstruction 


2.2.1 Tracks and Vertices 

Tracks, the trajectories of the charged particles in the tracker, propagate as helices between 
tracker layers. They are obtained by htting the tracker hits using the Kalman Filter method 


67 -69 , which takes into account multiple scattering, energy loss and the uncertainty of hit 


positions. A track’s initial momentum, its impact parameter with respect to the nominal 
collision vertex, and its charge result from the £t. Primary vertices are reconstructed by 
grouping tracks compatible with the region of primary interactions according to their distance 


in following a deterministic annealing (DA) algorithm 70 . The position for each vertex is 
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fitted from its corresponding tracks using adaptive vertex fitting 71 . A detailed description 


of track and primary vertex reconstruction is given in Reference 72 


2.2.2 Photons 


A photon, produced at the interaction point, hrst passes through the tracker, and then 
enters ECAL and loses all its energy through electromagnetic shower. There are two cases. 
In the first case, the photon traverses the tracker without interaction and deposits about 
94% (97%) of its energy into 3 x 3 (5 x 5) crystals in the ECAL. Such photon is called 
unconverted photon. In the second case, the photon converts to electron and positron pair 
before entering the ECAL, the electron and positron pair bend under the magnetic held and 
deposit their energies in a larger range in (j). Such photon is called converted photon. To 
include all the photon energy deposits, photons are reconstructed by clustering the energy 


deposits in the ECAL crystals into the so-called superclusters 54,73 


Superclusters in ECAL barrel and those in ECAL endcap are constructed following dif¬ 
ferent algorithms: 


• For ECAL barrel, a seed crystal is hrst located, which is the crystal with highest E-r 
above a certain threshold among the crystals not included in any other supercluster yet. 
Then 5x1 matrices of crystals (bars) each centered at the same rj with the seed crystal 
are built, within the range ±17 crystals in 0 from the seed crystal. The bars with total 
energy above certain threshold connected in cj) are further grouped into clusters called 
basic clusters. The basic clusters with the highest bar energy above certain threshold 
are hnally grouped to form a supercluster. 

• For ECAL endcap, a basic cluster, a 5 x 5 matrix of crystals centered at the seed 
crystal, is hrst built. The crystals at the boundary of the matrix are allowed to seed 
new basic clusters from the crystals not included in any cluster yet. A supercluster is 
then formed from the connected basic clusters. 


The raw photon energy Ejiaw is obtained by summing the energy deposits in the crystals 


of the supercluster calibrated as described in Section 2.1.2 and the energy deposited in the 
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preshower detector is added to it for photons in the endcap. The photon position in rf-cj) 
is obtained from the mean position of basic clusters weighted by energy, and the position 
of basic clusters is calculated from mean positions of crystals corresponding to the shower 


depth weighted by the logarithm of the crystal energy 74 


For a converted photon, if the conversion happens early enough in the tracker such 
that the tracks of the electron and positron pair are well reconstructed, the conversion is 


reconstructed by htting the conversion vertex from the pair of tracks 75 . For the vertex 
htting, the two tracks under consideration are required to have opposite charges, and to be 
parallel at the conversion vertex because the photon is massless, which removes the pollution 
from the random combination of two tracks from the primary interactions (prompt tracks). 
The reconstructed conversion tracks are then each matched to an ECAL supercluster, which 
completes the information about a converted photon. 


2.2.3 Electrons 

Electrons are reconstructed by matching an ECAL supercluster, the same as used for photon 
reconstruction, to a track [^. The candidate track is obtained by htting the tracker hits 
using the Gaussian-sum hlter (GSF) algorithm [^, which models the bremsstrahlung energy 
loss distribution by a weighted sum of Gaussians. 


2.2.4 Muons 


Muons used for this analysis are reconstructed following the so-called global muon recon¬ 


struction method 54,78 , which uses both information from the muon detector and the 


tracker. A muon track only using the muon detector information (standalone muon track) 
is constructed hrst. It starts from building short track traces (segments) from aligned hits 
in individual DT chambers and CSC chambers. These segments are then used for the ht of 
the standalone muon track, following the Kalman Filter method. The obtained standalone 
muon track is matched to a tracker track, and a global muon track is hnally htted from hits 
of both tracks again using Kalman Filter method. 







2.2.5 Jets and Transverse Missing Energy 

Jets and transverse missing energy ME'f’ are reconstructed from electron, muon, photon, 
charged hadron and neutral hadron candidates built from the particle-flow algorithm 

The particle-flow algorithm is designed to reconstruct and distinguish all the stable 
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particles by effectively grouping the information from the entire detector and associating the 
grouped information to each particle candidate, with no information double counted in two 
different candidates. This algorithm provides reconstructed particle candidates (particle-flow 
candidates) as ideal input to reconstruct higher level objects like jets and event level quantity 
like ME^, but not the optimal reconstructed photon. And so we still use the photons from 


more specialized reconstruction algorithm as described in Section |2.2.2| for the diphotons 
candidates in the analysis. 

Jets are built through clustering particle-flow candidates. The anti-fc^ algorithm is 
used and the size parameter is AR = 0.5. Jets from bottom quarks (b-jets) are identihed 
using Combined Secondary Vertex algorithm [^, which identihes the decay vertex displaced 
away from the primary vertex. ME^ is computed as the opposite '^t sum of all the particle- 
flow candidates. 
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Chapter 3 


Higgs Boson to Two Photons Analysis 
Overview 


We perform the analysis to observe the production and decay of the Higgs boson into two 
photons, in the Higgs mass hypotheses range 115 GeV < mn < 135 GeV. The basic flow 
of our analysis is to reconstruct diphotons from the events with at least two reconstructed 
photons, preselect a potentially signal-rich sample of diphoton events, obtain their masses, 
and £t the mass spectrum to search for an excess of signal over background. In the case 
an excess is observed, we further measure its corresponding Higgs mass, and the signal and 
coupling strengths to quantify its compatibility with the Standard Model Higgs boson. To 
optimize the Higgs search sensitivity and measurement precision, we classify the events into 
several categories according to the expected S/B under the mass peak. We use several 
BDTs, trained on the Monte Garlo simulated events, to both improve the diphoton mass 
reconstruction, and to combine all the rest of the diphoton information into a powerful 
diphoton event classifier, which provides a measure of S/B. We also use signatures of Higgs 
production processes to select events into high S/B classes. For the diphoton mass fit, we 
model the signal from Monte Garlo simulated Higgs events, and the background directly 
from the data. 


The main components for this analysis are summarized in Figure 3-1 Gonceptual de¬ 
scriptions of these components are provided in the following of this chapter. The data and 
Monte Garlo simulation samples used for the analysis are introduced afterwards. 
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3.1 Analysis Components 


3.1.1 Diphoton Reconstruction 

To reconstruct diphotons and their masses, we correct the single photon energies and select 
the diphoton production vertex for each diphoton from all the vertices in the same bunch 
crossing. 

Photon Energy Correction 

The raw energy of a reconstructed photon Ejiaw needs correction, as it is deviated from 
the true photon energy Ettus, niainly due to the combined effect of photon shower loss and 
the pileup contamination. The shower lost consists of the part outside of the supercluster 
window, especially for converted photons, and the part passing through the inter-crystal 
gaps or inter-module cracks within the window. The fraction of photon energy lost therefore 
depends on whether it is converted, and the location and detailed pattern of its shower in the 
ECAL. The fraction of energy contaminated depends on the energy density due to the pileup 
interactions in the event. We train a BDT (“photon energy correction regression BDT”) to 
regress the photon energy correction factor, taking the above factors into consideration. 
The target is the probability density of the ratio between the true photon energy and the 
reconstructed raw photon energy ETrue/Euaw, and the input variables are chosen such that all 
the relevant information is included: the supercluster energy, the global detector coordinates 
and local ECAL coordinates of the ECAL clusters, the shower shape variables as measures for 
photon conversion and shower pattern, and pileup information. The trained BDT provides 
an estimation of the probability density of Ettub/E naw for any given photon, and the most 
probable Ettus/E fiaw is used as the correction factor. 

Vertex Selection 

The diphoton production vertex needs to be selected from an average of 9 (21) pp collision 
vertices for 7 TeV (8 TeV) distributed in with an RMS of about 6 cm (5 cm). To keep the 
effect of the vertex selection on the diphoton mass resolution negligible with respect to the 
single photon energy resolution, the selected diphoton vertex is required to be within 1 cm 
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in 2 ; from the true diphoton vertex. For the discrimination between the diphoton production 
vertex and the pileup vertices, we use the knowledge that the total transverse momentum 
of the recoiling tracks, mainly from the underlying events associated with the diphoton 
production vertex, roughly balances the diphoton transverse momentum. The balance is 
not exact as we do not have the association between neutral particles and vertices, so the 
total transverse momentum of neutral particles recoiling against the diphoton for a given 
vertex is unknown. Nevertheless, comparing between the recoiling tracks of the diphoton 
production vertex and those for the pileup vertices, for the former, on average, the sum of 
their transverse momentum square is larger, the relative difference in the magnitude between 
their total transverse momentum and the diphoton transverse momentum is smaller, and the 
projection of their total transverse momentum onto the direction of the diphoton transverse 
momentum is larger. Besides the correlation between the kinematics of the recoiling tracks 
and that of the diphoton, in the case that at least one photon is converted, the position 
of the conversion vertex, together with either the direction of the conversion momentum 
or the position of the ECAL supercluster, provides an extrapolation of the position of the 
diphoton vertex, which is used for the vertex selection. We train a BDT (“vertex selection 
BDT”), using the above information, to distinguish between the prompt vertex and the 
pileup vertices. The BDT assigns scores to the vertices according to how likely it is a 
prompt vertex. The vertex with the highest score is selected. 


3.1.2 Signal to Background Separation 

The reconstructed diphoton events include potential Higgs signal events and a mixture of 
background events. In the background events, there are mainly “irreducible” prompt dipho¬ 
ton events, and “reducible” 7 -|- jet and dijet events with jets faking photons. The fake 
photons are majorly due to energetic neutral mesons, from jet fragmentation, decaying into 
two photons, which end up in the same supercluster and are reconstructed as a single photon. 
The task of the rest of the analysis is to maximize the separation between the signal and 
background. 
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Diphoton Event Preselection 


We first preselect a sample of diphoton events. We design the preselection mainly to select 
the maximum common phase space between the data and Monte Carlo simulated events, 
such that the BDTs trained on the Monte Carlo events are optimal for the data as well, 
the signal model derived from the Higgs Monte Carlo simulation is for the correct phase 
space in data, and the acceptance and efficiency for the Higgs signal is maintained as large 
as possible. We also apply an electron veto to distinguish electrons from photons. 

To select the common phase space, we apply geometric and kinematic acceptance cuts, 
and very loose photon identihcation cuts on the reconstructed photons to remove fake pho¬ 
tons. The photon identihcation depends on two different features between the fake photon 
and the prompt photon. First, the ECAL shower of the fake photon is expected to be wider 
than that of the prompt photon since it is supposed to be the combined shower of the two 
photons. Second, the fake photon is not isolated as other jet fragments leave traces in the 
detector around the photon supercluster. These fragments are reconstructed in the form 
of tracks, energy deposits in the ECAL and HCAL (detector isolation), or the particle-flow 
candidates (particle-flow isolation). We use a set of ECAL shower shape and isolation vari¬ 
ables for the discrimination between the prompt photon and fake photon, and choose the 
corresponding cut values to simulate the effects of the trigger cuts on data, and the generator 
level cuts on Monte Carlo simulated dijet and 7 -|- jet events. This removes most of dijet 
events and a signihcant amount of 7 -|- jet events, while remaining almost fully efficient for 
events with two prompt photons. 

Event Classification 

We then classify the preselected events into classes in the order of roughly S/B under the 
signal mass peak: 

• We hrst select the events into the exclusive tagged classes, based on the signatures of 
the Higgs production processes including vector boson fusion {VBF tag), associated 
production with a W or Z boson (VB tag), and associated production with ti {itH 
tag): 
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— VBF tag: it tags the VBF like events by identifying the additional pair of ener¬ 
getic jets with large separation in 77 . 

— VH tag: it tags the VFI like events by identifying the additional W 01 Z boson in 
its decays to lepton (electron or muon), dijet, or neutrino manifesting as transverse 
missing energy. 

— iiF[ tag: it tags the like events by identifying the additional pair of top quarks 
in their decays to lepton (electron or muon) or multijet. 

• We further classify the untagged events according to their diphoton quality, measured 
by the following elements: 

— Single photon energy resolution: it depends on the same factors as for the photon 
energy correction. 

— Diphoton opening angle resolution: it improves as the probability of selecting the 
right vertex (vertex probability) increases. The vertex probability depends on the 
transverse momentum of the diphoton, the total number of vertices, the number 
of converted photons, and how close the scores of the top ranked vertices and 
their distances to each other. 

— Photon identihcation: the further discrimination, between prompt photons and 
the more photon like fake photons passing the preselection, depends on hner 
photon shower and isolation information, which vary with respect to the energy 
density due to pileup interactions, the photon energy and the photon location. 

— Diphoton kinematics: the two photons from Higgs events have different kinematic 
distributions than those from the background events, because the former are the 
decays from scalar particles, and the initial states for Higgs production are differ¬ 
ent from those of the background events. This provides a way to distinguish the 
Higgs events from the “irreducible” diphoton background. We construct a set of 
variables, which contains full kinematic information of the two photons but with 
diphoton mass factorized as explained below. 

We train BDTs, optimally using the information for individual elements, to build a 
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single photon energy resolution estimator from the width of the probability density of 
ETrue/Enaw (“photou energy correction regression BDT”), a vertex probability estima¬ 
tor (“vertex probability BDT”), and a classiher between the prompt and fake photons 
(“photon identihcation BDT”). We hnally use a BDT (“diphoton BDT”) to construct 
an optimal diphoton event classiher, combining the outputs of all the BDTs for indi¬ 
vidual elements and the diphoton kinematic information. To maximize the expected 
Higgs sensitivity, we classify the events according to the diphoton event classiher. 

The training variables are built such that the diphoton BDT cannot reconstruct the 
diphoton mass to use it distinguishing the signal from the background. This is to 
achieve the same BDT performance for diherent Higgs mass hypotheses since the true 
Higgs mass is unknown. This is also to avoid the preference in selecting background 
events, with diphoton mass close to the Higgs mass of the signal training sample, into 
the high S/B event classes to produce an unwanted peak in the background diphoton 
mass spectrum. There is no loss of sensitivity for this “diphoton mass factorization” in 
the diphoton BDT because the diphoton mass information is used later in the diphoton 
mass £t for the signal extraction. 

3.1.3 Higgs Signal Extraction from Diphoton Mass Fit 

After the event classes are determined, we construct the diphoton mass spectrum for each 

event class, and the corresponding Higgs signal model and background model: 

• The expected diphoton mass spectrum of Higgs signal events is modeled by parametric 
functions, fitted from Monte Carlo simulated events with four Higgs production pro¬ 
cesses mixed according to their cross sections. The discrepancies between data and 
Monte Carlo simulation on photons are evaluated mainly using Z —)■ e’''e“ events with 
the electron reconstructed as the photon. The photons from Z —)■ events are 

used for the validation of Monte Carlo simulation as well, which, though the trans¬ 
verse momentum is on average lower and the statistical uncertainty is larger, provides 
a valuable cross check to the validation using electrons. The Monte Carlo simulation 
related to the vertex selection, which mainly depends on the number of interaction 
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vertices and recoiling tracks from the underlying events for a given vertex, is validated 
using Z —)■ events. The differences between data and Monte Carlo simulation 

are either corrected for or treated as systematic uncertainties for the signal model. 

• The expected diphoton mass spectrum of background events is modeled by parametric 
functions with a smoothly falling feature, htted directly from data. The background 
model in the signal region for any Higgs mass hypothesis under consideration is con¬ 
strained by the background events in the sidebands. The htting range is set as 100 GeV 
< < 180 GeV, to get the signal region well contained, and to get sufficient number 

of background events in the sidebands. The uncertainty on the Higgs signal extrac¬ 
tion and measurements due to the limited knowledge of the exact background shape 
is evaluated by prohling over a set of functions well describing the data and general 
enough to cover the true background function. 

The Higgs signal is hnally extracted by statistical procedures based on simultaneous 
likelihood £t to the diphoton mass spectra over all event classes. 
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Figure 3-1: Higgs boson to two photons analysis workflow. The blue circles represent the 
input elements to the analysis from the event reconstruction. The green boxes represent 
BDTs used for information processing. The yellow circles represent the quantities built from 
input information. The red box represents the process for Higgs signal extraction. 

3.2 Data and Monte Carlo Simulation Samples 

We analyze the full datasets collected by the CMS detector in 2011 and 2012 LHC run 
periods. The 2011 and 2012 datasets consist of pp collision events respectively at center of 
mass energy a/s = 7 TeV with a integrated luminosity L = 5.1 fb“^, and at ^/s = 8 TeV 
with L = 19.7 fb“^. An event only gets selected if it passes either of the following two classes 
of diphoton High-Level Triggers designed for 77 77 : 

• Trigger 1: the photon energy projected in the transverse plane > 26 GeV for the 
photon with the highest (leading photon), E^ > 18 GeV for the photon with the 
second highest Et (sub-leading photon), and both photons passing Level-1 trigger. 




























• Trigger 2: > 36 GeV for the leading photon, > 22 GeV for the sub-leading 

photon, and at least one photon passing Level-1 trigger. 

For both types of trigger, the leading and sub-leading photons are required to pass loose 
photon identihcation requirements based on shower shape and isolation. The trigger effi¬ 
ciency is 99.4% for events selected for the hnal statistical analysis, evaluated using the “Tag 


and Probe” method 83 on Z —)■ e’''e events. 


We use Monte Garlo simulation samples of iL —77 to train the BDTs, optimize the 
event classihcation, and build the signal model of the diphoton mass distribution. The 
Ff —)■ 77 samples are produced for all the four production processes ggH, VBF, VH and iiH 
at Higgs mass hypotheses ranging from 115 GeV to 135 GeV, at both ^/s = 7 TeV and a/s = 
8 TeV. To decrease the effect of the statistical fluctuation of any particular sample, samples 
at different Higgs masses are used in general for the BDT trainings, event class optimization 
and signal modeling, respectively. For ggH and VBF processes, POWHEG 84 88 is used 


for matrix element generation at next-to-leading order (NLO), and PYTHIA 89 is used for 
parton showering and hadronization. For VH and iiH processes, PYTHIA is used for both 
matrix element generation at leading order (LO), and parton showering and hadronization. 
The production cross sections for the Standard Model Higgs boson, and the branching ratio 
for its decay to two photons that are used are from the LHG Higgs boson Gross Section 


Working Group 45 . To describe the Higgs kinematics, we match the distribution of the 


transverse momentum of the Higgs boson from ggH process to the next-to-next-to-leading 


logarithmic resummation (NNLL) plus NLO calculations from HqT 90-92 , by reweighting 
the produced events at 7 TeV, and tuning POWHEG for event generation at 8 TeV according 


to Reference 93 respectively. To account for the effect of interference between the ggH 


process and the continuum gg —)■ 77 process, we reduce the the cross section for ggH process 


by 2.5% 94 


We use Monte Garlo simulation samples of background processes to train the BDTs 
and to optimize the event classihcation. For the “irreducible” diphoton background at 


7 TeV, the sample of diphoton Born process is generated using MADGRAPH 95 inter 


faced with PYTHIA, and the sample of Box process is generated using PYTHIA. The “ir¬ 
reducible” diphoton background at 8 TeV both Born and Box processes are generated using 
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SHERPA 96 , which provides a better description of the events in the phase space with ad¬ 
ditional jets from Initial State Radiation(ISR). For the “reducible” background, the samples 
of 7 -I- jet process and dijet process are generated using PYTHIA. A “double EM-enriched 
hlter” including loose isolation cuts is applied to select the events which are likely to pass 
the later diphoton selection of the analysis, in order to save computing power for the further 
simulation of interactions between the particles and the detector. The background cross 
sections are calculated at LO and corrected by a scale factor from 1.0 to 1.3 obtained from 


CMS measurements 97,98 


For both signal and background Monte Carlo samples pileup interactions are simulated 
using PYTHIA. For event class optimization and signal modeling, the Monte Carlo events 
are reweighted to match the pileup distribution in data. The detector response is simulated 
using GEANT4 |^. The discrepancy between data and Monte Carlo simulation is evaluated 
using events from Z —>■ Z —/i+p“ 7 , and Z —)■ data and Monte Carlo simulation 

generated using POWHEG. Comparison between data and Monte Carlo distributions of the 
number of reconstructed vertices in Z —events after pileup reweighting for 7 TeV 


and 8 TeV are shown on the left and right of Figure 3-2 39 . Good agreement is observed. 
The data and Monte Carlo samples including the additional background samples used for 


the VH and ttH tags are listed in detail in the analysis note of Reference 39 
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Vs = 7 TeV L = 5.1 fb 1 VS = 8 TeV L = 19.7 fb ^ 



Figure 3-2: The comparison between data and Monte Carlo distributions of the number of 
reconstructed vertices in Z —>■ events after pileup reweighting for 7 TeV (left) and 

8 TeV (right) are shown. 
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Chapter 4 


Diphoton Reconstruction and 
Selection 


We preselect a potentially signal-rich sample of diphoton events as described in Section |4.1 
The diphotons are reconstructed with corrected photon energy and selected diphoton vertex. 
The photon energy correction and the diphoton vertex selection are described in Section |4.2| 


and Section 4.3.1, respectively. To further classify the diphoton events according to S/B 


under the diphoton mass peak, individual BDTs are hrst trained to provide a single photon 
energy resolution estimator as described in Section 14.2. a diphoton vertex probability esti¬ 
mator as described in Section 14.3.21 and a single photon identihcation classiher as described 


in Section 4.4 A diphoton BDT is then trained to combine the outputs of the above BDTs 


into a single diphoton event classiher as described in Section |4.5[ which provides a measure 
of expected S/B for each diphoton event and is used for event classihcation later. 


4.1 Diphoton Event Preselection 

For each event with at least two reconstructed photons, the diphoton pairs are hrst recon¬ 
structed by grouping the reconstructed photons into all possible two photon combinations. 
For each diphoton pair, a primary vertex is selected as described in Section 4.3.1 The 


momentum of each photon is constructed with its magnitude obtained from the corrected 


photon energy as described in Section 4.2 and its direction pointing from the selected ver- 
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tex to the supercluster. A preselection is then applied to the diphotons, which consists of 
a so called single photon preselection on each photon, and a set of cuts on the diphoton 
kinematic acceptance. If more than one diphoton pair pass the preselection, the diphoton 
with the maximum scalar sum of photon transverse momentum is used for the analysis. A 
detailed description of the preselection is given below. 

4.1.1 Single Photon Preselection 

The single photon preselection includes a cut on the acceptance of supercluster pseudorapid¬ 
ity measured with respect to the origin of the detector coordinate rjsc, a set of loose photon 
identihcation cuts against jets faking photons and an electron veto. 

Acceptance of Supercluster Pseudorapidity 

The acceptance on the supercluster pseudorapidity is determined to exclude the transition 
region between the ECAL barrel and endcap, and the region outside the tracker acceptance, 
which is < 1-4442 in the barrel or 1.566 < \f]sc\ < 2.5 in the endcap. 

Loose Photon Identification Cuts 

The loose photon identihcation cuts are applied to a set of ECAL shower shape and isolation 
variables dehned as following: 

• Shower shape variables 

— Rg: the ratio between the energy in the 3x3 crystals centered at the seed crystal 
and supercluster energy. 

— airjirj'. the log-energy weighted standard deviation of single crystal t] in crystal 
index within the 5x5 crystals centered at the seed crystal. The weight per- 
crystal is 4.7 plus the logarithm of the ratio between the energy in the crystal to 
the energy in the 5x5 crystals. If the weight is negative then 0 is used instead. 

• AR\ the separation in the rj-cp plane a/ At]'^ + A(j?. 

• Detector isolation variables 
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— H/E: the ratio between the sum of energies of deposits in HCAL within AR < 0.15 
from the ECAL supercluster, and the ECAL supercluster energy. 

— ISOTrk'- the sum of pt of tracks within 0.04 < AR < 0.3 from the photon mo¬ 
mentum direction. The photon momentum direction used in this case is obtained 
with respect to the vertex with the maximum sum of track and only tracks 
matching this vertex are included in the isolation computation. 

— ISOTrkPtCorr- ~ 0.002p'r 

— ISO/zcal: the scalar sum of transverse energies of deposits in HCAL within 0.15 
< AR < 0.3 from the ECAL supercluster. 

— ISOHCALPtCorr- ^SOhCAL “ 0.005pT 


• Particle-flow isolation variable 
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— ISO PFChargedSeivtx02- the sum of Pt of particle-flow charged hadron within 0.02 < 
AR < 0.2 from the photon momentum direction. Only the particle-flow charged 
hadrons with impact parameter along z direction \dz\ < 0.2 cm and transverse 
impact parameter \dxy\ <0.1 cm with respect to the selected photon vertex are 
included for the isolation computation. 


To apply the loose photon identihcation cuts, the photons are classihed into four categories 
according to the photon supercluster location in the ECAL (barrel or endcap) and the value 
of Rg (> 0.9 or < 0.9). The photons in the barrel and endcap are treated separately because 
the geometry of the crystals and the amount of tracker materials in front are different for 
ECAL barrel and endcap. The value of Rg is used as a measure of the shower width, and 
the photons with higher Rg are more likely to be prompt photons. The cut values are in 
Table lO 


Electron Veto 

The electron veto is used to distinguish electrons from photons. The photon candidates 
having the same supercluster with a GSF electron candidate are removed. To avoid rejecting 
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Table 4.1: The loose photon identihcation cuts for single photon preselection. The photons 
are divided into four categories according to the photon supercluster location in the ECAL 
(barrel or endcap) and the value of Rg (> 0.9 or < 0.9). The cut values vary with the photon 
categories. 


Rg > 0.9 

Barrel 

Endcap 

H/E 

< 0.082 

< 0.075 


< 0.014 

< 0.034 

ISOnCALPtCorr 

< 50 GeV 

< 50 GeV 


< 50 GeV 

< 50 GeV 

ISO pFChargedSelVtx02 

< 4 GeV 

< 4 GeV 

Rg < 0.9 

Barrel 

Endcap 

H/E 

< 0.075 

< 0.075 


< 0.014 

< 0.034 

ISOnCALPtCorr 

< 4 GeV 

< 4 GeV 

^SOprkPtCorr 

< 4 GeV 

< 4 GeV 

ISO pFChargedSelVtx02 

< 4 GeV 

< 4 GeV 


the converted photons, the electron track is required to have no missing hits in the tracker 
before its hrst hit, and not to match an identihed conversion. 


4.1.2 Diphoton Kinematic Acceptance 

The cuts of diphoton kinematic acceptance are determined to select the phase space right 
above the trigger threshold and to dehne a region for the diphoton mass £t. The cuts 
include > 1/3 and >1/4, for leading photon 7I and sub-leading photon 

72 respectively, and 100 GeV < < 180 GeV. The threshold of the transverse momentum 

for photons entering the analysis is thus 100 GeV/4 = 25 GeV. 


4.1.3 Selection Efficiencies and Scale Factors Between Data and 
Monte Carlo Simulation 


The efficiencies of the loose photon identihcation cuts for the prompt photons in the four 
photon categories are evaluated using electrons from Z —)■ e’''e“ events, for which the electron 

is 


Rg is rescaled to match the photon Rg distribution. The “Tag and Probe” method 83 


used to evaluate the efficiencies on data and Monte Garlo simulation at 7 TeV and 8 TeV 
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respectively. The efficiencies as well as the corresponding efficiency scale factors, ratios 


between efficiencies on data and Monte Carlo simulation, are in Table 4.2 


Table 4.2: The loose photon identihcation efficiencies for prompt photons from data and 
Monte Carlo simulation at 7 TeV and at 8 TeV as well as the corresponding efficiency scale 
factors between data and Monte Carlo simulation. The photons are classihed into four 
categories according to the photon supercluster location in the ECAL (barrel or endcap) 
and the value of Rg (> 0.9 or < 0.9). The efficiencies are evaluated using electrons from 
Z —)■ events. 


ys = 7 TeV 

Data 

Monte Carlo 

Data/Monte Carlo Scale Factor 

> 0.9 Barrel 
i ?9 < 0.9 Barrel 
i ?9 > 0.9 Endcap 
i ?9 < 0.9 Endcap 

0.9872 ± 0.0025 

0.9619 ± 0.0050 

0.9906 ± 0.0085 

0.9606 ± 0.0150 

0.9908 ± 0.0002 

0.9670 ± 0.0005 
0.9824 ± 0.0004 

0.9560 ± 0.0011 

0.996 ± 0.003 

0.995 ± 0.006 

1.008 ± 0.009 

1.005 ± 0.018 

ys = 8 TeV 

Data 

Monte Carlo 

Data/Monte Carlo Scale Factor 

i ?9 > 0.9 Barrel 
i ?9 < 0.9 Barrel 
i ?9 > 0.9 Endcap 
i ?9 < 0.9 Endcap 

0.9879 ± 0.0030 

0.9566 ± 0.0055 

0.9838 ± 0.0090 

0.9545 ± 0.0170 

0.9864 ± 0.0001 

0.9610 ± 0.0002 

0.9789 ± 0.0002 
0.9445 ± 0.0005 

0.999 ± 0.003 

0.995 ± 0.006 

1.005 ± 0.009 

1.011 ± 0.018 


The electron veto efficiencies for the prompt photons are evaluated using photons from 
Z —)■ /r+/i “7 events. The photons are classified into four categories according to the photon 
supercluster location in the ECAL (barrel or endcap) and the value of i ?9 (> 0.94 or < 0.94). 
The value of Rg is used as a measure of the likelihood of photon conversion. The photons 
with i ?9 > 0.94 are dominated by unconverted photons while photons with R^ < 0.94 are 
dominated by converted photons. The efficiencies on data and Monte Carlo simulation at 
8 TeV as well as the corresponding efficiency scale factors between data and Monte Carlo 
simulation are shown in Table 14.31 The efficiencies on data and Monte Carlo simulation at 
7 TeV and the corresponding scale factors come out to be 1. 
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Table 4.3: The electron veto efficiencies for prompt photons from data and Monte Carlo at 
8 TeV as well as the corresponding efficiency scale factors between data and Monte Carlo. 
The photons are classihed into fonr categories according to the photon snpercluster location 
in the ECAL (barrel or endcap) and the value of (> 0.94 or < 0.94). The efficiencies are 
evaluated using photons from Z —)■ /i+/r “7 events. 


ys = 8 TeV 

Data 

Monte Garlo 

Data/Monte Garlo Scale Factor 

Rg > 0.94 Barrel 

Rg < 0.94 Barrel 

Rg > 0.94 Endcap 
Rg < 0.94 Endcap 

0.9984 ± 0.0003 
0.9867 ± 0.0012 

0.9893 ± 0.0016 

0.9639 ± 0.0033 

0.9991 ± 0.0003 

0.9930 ± 0.0009 

0.9938 ± 0.0012 

0.9738 ± 0.0030 

0.9994 ± 0.0004 

0.9937 ± 0.0014 

0.9955 ± 0.0020 

0.9899 ± 0.0045 


4.2 Photon Energy Correction 


4.2.1 Photon Energy Correction Regression BDT 


The photon energy correction regression BDT is trained to provide each photon 
factor to its raw energy, and a per-photon energy resolution estimator, which is 


diphoton BDT as described in Section 4.5 


a correction 
used for the 


Training Samples 

The training sample for the BDT is composed of reconstructed photons from Monte Carlo 
7 + jet events. Each photon is required to match a prompt photon at the generator level, 
and the generated energy of the prompt photon is used as the true photon energy ETme- In 
addition, the photon is required to pass the single photon preselection with px > 15 GeV, 
looser than the analysis threshold of 25 GeV to increase the size of training sample. The 
trainings are performed separately for photons from pp collisions with different center-of- 
mass energies (7 TeV or 8 TeV) and in different EGAL locations (barrel or endcap). 


Input Variables 

The input variables are summarized as following: 
• Snpercluster variables: 
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— Esc'- energy deposit in the ECAL supercluster. 

“ Vsc- pseudorapidity of ECAL supercluster measured with respect to the origin of 
the detector coordinate. 

— Rq: the ratio between the energy in the 3x3 crystals centered at the seed crystal 
and supercluster energy. 

— H/E: the ratio between the sum of energies of deposits in HCAL within AR < 0.15 
from the ECAL supercluster, and the ECAL supercluster energy. 

— SC ? 7 -Width: the energy-weighted standard deviation of single crystal eta in de¬ 
tector coordinate within supercluster. The weight per-crystal is the ratio of the 
single crystal energy to the supercluster energy. 

— SC 0-Width: the energy-weighted standard deviation of single crystal phi in de¬ 
tector coordinate within supercluster. The weight per-crystal is the ratio of the 
single crystal energy to the supercluster energy. 

— The number of basic clusters. 

— The supercluster azimuthal angle 4>sc- (This is only used for the barrel since 
its inclusion does not improve the resolution for electrons in the endcap from 
Z —)■ events in data.) 

— Ratio between preshower energy and supercluster energy (endcap only). 

• Seed basic cluster variables: 

— Ratio between seed basic cluster energy and supercluster energy. 

— Seed basic cluster rj and 0 relative to the supercluster 

— airfir^: the log-energy weighted standard deviation of single crystal t] in crystal 
index within the 5x5 crystals centered at the seed crystal. The weight per- 
crystal is 4.7 plus the logarithm of the ratio between the energy in the crystal to 
the energy in the 5x5 crystals. If the weight is negative then 0 is used instead. 

“ o'itjji,!,'- the log-energy weighted standard deviation of single crystal 0 in crystal 
index within the 5x5 crystals centered at the seed crystal. 
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— coYirficj)'. the log-energy weighted covariance of single crystal rf-cj) in crystal index 
within the 5x5 crystals centered at the seed crystal. 

— Ratios between energies of various combinations of crystals within the seed basic 
cluster and seed basic cluster energy. 

• Seed crystal variables: 

— Seed crystal rj and 0 relative to the seed basic cluster. 

• Pileup variables: 

— pEvent- the estimate of transverse energy per unit area in the rj-cj) plane contributed 
by the pileup interactions and underlying-event effects in the event. It is the 


— Nytx- the number of reconstructed vertices. 


median of jets constructed using the kx algorithm 


inn 


Output and Performance 

The target is the probability density of ETme/ERaw for any photon with input variable . 
It is parametrized empirically using a modihed Crystal Ball function (CBdoubie-sided) 
consisting of a Gaussian core and power law tails on both sides: 


Target = CBdoubie-sided| , cr(^), 0 ^(^), (^), , nR(^)), (4.1) 


where piF^) and are the mean and standard deviation of Gaussian core, and aL{R)iF^) 

and nE{R)(F^) are the cut off and power of left (right) tail. The parameters are functions of 
the input variables iF estimated by BDT and are determined by maximum likelihood £t. 

The trained BDT estimates the probability density of ETrue/Enaw for each photon ac¬ 
cording to its input variable The performance of the estimation is evaluated on a testing 
Monte Carlo sample of photons, independent from the training sample. As shown on the 


left (right) in Figure 4-1 [^, for photons in the barrel (endcap), the normalized sum of 
the estimated ETme!ERaw distribution for each photon (blue line) agrees well with the true 
Ettuc/ER aw distribution of the sample (points). 
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Figure 4-1: The normalized sum of the individual photon ETme/E^aw distributions estimated 
by the regression BDT (blue line), compared to the true ExmelEuaw distribution (points), 
for photons in the barrel (left) and in the endcap (right) of a Monte Carlo sample independent 
from the training sample. 


For each photon, its energy is corrected to the most probable value of the true energy 
E^E^aw, '^) by multiplying the correction factor as: 

E{ERaw, = fi{'^)E Raw (^*^) 

The per-photon energy resolution estimator [aE/E)(^) is assigned as: 

(crs/F;)(^)=a(^)/M^). (4.3) 

4.2.2 Energy Correction Between Data and Monte Carlo Simula¬ 
tion 

The imperfect simulation of detector effects causes discrepancies in the scale and resolution 
of regression photon energy between data and Monte Carlo simulation. The discrepancies 
are corrected for building the model of diphoton mass spectrum for the Higgs boson from 
Monte Carlo simulation, which is used in the signal extraction and is crucial for the Higgs 
mass measurement. The corrections are derived from Z —)■ e+e“ events from data and Monte 
Carlo simulation with the electron ECAL supercluster reconstructed the same way as the 
photon supercluster, and performed in a three-step procedure. The hrst step corrects the 
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energy scale difference mainly due to the imperfect correction for the crystal transparency 
loss in data, which varies with time and photon location. The second step mainly corrects 
the underestimation of the energy resolution in Monte Carlo, which varies with the photon 
location and whether a photon is converted. The third step corrects the residual energy scale 
difference as a function of photon energy for photons from 8 TeV data, for which enough 
statistics is available for the derivation of this hne-grained correction. A detailed description 
of these three steps are as follows. 

First, the energies of photons from data are scaled to match the energy scale of Monte 
Carlo simulated photons. The scale factors are derived separately for photons from different 
LHC run ranges and located in different pseudorapidities, which are classihed into 59 run 
ranges x 4 \risc\ ranges, 2 \r]sc\ ranges for barrel and 2 \risc\ ranges for endcap, single photon 
categories. To derive the scale factor Rgtepi for each category, the mass spectra are built for 
Z —)■ e’''e“ events from both data and Monte Carlo simulation, with both electrons from 
the same single photon category. Each spectrum is htted by an expected mass distribution 
p(mee), parametrized as a Breit-Wigner (BW) function convoluted with a Crystal Ball (CB) 
function: 

p(mee) = BW(mee | mz,Tz) X CB(mee I AM,Aa,a,n), (4.4) 


where the Breit-Wigner function models the intrinsic distribution of Z —>■ e’''e , with the 


peak mass mz and width T^ parameters hxed to the Particle Data Group values 102 . The 
Crystal Ball function models effects from the detector measurement with the parameters 
mean and standard deviation of Gaussian core AM and Aa, and cut off and power of the 
power law tail a and n, floating during the ht. The scale factor Rgtepi is then obtained as the 
ratio between the measured mass mpeak mc of Monte Carlo (MC) simulation and mpeak Data 
of data as: 


R 


stepl 


'^peak MC 
^^peak Data 


mz + AMmc 
mz + AM Data ’ 


(4.5) 


where AMMc(Data) is the htted mean of the Gaussian core of the Crystal Ball function for 
Monte Carlo simulation (data). 


Second, the energies of photons from Monte Carlo simulation are smeared to match the 
energy resolution of photons from data, while the energies of photons from data are further 
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scaled to correct the residual scale difference with Monte Carlo simulation. The corrections 
are derived for 2 i ?9 x 4 \risc\ single photon categories, with a smearing parameter a smear 
and a scale factor Rstep 2 for each category. To derive the corrections, the double electrons 
from data and Monte Carlo simulation are classified into 36 categories according to the 
single photon categories of the two electrons. For each double electron category, the energy 
of each Monte Carlo electron is scaled by a random factor from a Gaussian distribution with 
mean at 1 and standard deviation agmear corresponding to its single photon category. The 
histogram of the smeared double electron mass rriee for Monte Carlo events is constructed 
correspondingly, which is a function of the smearing values of both electrons 
where i and j represent the related single photon category numbers. The energy of each 
electron from data is scaled by Rstep 2 corresponding to its single photon category, and the 
scaled double electron mass spectrum is built as a function of {Rltep 2 j Ritep 2 )- The smeared 
Monte Carlo histogram of rugg is then fitted to the scaled data. The 8 pairs of {asmear, Rstep 2 ) 
for each single photon category are determined by maximizing the total likelihood of the 36 
double electron categories. 

The smearing parameter Cgmear is parametrized as a constant for electrons in the endcap 
at 7 TeV and 8 TeV. An additional energy dependent term is added in quadrature for 
electrons in the barrel at 8 TeV, where more statistics are available to fit the improved 
parametrization: 


(Js 


Cl (7 TeV and 8 TeV Endcap), 

yCf + {C,Js/BtY (8 TeV Barrel), 


(4.6) 


where Ci and C 2 are constants for each single photon category, and Et is the photon trans¬ 
verse energy. 

Third, a residual energy dependent scale factor is applied to the energies of photons in the 
barrel from data at 8 TeV. The scale factors are derived for 20 R^ x \risc\ x Et categories 
following the same method as in the second step. 

The comparison between 8 TeV data and Monte Carlo simulated Z —>■ e+e“ mass dis¬ 


tributions after energy corrections are shown in Figure 4-2 40 . The events are requested 


to pass the preselection with inverted electron veto. The distributions for the events with 
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both electrons in the barrel are shown on the left, and the distribntions for the events with 
at least one electron in the endcap are shown on the right. Good agreement is observed and 
the residnal difference is taken into acconnt as the systematic nncertainties dne to correction 


method in the signal modeling for the hnal statistical analysis as described in Section 


7.3 




rriee (GeV) 


rriee (GeV) 


Figure 4-2: The Z —)■ e’''e“ mass distributions of data (points) and Monte Carlo simulated 
events (histogram) at 8 TeV after energy corrections with both electrons in the barrel (left) 
and at least one electron in the endcap (right). The electron ECAL superclusters are recon¬ 
structed in the same way as the photon superclusters, and the events are requested to pass 
the preselection with inverted electron veto. 


The relative resolution estimator (Je/E for each single photon is smeared as well for both 
data and Monte Carlo simulation, by adding in quadrature the smearing parameter agmear 
for the corresponding photon category. The discrepancy between the (Te/E distributions of 
Monte Carlo simulation and data due to the imperfect simulation of detector response are 
evaluated using electrons from Z —)■ events and photons from Z —)■ events. A 

scaling of ±10% of the Monte Carlo ue/E is shown to cover the discrepancy. 
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4.3 Vertex Selection 


4.3.1 Vertex Selection BDT 

The vertex selection BDT is trained to select the diphoton production vertex in an event, 
which has an average of 21 (9) pp collision vertices for 8 TeV (7 TeV) as a result of pileup 
interactions. 

Training Samples 

The training is performed on a Monte Carlo simulation of if —t 77 events. The signal sample 
consists of the reconstructed vertices of diphotons from Higgs decays, while the background 
sample consists of the pileup vertices. 

Input Variables 

The input variables are the following: 

• of square of transverse momentum of each track associated 

with the vertex, This quantity is expected to be larger for the diphoton vertex 

than for pileup vertices. 

• ■ {~^'t IPt)'- projection of the sum of transverse momenta of tracks 

associated with the vertex onto the opposite direction of the diphoton trans¬ 

verse momentum . This quantity is expected to be near 0 for the pileup vertices 
while near for the diphoton vertex. 

• (I Xli ~ Pt^)/(I Sj + Pt)'- asymmetry between the magnitude of the 

vector sum of transverse momenta of tracks associated with the vertex, | 

and the magnitude of diphoton transverse momentum, . This quantity is expected 
to be near —1 for the pileup vertices while near 0 for the diphoton vertex. 

• \zvtx — zconv\lO'zconv cases with at least one converted photon: the distance 

between the 2 ; position of the vertex, zytx) and the estimated position of diphoton 
vertex, zconv, from the conversion and normalized by the uncertainty of the estimation. 
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(Tzconv- This quantity is expected to be near 0 for the diphoton vertex while larger for 
pileup vertices. 

Output 

The BDT output is a score assigned to each vertex which ranges from —1 to 1. The higher 
the score assigned to a vertex, the more likely the vertex is the diphoton production vertex. 
The vertex with the highest BDT score is selected as the diphoton vertex. 


4.3.2 Vertex Probability BDT 

The vertex probability BDT is trained to estimate the probability that the selected vertex is 
the correct diphoton vertex for each event. The criteria for being correct is that the distance 
between the selected vertex and the true diphoton vertex is within 1 cm in the z direction, 
in which case the diphoton mass resolution is insensitive to the exact position of the vertex. 
The vertex probability is a measure of the diphoton opening angle resolution. It is used for 


the diphoton BDT as described in Section 4.5 


Training Samples 

The training is performed on a Monte Carlo simulation of if —t 77 events. The signal sample 
consists of events with correct vertex selected, and the background sample consists of events 
with wrong vertex selected. 


Input Variables 

The input variables are the following: 

• Nvtx'- the number of reconstructed vertices. 

• piP'. the magnitude of diphoton transverse momentum. 

• The top three vertex selection BDT scores for the vertices in the event. 

• The distances in z between the selected vertex and the vertices with the second and 


the third highest BDT scores. 



• The number of conversions in the diphoton (0, 1 or 2). 

Output 

The BDT output is a score assigned to each event, which ranges from —1 to 1. Events 
are binned according to BDT score, and the diphoton vertex selection efficiency in each 
bin, dehned as the fraction of events with diphoton vertex selected correctly in the bin, is 
measured. A linear relation between the vertex selection efficiency and the BDT score is 
derived, which is used to transform a BDT score to a per-event vertex probability between 
0 and 1 . 


4.3.3 Performance 


To measure the performance of both the vertex selection BDT and the vertex probability 
BDT, the diphoton vertex selection efficiency and the average vertex probability are eval¬ 
uated on Monte Carlo simulated if —)■ 77 events at a Higgs mass of 125 GeV, in bins of 


p'p. As shown in Figure 4-3 40 , the average vertex probability along with uncertainty 


(blue band) predicts well the measured vertex selection efficiency (data points), and both 
increase with the increasing . The total vertex selection efficiency is 79.6% (85.4%) for 
the if —77 events at a Higgs mass of 125 GeV at 8 TeV (7 TeV). The efficiency at 8 TeV 
is lower than that at 7 TeV because of higher number of pilup interactions. 
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8TeV 



Figure 4-3: The measured vertex selection efficiency (points) and the average vertex proba¬ 
bility along with uncertainty (blue band) evaluated from the BDT on Monte Carlo simulated 
Ff —)■ 77 events at a Higgs mass of 125 GeV at 8 TeV, in bins of pp. 


4.4 Photon Identification BDT 


The photon identification BDT is trained to provide each photon with a score, measuring 


how likely it is a prompt photon rather than a jet faking a photon (fake photon), which is 
used as an input to the diphoton BDT as described in Section 4.5[ 


4.4.1 Training Samples 

The training is performed on a Monte Carlo simulation of 7 -|- jet events passing the pre¬ 
selection with pt > 15 GeV, looser than the analysis threshold of 25 GeV to increase the 
size of the training sample. The training is tested with another training on a sample passing 
the preselection with px > 25 GeV. The performances of BDTs from both trainings agree 
well. The signal sample consists of the reconstructed photons which match prompt pho¬ 
tons at the generator level, while the background sample consists of the ones that do not 
match. The trainings are performed separately for photons from pp collisions with different 
center-of-mass energies (7 TeV or 8 TeV) and in different ECAL locations (barrel or endcap). 
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4.4.2 Input Variables 

The input variables for the photon identification BDT are listed as following: 

• Shower shape variables: 

— /? 9 : the ratio between the energy in the 3x3 crystals centered at the seed crystal 
and the supercluster energy. 

— SC //-Width: the energy-weighted standard deviation of single crystal // in detector 
coordinate within supercluster. The weight per-crystal is the ratio of the single 
crystal energy to the supercluster energy. 

— SC 0-Width: the energy-weighted standard deviation of single crystal 0 in de¬ 
tector coordinate within supercluster. The weight per-crystal is the ratio of the 
single crystal energy to the supercluster energy. 

— aijjirj- the log-energy weighted standard deviation of single crystal rj in crystal 
index within the 5x5 crystals centered at the seed crystal. The weight per- 
crystal is 4.7 plus the logarithm of the ratio between the energy in the crystal to 
the energy in the 5x5 crystals. If the weight is negative then 0 is used instead. 

— covirjif the log-energy weighted covariance of single crystal rj-cf) in crystal index 
within the 5x5 crystals centered at the seed crystal. 

“ E 2 X 2 /E 5 x 5 - the ratio of the energy in the 2x2 crystal array containing the 
seed crystal (the 2x2 crystal array with the highest energy in all the possible 
combinations) to the energy in the 5x5 crystals centered at the crystal. 

— Preshower (endcap only): the sum in quadrature of the energy-weighted 
standard deviation of the strip index in the x and y planes of the preshower 
detector. 


Particle-flow based isolation variables 79,80 : 


- I'^OpFChargedSeivtxo?.- defined in the same way as with ppchargedSeivtxw, which 


is defined in Section 4.1 but using here a different annulus of 0.02 < Ai7 < 0.3. 






— ISOpFChargedWorstvtxos- defined in the same way as with ISOppchargedSeivtxos but 
using the vertex with the maximum isolation for photon momentum direction and 
isolation computation. 

— ISOpFPhoton- the pt sum of particle-flow photon within annulus 0.07 < AR < 0.3 
{AR < 0.3 and \Ap\ > 0.015) from the photon momentum direction for photon 
in the endcap (barrel). The photon momentum direction used in this case is 
obtained with respect to the vertex associated with each particle-flow photon. 

• Auxiliary variables: 

— pEvent- the estimate of transverse energy per unit area in the p-cf) plane contributed 

by the pileup interactions and underlying-event effects in the event. It is the 
median constructed using the kt algorithm 

— psc- pseudorapidity of the ECAL supercluster measured with respect to the origin 
of the detector coordinate. 

— Esc- energy deposit in the ECAL supercluster. 
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The shower shape variables and the isolation variables are used as they are related to the 
two intrinsic differences between a prompt photon and a fake photon, respectively. One is 
that the shower of a fake photon is wider on average since it is the combined shower of 
the two photons from a neutral meson decay. The other is that the isolation for a fake 
photon is larger due to the traces of other fragments of the associated jet leaving in the 
detector around the photon supercluster. The auxiliary variables are included such that the 
distributions of shower shape and isolation variables are used differentially as functions of 
pileup contamination measured by pEvent and photon kinematics measured by psc and Eso¬ 
la order to reduce the photon kinematic dependence of the photon identification BDT and 
the associated mass dependence in the diphoton BDT, explicit use of kinematic differences 
between prompt photons and fake photons in the training sample is avoided, by reweighting 
the 2D PT-psc distribution of the signal to that of the background. 

The distributions of the input variables for the signal and background training samples 
after the reweighting are shown in Figure 4-4, Figure |4^ Figure 4-6| and Figure 4-7 The 
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discontinuities in the V^OppchargedSeivtxoz distribution and in the PFChargedWovstvtxo?. dis¬ 
tribution for the background, shown in Figure [4^ are due to the cut ppchargedSeivtxm < 
4 GeV in the preselection. The signihcant drops in the rjsc distribution around the transition 
regions between ECAL barrel and endcap for both signal and background, shown in Figure 


4-7, are due to the acceptance cut which removes the photons in the region 1.4442 < \risc\ 


< 1.566. 

The reweighting is only done for the training process, but not for the evaluation of BDT 
output and performance as introduced below. 


4.4.3 Output and Performance 

The photon identihcation BDT output is a score named IDBDT assigned to each photon 
which ranges from —1 to 1. The higher the score assigned to a photon, the more likely 
the photon is a prompt photon rather than a fake photon. Figure 4-8| shows the IDBDT 
distributions of the signal (blue) and background (red) training samples (solid circles), and 
of the corresponding testing samples (hollow circles), separately for photons in the barrel 
(left) and for photons in the endcap (right), at 7 TeV (up) and 8 TeV (down). The testing 
signal sample consists of prompt photons from a Monte Carlo simulation of 77 —)■ 77 events 
at a Higgs mass of 121 GeV (124 GeV) at 7 TeV (8 TeV). The testing background sample 
consists of fake photons from Monte Carlo 7 -|- jet events not used for training. Both training 
and testing samples of photons for the plots pass the preselection with pp > 25 GeV. Good 
agreement between the distributions of the testing samples and those of the training samples 
is shown, which verihes the statistical stability of the IDBDT. 

The photon identihcation BDT performance is evaluated using the testing samples. The 
curves of overall background efficiency versus signal efficiency, corresponding to IDBDT 
cuts, for photons in the barrel and endcap at 7 TeV (8 TeV) are shown on the left (right) in 


Figure 4-9 As a reference, the background efficiency for photons in the barrel (endcap) at 


7 TeV (8 TeV), at 80% signal efficiency, is listed in Table 4.4 The corresponding differential 


signal and background efficiencies versus psc, Pt and Nytx are shown in Figure |4-10[ The 
efficiencies versus psc and pp are reasonably hat for photons in the barrel (endcap). This is 
a desirable feature due to the inclusion of psc and Esc info the input variables, and the 2D 
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Figure 4-4: The distributions of photon identihcation BDT input variables (hrst row), 
coYirji^p (second row), and E 2 X 2 /E 5 X 5 (third row) for signal prompt photons (blue) and back¬ 
ground fake photons (red) in the barrel (left) and in the endcap (right) from pp collisions at 
7 TeV (hollow) and 8 TeV (solid). The photons are from the training samples passing the 
preselection with pt > 15 GeV and after pr-Vsc reweighting. 
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Figure 4-5: The distributions of photon identihcation BDT input variables i ?9 (hrst row), 
SC r^-Width (second row), and SC 0-Width (third row) for signal prompt photons (bine) and 
background fake photons (red) in the barrel (left) and in the endcap (right), along with the 
distribntion of Preshower (fourth row) for photons in the endcap only, from pp collisions 
at 7 TeV (hollow) and 8 TeV (solid). The photons are from the training samples passing the 
preselection with pt > 15 GeV and after pr-Vsc reweighting. 

93 
































Figure 4-6: The distributions of photon identihcation BDT input variables ISOppphoton 
(hrst row), ISOppchargedSeivtxos (secoud row), and ISOppchargedWorstvtxos (third row) for 
signal prompt photons (blue) and background fake photons (red) in the barrel (left) and 
in the endcap (right) from pp collisions at 7 TeV (hollow) and 8 TeV (solid). The photons 
are from the training samples passing the preselection with pp > 15 GeV and after pp-Vsc 
reweighting. 
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Figure 4-7: The distributions of photon identification BDT input variables rjsc (first row), 
Esc (second row), and pcvent (third row) for signal prompt photons (blue) and background 
fake photons (red) in the barrel (left for Esc and pEvent) and in the endcap (right for Esc 
and pEvent) from pp collisions at 7 TeV (hollow) and 8 TeV (solid). The photons are from the 
training samples passing the preselection with px > 15 GeV and after px-psc reweighting. 
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Figure 4-8: The distributions of IDBDT for signal prompt photons (blue) and background 
fake photons (red) in the barrel (left) and endcap (right) from 7 TeV (top) and 8 TeV 
(bottom) pp collisions. The photons are from the training samples (solid) and the testing 
samples (hollow) passing the preselection with px > 25 GeV. 
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Pt-Vsc reweighting in the training. The efficiencies are also reasonably flat as a function of 
Nvtx, which is expected as a result of using pEvent as an input variable. 




Signal Efficiency 


Signal Efficiency 


Figure 4-9: The efficiency for background fake photons versus the efficiency for signal prompt 
photons in the barrel (magenta) and in the endcap (black) from 7 TeV (left) and 8 TeV 
(right) pp collisions. The photons are from testing samples passing the preselection with 
Pt > 25 GeV. 


Table 4.4: The efficiency of background fake photons at the signal prompt photon efficiency 
80% . The photons are from testing samples passing the preselection with px > 25 GeV. 



7 TeV (%) 

8 TeV (%) 

Barrel 

9.0 ± 0.2 

10.0 ± 0.2 

Endcap 

12.4 ± 0.2 

13.0 ± 0.2 


The discrepancy between the IDBDT distributions of Monte Garlo simulation and data 
due to the imperfect simulation of detector response is evaluated using electrons from Z —)■ 
events and photons from Z —)■ events, and is treated as a systematic uncertainty 

in the Higgs signal extraction. Photons are required to pass the cut IDBDT > —0.2, since 
in IDBDT region below —0.2 the agreement between data and Monte Garlo simulation is 
relatively poor, and signal to background ratio is very small. The efficiency for the cut on 
the prompt photons is rounded to 1, as well as the efficiency scale factor from Monte Garlo 
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Figure 4-10: The differential efficiencies versus rjsc (top), px (middle) and Nytx (bottom) 
for signal prompt (left) and background fake (right) photons at 7 TeV (hollow) and 8 TeV 
(solid) pp collisions. The efficiencies are evaluated at overall signal efficiency 80%. The 
efficiencies versus psc are in blue for both photons in the barrel and endcap. The efficiencies 
versus px and Nytx are in magenta for the photons in the barrel and in black for photons 
in the endcap. The photons are from testing samples passing the preselection with px > 
25 GeV. 
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simulation to data. A shift of ±0.01 of the Monte Carlo IDBDT is shown to cover the 


discrepancy in the region IDBDT > —0.2. Figure 4-11 ^ shows the IDBDT distributions 
for electrons in the barrel, from Z —>■ e’''e“ events from data (points) and Monte Carlo 
simulation (histogram) at 8 TeV with < 15 (left) and A’uta, > 15 (right), passing the 
preselection with inverted electron veto and IDBDT > —0.2. Good agreement between data 
and Monte Carlo simulation is observed within the ±0.01 variation band of the Monte Carlo 
simulation. 



Figure 4-11: IDBDT (Photon ID BDT score) distributions for electrons in the barrel, from 
Z —)■ e’''e“ events from data (points) and Monte Carlo simulation (histogram) at 8 TeV 
with Nytx ± 15 (left) and Nytx > 15 (right) are shown. Electrons are required to pass the 
preselection with inverted electron veto and IDBDT > —0.2. The ±0.01 shift of the Monte 
Carlo distribution is shown as the red band. 


4.5 Diphoton BDT 

The diphoton BDT is trained to provide each preselected diphoton pair a score, as a measure 
of its expected SjB under the signal diphoton mass peak in the existence of the Higgs boson, 
which is used for event classification as described in Chapter 
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4.5.1 Training Samples 

The training is performed on Monte Carlo simulation of reconstructed diphoton events which 
pass the preselection and IDBDT > —0.2 for both photons. The signal sample consists of 

—)■ 77 events at a Higgs mass of 123 GeV, with all four production processes weighted by 
cross section. The background sample consists of a proper mixture of prompt diphoton, 7 + 
jet and dijet events. The trainings are performed separately for events at 7 TeV and 8 TeV. 

4.5.2 Input Variables 

The input variables are described as following: 

• Diphoton mass resolution variables: 

— the diphoton mass resolution estimator assuming the correct vertex is 
selected. It is the sum in quadrature of the per-photon energy resolution estima¬ 
tors of the leading and sub-leading photon {aE/Ey^ and {aE/E)'^'^ as: 

(<r^hn)a = {[o E + {(a e I (4.7) 

— the diphoton mass resolution estimator assuming the wrong vertex 
is selected. It is the sum in quadrature of {oml'^R and the mass resolution 
contributed by vertex selection /m) as: 

{a„^/m)w = \J -h (4.8) 

where is computed by propagating the uncertainty of the distance be¬ 

tween the selected vertex and true vertex, approximated by ^/2 times the average 
standard deviation of the pp interaction region in ; 2 . 

“ Pvtx'- the probability of the selected vertex being the right vertex estimated from 
vertex probability BDT. 

• Photon identihcation variables: 
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— IDBDT'^^: the score assigned to the leading photon from the photon identihcation 
BDT. 

— IDBDT^^: the score assigned to the snb-leading photon from the photon identih¬ 
cation BDT. 


• Diphoton kinematics variables: 

— the transverse momentnm of the leading photon divided by the diphoton 
mass. 

— /m^^: the transverse momentnm of the snb-leading photon divided by the 
diphoton mass. 

— the psendorapidity of the leading photon momentnm. 

— the psendorapidity of the snb-leading photon momentnm. 

— cos{A(j)^^): the cosine of the separation in the azimnthal angle between the leading 
and snb-leading photon. 


The diphoton mass resolntion variables are not nsed directly in the training bnt combined 


into a weight 1/aEff as in Eqnation 4.9, where UEff is the effective diphoton mass resolntion 
estimator. The if —?• 77 events in the training sample are weighted by I/cte//, snch that 
the events with better resolntion get higher weights and appear more signal like. 




PVtx 


+ 


1 - Pvtx 


C^m/m)R {(Tm/m)w 


(4.9) 


The distribntions of the inpnt variables for data and Monte Carlo signal and backgronnd 
events at 8 TeV, which pass the preselection and IDBDT > —0.2 for both photons, are shown 


in Fignre |4-12| and Fignre |4-13[ Photon energy corrections are applied to both data and 
Monte Carlo events, and additional corrections are applied to Monte Carlo events inclnding 
efficiency scaling and pilenp reweighting. The signal consists of if —>■ 77 events at a Higgs 
mass of 124 GeV, which is later used for the event classihcation optimization. The data 
and Monte Carlo background events in the signal region 120 GeV < < 130 GeV are 

removed for the data and Monte Garlo background comparison. Good agreement between 
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data and Monte Carlo simnlation of backgronnd is shown. The remaining discrepancy makes 
the performance of the diphoton BDT snb-optimal bnt does not affect the correctness of the 
analysis as the backgronnd is evalnated from data for the hnal statistical analysis. 


4.5.3 Output and Performance 


The distribntions of the diphoton BDT ontpnt named DiphotonBDT for data, Monte Carlo 
signal and backgronnd events, which pass the preselection and IDBDT > —0.2 for both 
photons, are shown on the left in Figure |4-14 The DiphotonBDT ranges from —1 to 1. 
The number of if —)■ 77 events over the number of background events in each bin increases 
with the DiphotonBDT as expected. The if —>■ 77 events from the production processes 
VBF, VH and itH tend to have higher score than the events from ggH production process. 
This is due to the fact that VBF, VH and itH events have higher Higgs pt and so higher 
cos(A0.y.y) on average than ggH events. Also the BDT assigns higher score on the events 
with higher as the if —)■ 77 events have higher cos(A 0 .y.y) on average than the 


background events, as shown at the bottom plot in Figure 4-13 The data and Monte Carlo 
background DiphotonBDT distributions are compared with signal region 120 GeV < < 

130 GeV removed. The Monte Garlo background in general describes the data well. The 


contributions from each background component are shown on the right in Figure [4-14[ The 
average score increases with the number of prompt photons in the background as expected. 
The discrepancy between data and Monte Garlo background in the high score region is due 
to the discrepancy between the actual and simulated kinematics for the prompt diphoton 
background, but this does not influence the correctness of the analysis as explained above. 

The uncertainty of IDBDT is propagated to DiphotonBDT by shifting the IDBDT of 
both photons by ±0.01. The differences between the two varied DiphotonBDT distributions 
corresponding to the IDBDT shifts and the original DiphotonBDT distribution are shown 
as the red error bands for Monte Garlo signal and background on the left and right of Figure 


4-15 The uncertainty of (Te/E is propagated to the DiphotonBDT by scaling the (Te/E 
of both photons by ±10%, and the corresponding error bands for Monte Garlo signal and 


background are shown on the left and the right in Figure 4-16 The uncertainty of the 
DiphotonBDT due to diphoton kinematics is taken into account by varying the Higgs pt 
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^=8TeV L= 19.7fb"' ^ = 8TeVL = 19.7 fb’ 






Figure 4-12: The distributions of the diphoton BDT input variables IDBDT'>^^ (top left), 
IDBDT'’'^ (top right), (middle left), {am/m)w (middle right), andp^ta; (bottom) for 

data (points), Monte Carlo background (histogram with blue band for statistical uncertainty) 
consisting of prompt diphoton, 7 -|- jet and dijet events weighted by cross section, and Monte 
Carlo signal (red line) consisting of if —)■ 77 events at a Higgs mass of 124 GeV with all 
four production processes weighted by cross section, at 8 TeV. The data and Monte Carlo 
background events in the signal region 120 GeV < < 130 GeV are removed. All events 

pass the preselection with IDBDT > —0.2 for both photons. Photon energy corrections 
are applied to both data and Monte Garlo events, and additional corrections are applied to 
Monte Garlo events including efficiency scaling and pileup reweighting. 
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Figure 4-13: The distributions of the diphoton BDT input variables (top left), 

(top right), (middle left), (middle right), and cos{A(j)^^) (bottom) for data 
(points), Monte Carlo background (histogram with blue band for statistical uncertainty) 
consisting of prompt diphoton, 7 -|- jet and dijet events weighted by cross section, and 
Monte Carlo signal (red line) consisting of if —)■ 77 events at a Higgs mass of 124 GeV with 
all four production processes weighted by cross section, at 8 TeV. The data and Monte Carlo 
background events in the signal region 120 GeV < < 130 GeV are removed. All events 

pass the preselection with IDBDT > —0.2 for both photons. Photon energy corrections 
are applied to both data and Monte Garlo events, and additional corrections are applied to 
Monte Garlo events including efficiency scaling and pileup reweighting. 
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= 8 TeV L = 19.7 fb“' = 8 TeV L = 19.7 fb ' 



Figure 4-14: Left: the distributions of DiphotonBDT for data (points), Monte Carlo back¬ 
ground (blue line) consisting of prompt diphoton, 7 -|- jet and dijet events weighted by cross 
section, and Monte Carlo signal (stacked histogram) consisting of iL —>■ 77 events at a Higgs 
mass of 124 GeV with all four production processes weighted by cross section, at 8 TeV. 
Right: the distributions of DiphotonBDT for data (points) and Monte Carlo background 
(stacked histogram) consisting of prompt diphoton, 7 -|- jet and dijet events weighted by 
cross section, at 8 TeV. The data and Monte Carlo background events in the signal region 
120 GeV < < 130 GeV are removed. All events pass the preselection with IDBDT > 

—0.2 for both photons. Photon energy corrections are applied to both data and Monte Garlo 
events, and additional corrections are applied to Monte Garlo events including efficiency 
scaling and pileup reweighting. 
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and rapidity by their theoretical uncertainties. 


is = 8 TeV L = 19.7 fb“' is = 8 TeV L = 19.7 fb“' 



Figure 4-15: Left: the distribution of DiphotonBDT for Monte Carlo signal (stacked his¬ 
togram) consisting of iL —)■ 77 events at a Higgs mass of 124 GeV with all four production 
processes weighted by cross section at 8 TeV. The variation of the DiphotonBDT distribu¬ 
tion under the shift of IDBDT of both photons by ±0.01 is shown as the red band on top 
of the stacked histogram. The corresponding ratios between the varied distributions and 
the original distribution are shown as the red band at the bottom. Right: the distributions 
of DiphotonBDT for data (points), Monte Carlo background (stacked histogram) consisting 
of prompt diphoton, 7 ± jet and dijet events weighted by cross section at 8 TeV, and the 
variation of Monte Carlo background DiphotonBDT distribution under the shift of IDBDT 
of both photons by ±0.01 (red band). At bottom, the ratio between the DiphotonBDT 
distributions of data and original Monte Carlo background (points), and the ratios between 
the varied Monte Carlo background distributions and the original Monte Carlo background 
distribution (red band) are shown. The data and Monte Carlo background events in the 
signal region 120 GeV < m-y-y < 130 GeV are removed. All events pass the preselection with 
IDBDT > —0.2 for both photons. Photon energy corrections are applied to both data and 
Monte Carlo events, and additional corrections are applied to Monte Carlo events including 
efficiency scaling and pileup reweighting. 
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Js = 8 TeV L = 19.7 fb' ^ = 8 TeV L = 19.7 fb ' 





Figure 4-16: The same with Figure [4-15| except that the red bands represent the variation 
under the scaling of (Te/E of both photons by ±10% instead. 
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Chapter 5 


Tags of Higgs Production Processes 


We assign tags of Higgs production processes to each diphoton event passing the preselection 
and IDBDT > —0.2 for both photons. The tags are determined by identifying the signatures 
associated with VBF, VH and itH processes, which are the presence of additional objects 
such as jets, leptons or transverse missing energy reconstructed following the descriptions 


in Section 2.2 Further energy correction and identihcation of these objects are provided in 


Section 5.1 The criteria for the VBF, VH and ttH tags based on these objects are given 


in Section |5.2[ Section |5.3| and Section |5.4[ respectively. The variables of these objects and 
photons used for the tagging are dehned in Appendix If none of the tagging criteria is 
satished, the diphoton event is labeled as “untagged”, equivalent to the tag of ggH process. 


5.1 Objects for Higgs Production Tagging 


5.1.1 Jets 

For jets, their energies are corrected following a multi-step procedure. A set of selection cuts 
are applied afterwards to jets in the events at 8 TeV, in order to remove the “fake jets” due 
to clustering of random particles from pileup interactions, which are negligible for the events 


at 7 TeV. The corrections 100,103-105 are listed below, and the selection cuts are listed in 
Table 15. II 


The subtraction of the pileup contamination estimated by pEvent times the jet area. 
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• A relative scale correction to achieve a uniform jet energy response as a function of jet 
pseudorapidity rjK 

• An absolute scale correction to match the original parton energy as a function of jet 
transverse momentum as derived from dijet, 7 + jet and Z + jet events. 

• A residual scale correction between data and Monte Carlo simulation. 


Table 5.1: The selection for jet identihcation. 


w\ 



< 2.5 

< 0.21og(Ayte-0.64) 

< 0.06 

2.5 < 17^1 < 2.75 

< 0.31og(Ayte-0.64) 

< 0.05 

2.75 < IpJj < 3 

- 

< 0.05 

3 < y\ < 4.7 

- 

< 0.055 


5.1.2 Electrons 


For electrons, a set of selection cuts are applied to remove the electrons from jets or electrons 


from converted photons, as listed in the Table 5.2 


Table 5.2: The selection for electron identihcation. 


d%v 

< 0.2 mm 

dl 

< 2 mm 

PconvVtx 

< 10-^ 

^Miss 

< 1 

EleMVA 

> 0.9 

ISO/^e/Pf /CorrPFCombineOS 

< 0.15 


5.1.3 Muons 

For muons, a set of selection cuts are applied to reject the background muons including 
muons from hadron decays, beam-halo muons induced by the accelerator and cosmic muons, 
as listed in Table 15.31 
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Table 5.3: The selection for muon identification. 


^Pixel 

> 0 

^T RK Layper 

> 5 

^MuonCharaber 

> 0 

^Matching 

> 1 

^xy 

< 2 mm 


< 5 mm 

X^NDF 

< 10 

^SO JlQlBetaPuCorrPFCombineOA 

< 0.2 


5.1.4 Transverse Missing Energy 

For transverse missing energy, the difference in the magnitude MET between data and Monte 
Carlo, due to the imperfect simulation of detector effects, is corrected by smearing the 
Monte Carlo jet energy to data in addition to the jet energy correction mentioned above. A 
further correction is applied to both data and Monte Carlo simulation to achieve a uniform 
distribution of the azimuthal angle of transverse missing energy. 


5.2 VBF Tag 

The criteria for the VBF tag are based on the feature of two energetic jets with large 
separation in rj. The VBF candidates are hrst preselected from the diphoton events by 
applying a set of loose cuts on dijet kinematics. Each VBF candidate is assigned a score 
from a so called combined BDT, measuring how likely it is a real VBF event rather than 
a background event or a ggH event with two jets from ISR. The VBF tagged events are 
identihed from the VBF candidates and further classihed according to the combined BDT 
score as described in Chapter To train the combined BDT, a dijet-diphoton kinematic 
BDT is trained beforehand, which provides a kinematic discriminator between the VBF 
events and both the background events and the ggH events. The combined BDT is then 
trained using the output of the dijet-diphoton kinematic BDT, along with the output of the 
diphoton BDT as a measure of the diphoton quality, and pJJ correlated to the outputs 
of both BDTs. The dijet preselection cuts for the VBF candidates, and the details about 
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the dijet-diphoton kinematic BDT and the combined BDT are provided below. 


5.2.1 Dijet Preselection 


The dijet kinematic cuts to select the VBF candidates are summarized in Table 5.4 


Table 5.4: The dijet kinematic cuts for the selection of VBF candidates. 


Pt 

> 30 GeV 

Pt 

> 20 GeV 

\gH\ 

< 4.7 

\gF\ 

< 4.7 

^jj 

> 250 GeV 


5.2.2 Dijet-Diphoton Kinematic BDT 

The dijet-diphoton kinematic BDT is trained separately on Monte Carlo events at 7 TeV 
and 8 TeV, which pass a looser diphoton and dijet preselection to increase the number of 
training events. The requirements on the transverse momenta and IDBDTs of photons as 
well as dijet kinematics are loosen as: 

• /m.y.y > 1/4, > 1/5, IDBDT'’'^ > —0.3, IDBDT'^^ > —0.3. 

• p^^ > 15 GeV, p!^ > 10 GeV, rrijj > 75 GeV. 

The signal sample consists of VBF events at a Higgs mass of 123 GeV. The background 
sample consists of background diphoton, 7 -|- jet and dijet events weighted by cross section, 
and ggH events at a Higgs mass of 123 GeV. The contribution of the ggH events in the 
background sample is inflated by applying a weighting factor about 200. The weighting 
factor is chosen such that the BDT uses more features discriminating between the VBF and 
ggH events, while keeps good distinguishment between the Higgs and background events. 
The input variables are the following: 

• All variables for the dijet kinematic cuts for VBF candidate selection. 
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p7p/m^^\ diphoton transverse momentum divided by the diphoton mass. 


1^77 ~ separation between the diphoton pseudorapidity and the average 


pseu¬ 


dorapidity of the dijet 106 


separation in the azimuthal angle between dijet and diphoton. The value is 
set as the maximum between and vr— 0.2 to avoid large theoretical uncertainty 

on the cross section of ggH events with two jets from initial state radiation in the phase 


space where A(j)jj^^^ is close to tt 45,107 . 


The output is a score assigned to each event ranging from —1 to 1, which increases with the 
compatibility of the event kinematics to the VBF kinematics. 


5.2.3 Combined BDT 

The combined BDT is trained separately for events at 7 TeV and 8 TeV. The signal sample 
is the same for the dijet-diphoton kinematic BDT. The background sample consists of the 
same background events for the dijet-diphoton kinematic BDT but not the ggH events, 
to achieve good discrimination between the Higgs events and the background events. The 
training variables are the outputs of the dijet-diphoton kinematic BDT and the diphoton 
BDT, along with p7p/m^^. 

The output of the combined BDT, named CombinedBDT, ranges from —1 to 1. The 
corresponding distributions for Monte Carlo 77 —?■ 77 events at a Higgs mass of 124 GeV 
selected as VBF candidates are shown on the left of Figure |5-1[ and the corresponding 
distributions for data versus Monte Carlo background events, with events in the signal region 
120 GeV < < 130 GeV removed, are shown on the right. Photon and jet energy 

corrections are applied to both data and Monte Carlo events, and additional corrections 
are applied to Monte Carlo events including efficiency scaling and pileup reweighting. The 
number of VBF events over the number of background events or the number of ggH events 
in each bin increases with the CombinedBDT as expected. The background Monte Carlo 
in general describes the shape of the data. The granularity of the comparison is limited by 
the number of Monte Carlo 7 -|- jet and dijet events, and the spikes are due to the Monte 
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Carlo events with large weight. As explained previously, though the discrepancy between 
the Monte Carlo background events and data makes the CombinedBDT sub-optimal, it does 
not influence the correctness of the analysis. 



-1.0 -0.5 0.0 0.5 1.0 

CombinedBDT 



Figure 5-1: Left: the distribution of CombinedBDT for Monte Carlo signal (stacked his¬ 
togram) consisting of iL —)■ 77 events at a Higgs mass of 124 GeV with all four production 
processes weighted by cross section at 8 TeV. Right: the distributions of CombinedBDT 
for data (points) and Monte Carlo background (stacked histogram) consisting of prompt 
diphoton, 7 -|- jet and dijet events weighted by cross section at 8 TeV. The data and Monte 
Carlo background events in the signal region 120 GeV < < 130 GeV are removed. All 

events pass the preselection and IDBDT > —0.2 for both photons. Photon and jet energy 
corrections are applied to both data and Monte Garlo events, and additional corrections are 
applied to Monte Garlo events including efficiency scaling and pileup reweighting. 


5.3 VH Tag 


The criteria for the VH tag are based on the signatures from the decays of the W oi Z 
boson. There are three sub-tags for different decay modes: 


• Lepton (electron or muon) tag for the leptonic W decay or Z decay. 


• Dijet tag for the hadronic W decay or Z decay. 
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• MET tag for Z decaying into two neutrinos, or leptonic W decay with the lepton lost 
from reconstruction or outside of acceptance. 

The tagging criteria optimized for the sensitivity of the VH signal are introduced below. 


5.3.1 VH Lepton Tag 

The VH lepton tag is further divided into tight and loose tags according to the number of 


leptons Niep and MET as dehned in Table 5.5 


Table 5.5: The dehnition of tight and loose VH lepton tags. 



Aflep 

MET 

Tight 

2 

- 

1 

> 45 GeV 

Loose 

1 

< 45 GeV 


The requirements for these tags are summarized in Table 




which include: 


• A set of kinematic cuts on leptons. The dilepton mass mu is required to be close to 
the Z mass since the leptons are supposed to come from the Z decay. 


• Requirements for photons: 

— A pt cut on the leading photon higher than the preselection. This is due to the 
higher Higgs pt and so higher leading photon pt on average for the VH events 
than for the ggH events. 

— Large enough distance in AR from the photon to electron, electron track and 
muon. This is to reject the photon from lepton radiation and electron faking 
photon, which reduces the dominant background from PE + 7 and Z + 7 . 

— The photon-electron mass m^^^, away from the Z mass. This is to reject the 
electron faking photon from Z —)■ 

• Less than three jets with > 20 GeV, \r]^\ < 2.4 and Ai? > 0.5 from any lepton or 
photon. This is to reject contamination from the itH process. 
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Table 5.6: The requirements for VH lepton tag. 


Nlep 

1 

2 

Pt 

> 20 GeV 

> 10 GeV 

Pt 

> 20 GeV 

> 10 GeV 

\v^\ 

< 2.4 


|77"| < 1.4442 or 1.566< < 2.5 

mu 

- 

70 GeV < mu < 110 GeV 


> 1 

> 0.5 

^R-y,e 

> 1 

^R'y^etrk 

> 1 

\m^^e - Mz\ 

> 10 GeV 

p'^ 

> 3/8 

N, 

< 3 

- 


5.3.2 VH Dijet Tag 


The VH dijet tagged events are selected from the diphoton events with two jets, following the 


requirements summarized in Table 5T For the events with more than two jets, the leading 
and sub-leading jets in px are considered. Among the requirements, the dijet mass rrijj is 
required to be close to the Z{W) mass since the jets are supposed to come from the Z{W) 
decay. The cosine of the angle 9* between the diphoton momentum in the center-of-mass 
frame of diphoton-dijet and the total momentum of diphoton-dijet in the lab frame is used. 
Its distribution is flat for VH events while peaking at |cos(6'*)| = 1 for background events. 


Table 5.7: The requirements for VH dijet tag. 


Ji 

Pt 

> 40 GeV 

47 

> 40 GeV 

W^\ 

< 2.4 

w \ 

< 2.4 

mjj 

60 GeV < mjj <120 GeV 


< 0.5 


> 1/2 

p7p / 

> 13/12 
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5.3.3 VH MET Tag 


The VH MET tagged events are selected from the diphoton events with large MET. The 


selection requirements are summarized in Table 5.8 Among the requirements, large separa¬ 
tion in the azimuthal angle between the diphoton and MET (j)-y^^MET is required, because of 
the momentum balancing between the diphoton and MET in the VH events with Z decaying 
into two neutrinos, or W leptonic decay with the lepton lost from reconstruction or outside 
of acceptance. An upper bound is put on the separation in the azimuthal angle between the 
diphoton and the leading jet, in order to reduce the contamination from the MET caused by 
the inaccurate measurement of jet energy when the jet and the diphoton is back to back. 


Table 5.8: The requirements for VH MET tag. 


MET 

> 70 GeV 


> 2.1 


< 2.7 


> 3/8 


5.4 ttH Tag 

The criteria for the iiH tag are based on the signatures from the decays of the ti. There are 
two sub-tags for different decay modes: 

• Lepton (electron or muon) tag for two or one leptonic W decay. 

• Multijet tag for two hadronic W decays. 

The tagging criteria optimized for the sensitivity of the itH signal are introduced below. 


5.4.1 ttH Lepton Tag 


The UH lepton tagged events are selected from the diphoton events with at least one lepton. 


following the requirements summarized in Table 5.9 Among the requirements, the pt cut 
on the leading photon is increased with respect to the preselection because of the higher 


117 












Higgs pt and so the higher leading photon px on average for the iiH events than for the 
ggH events. For the requirement of the number of jets (b-jets), the jets (b-jets) with plp > 
25 GeV, \p^\ < 2.4, AR > 0.5 from any lepton are counted. 


Table 5.9: The requirements for ttH lepton tag. 


Pt 

> 20 GeV 

Pt 

> 20 GeV 

\p^\ 

< 2.4 

w\ 

\p^\ < 1.4442 or 1.566< \p^\ < 2.5 


> 0.5 


> 1 

^f^'y,etrk 

> 1 

p'^ / 

> 1/2 

N, 

> 1 

Nb-, 

> 0 


5.4.2 ttH Multijet Tag 


The ttH Multijet tagged events are selected from the diphoton events with multiple jets, 


following the requirements summarized in Table |5.10[ For the requirement of the number of 
jets (b-jets), the jets (b-jets) with plp > 25 GeV, l < 2.4 are counted. 


Table 5.10: The requirements for ttH multijet tag. 


TV, 

> 4 

Nb-j 

> 0 

plj^ / 

> 1/2 
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Chapter 6 


Event Classification 


We classify the diphoton events passing the preselection and IDBDT > —0.2 for both photons 
into the tagged and the untagged Higgs production process classes. The events in VBF 
tagged classes are selected from VBF candidates passing the dijet kinematic selection, and 
are classified in terms of the CombinedBDT. The corresponding class boundaries are chosen 
to minimize the expected uncertainty of the signal strength for the VBF + VH processes 
fiVBPyH, sensitive to the Higgs coupling to bosons. The events in the VH and itH tagged 
classes are the VH and iiH tagged events, which pass the additional DiphotonBDT cuts to 
improve the VH and iiH sensitivity. The untagged events are classified into the untagged 
classes in terms of the DiphotonBDT, and the corresponding class boundaries are chosen 
to minimize the expected uncertainty of the overall signal strength hh- The optimization 
of the class boundaries for the VBF tagged classes and the untagged classes is described in 


Section 6.1 The final tagged and untagged event classes are summarized in Section 6.2 


6.1 Boundary Optimization for VBF Tagged Classes 
and Untagged Classes 

The boundaries on CombinedBDT and DiphotonBDT for the VBF tagged classes and the 
untagged classes are optimized separately for events at 7 TeV and 8 TeV, using Monte Carlo 
diphoton events passing the preselection and IDBDT > —0.2 for both photons, independent 
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from the training sample. The signal sample consists of iif —)■ 77 events for a Higgs mass of 
124 GeV (121 GeV) with all fonr production processes weighted by cross section at 8 TeV 
(7 TeV). The background sample consists of prompt diphoton, 7 + jet and dijet events 
not used for BDT training and weighted by cross section. The total number of Monte 
Garlo events are weighted to match the luminosity in data, and corrections on the Monte 
Garlo simulation including photon and jet energy corrections, efficiency scaling and pileup 
reweighting are applied. The background GombinedBDT and DiphotonBDT distributions 


are smoothed using adaptive Gaussian kernel estimations 108 


6.1.1 VBF Tagged Class Optimization 

The GombinedBDT boundaries are hrst optimized on events passing VBF dijet kinematic 
selection. The number of boundaries and the corresponding values for the boundaries are 
adjusted interactively until the decrease of the expected uncertainty is less than 1%. The 
evaluation of the expected uncertainty is based on the prohle likelihood £t on the diphoton 
mass spectra cross all the classes, which follows the procedure as described in Ghapter 
with a simplihed signal model and background model. For each class, the histogram of the 
diphoton mass for Monte Garlo Higgs events is used as the signal model while a power law 
funtion htted from the Monte Garlo background events is used as the background model. The 
variation of the signal model due to systematic uncertainties is not considered for simplicity. 

For the VBF tagged classes, 3 classes for events at 8 TeV and 2 classes for events at 7 TeV 
are determined. Due to the limited number of Monte Garlo prompt diphoton background 
events passing the dijet kinematic selection, the VBF tagged classes for events at 7 TeV are 
determined by matching the acceptance times efficiency for VBF events to those of the VBF 
tagged classes of events at 8 TeV, instead of using the optimization procedure as described 


above. Figure [ 6 M] shows the GombinedBDT distributions, along with the class boundaries 
(dashed lines), in the range of GombinedBDT > 0 , of Monte Garlo signal events (left), Monte 
Garlo background events and data (right) at 8 TeV, passing the preselection and IDBDT 
> —0.2 for both photons, and dijet kinematic cuts. The events in the shaded region below 
the dashed line with the lowest GombinedBDT value is taken away from the VBF tagged 
classes, and used for the selection for the rest of classes. 
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6.1.2 Untagged Class Optimization 

The events with CombinedBDT below the lowest boundary and the events not passing the 
dijet kinematic cuts are used for the optimization of the DiphotonBDT boundaries of the 
untagged classes. The procedure is the same as for the determination of the boundaries of 
the CombinedBDT. The events with DiphotonBDT below the lowest boundary are removed, 
which includes few signal but large number of background events. The removal of these 
events causes a negligible loss in the sensitivity for the Higgs signal and largely simplifies the 
final statistical analysis. 

For the untagged classes, 4 classes for events at 7 TeV and 5 classes for events at 8 TeV 


are determined. Figure |6-2| shows the DiphotonBDT distributions, along with the class 
boundaries (dashed lines), of Monte Carlo signal events (left), Monte Carlo background 
events and data (right) at 8 TeV, passing the preselection and IDBDT > —0.2 for both 
photons. The region below the dashed line with the lowest DiphotonBDT value is removed. 


6.2 Final Event Classes 


The diphoton events, passing the preselection and IDBDT > —0.2 for both photons, are 
first selected into the tagged classes and the rest are selected into the untagged classes. 
The classes are mutually exclusive. In the case that an event satisfies the criteria of more 
than one tagged classes, the class with higher fraction of events from the corresponding 
tagged production process among the selected signal events is chosen in general. The hnal 
event classes including 11 classes for events at 7 TeV and 14 classes for events at 8 TeV are 


summarized in Table 6.1 There are very few ttH lepton tagged events and multijet tagged 


events at 7 TeV, and so these events are combined into a single class. 
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^ = 8 TeV L = 19.7 fb‘ 




Figure 6-1: Left: the distribution of CombinedBDT in the range of CombinedBDT > 0 for 
Monte Carlo signal (stacked histogram) consisting of Ff —)■ 77 events at a Higgs mass of 
124 GeV with all four production processes weighted by cross section at 8 TeV. Right: the 
distributions of CombinedBDT in the range of CombinedBDT > 0 for data (points) and 
Monte Carlo background (stacked histogram) consisting of prompt diphoton, 7 -|- jet and 
dijet events weighted by cross section at 8 TeV. The data and Monte Carlo background events 
in the signal region 120 GeV < m-f^ < 130 GeV are removed. All events pass the preselection 
with IDBDT > —0.2 for both photons and dijet kinematic cuts. Photon and jet energy 
corrections are applied to both data and Monte Garlo events, and additional corrections are 
applied to Monte Garlo events including efficiency scaling and pileup reweighting. The class 
boundaries are labeled as dashed lines. The events in the shaded region below the dashed 
line with the lowest GombinedBDT value are taken away from the VBF tagged classes and 
used for the selection for the rest of classes. 
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\fs = 8 TeV L = 19.7 fb' \fs = 8 TeV L = 19.7 fb^ 



Figure 6-2: Left: the distribution of DiphotonBDT for Monte Carlo signal (stacked his¬ 
togram) consisting of iL —)■ 77 events at a Higgs mass of 124 GeV with all four production 
processes weighted by cross section at 8 TeV. Right: the distributions of DiphotonBDT for 
data (points), Monte Carlo background (stacked histogram) consisting of prompt diphoton, 
7 -|- jet and dijet events weighted by cross section at 8 TeV. The data and Monte Carlo 
background events in the signal region 120 GeV < < 130 GeV are removed. All events 

pass the preselection with IDBDT > —0.2 for both photons. Photon energy corrections 
are applied to both data and Monte Garlo events, and additional corrections are applied to 
Monte Garlo events including efficiency scaling and pileup reweighting. The class boundaries 
are labeled as dashed lines. The events in the shaded region below the dashed line with the 
lowest DiphotonBDT value are dropped. 
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Table 6.1: The event classes listed in the event selection order. The events for each class are 
selected from the preselected diphoton events with IDBDT > —0.2 for both photons. 


Event classes 

Tag 

DiphotonBDT 

CombinedBDT 


ttH Lepton + 
Multijet 

itH Lepton 

or iiH Multijet 

> 0.6 

- 

1 

VH Lepton 

Tight 

VH Lepton Tight 

> 0.1 

- 

£ 
T —1 

id 

VH Lepton 

Loose 

VH Lepton Loose 

> 0.1 

- 

> 

VBF Dijet 0 

VBF Candidate 

- 

> 0.995 


VBF Dijet 1 

VBF Candidate 

- 

> 0.917 && < 0.995 

t- 

VH MET 

VH MET 

> 0.8 

- 


VH Dijet 

VH Dijet 

> 0.6 

- 


Untagged 0 

- 

> 0.93 

- 


Untagged 1 

- 

> 0.85 && < 0.93 

- 


Untagged 2 

- 

> 0.7 && < 0.85 

- 


Untagged 3 

- 

> 0.19 && < 0.7 

- 


iiH Lepton 

iiH Lepton 

> - 0.6 

- 


VH Lepton 

Tight 

VH Lepton Tight 

> - 0.6 

- 

7 

VH Lepton 

Loose 

VH Lepton Loose 

> - 0.6 

- 


VBF Dijet 0 

VBF Candidate 

- 

> 0.94 

S 

VBF Dijet 1 

VBF Candidate 

- 

> 0.82 && < 0.94 

> 

VBF Dijet 2 

VBF Candidate 

- 

> 0.14 && < 0.82 


VH MET 

VH MET 

> 0 

- 

00 

iiH Multijet 

UH Multijet 

> - 0.2 

- 


VH Dijet 

VH Dijet 

> 0.2 

- 


Untagged 0 

- 

> 0.76 

- 


Untagged 1 

- 

> 0.36 && < 0.76 

- 


Untagged 2 

- 

> 0 && < 0.36 

- 


Untagged 3 

- 

> -0.42 && < 0 

- 


Untagged 4 

- 

> -0.78 && < -0.42 

- 
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Chapter 7 


Statistical Procedure for the 
Extraction of the Higgs Signal 


We extract the signal of Higgs boson from the observed diphoton mass spectra of all the 
event classes. For each event class, the models of the expected diphoton mass spectrum of 
the Higgs signal events and that of background events are constructed. The signal model 
is built for each Higgs mass hypothesis niH in the search range [115,135] GeV, using para¬ 


metric functions htted from Monte Carlo simulated events, as described in Section 7.1 The 
background model is built using a set of parametric functions htted directly from data, and 
the uncertainty due to the limited knowledge of the true background function is taken into 
account by prohling the choice of functions in the signal extraction, as described in Section 


7.2 The main systematic uncertainties related to the signal model are discussed in Section 


7.3 The statistical procedure for the hnal Higgs signal extraction, based on the simultaneous 


likelihood ht to the diphoton mass spectra over all event classes, is described in Section 7.4 


7.1 Signal Model 

For each event class, there are three steps in the signal model construction. First, the signal 
models of the Standard Model (SM) Higgs boson are built for hve reference Higgs mass 
hypotheses separated by a 5 GeV step, mji^G{115,120,125,130,135} GeV, on Monte Garlo 
simulations with the resolution correction, preselection efficiency scale factors and trigger 
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efficiency applied to match data. Second, the signal model as a function of Higgs mass 
is built through interpolation between the neighboring reference masses of Monte Carlo 
models. Finally, the variations of the signal model are constructed based on the SM model 
with signal strengths or coupling strengths included as free parameters (their values equal 
to one for the SM Higgs boson), which are used in the hnal signal extraction. These three 
steps are described below. 


7.1.1 Signal Model for a Reference Higgs Mass 

For a reference Higgs mass models for the four Higgs production processes are hrst built 
respectively, and then combined according to their cross sections. The model for a particular 
Higgs production process XH, any of ggH, VBF, VH and ttH, is described as follow. The 
combined model is described afterwards. 


Model for a Higgs Production Process 

The model for a Higgs production process XH consists of the expected yield and the diphoton 
mass distribution. 

The expected yield, is the product of luminosity, L, SM Higgs production 

cross section for the process, SM branching ratio of the Higgs decaying to two 

photons, acceptance, X{rri^), and efficiency, 

L (7.1) 


The luminosity is taken from the experimental measurement described in References 109 


110 . The cross section and branching ratio are taken from the LHC Higgs boson Cross 


Section Working Group 45 . The acceptance times efficiency is evaluated on the Higgs Monte 
Carlo sample for the process XH with all the corrections applied, which is the fraction of 
events passing all the selection. 

The expected diphoton mass distribution is modeled by an empirical parametric function, 
which consists of two components, one for the events with right vertex selected while the 
other for the events with wrong vertex selected. The right (wrong) vertex component. 
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with the set of parameters, is parametrized as a 

Gaussian distribution, or a sum of two Gaussians, with the set of parameters 
determined by maximum likelihood £t to the Monte Garlo events: 




=m^ +Amf^^(m^),cr2^^(m^)), (7.2) 


where = rrijj + represents the ith Gaussian 

with mean, standard deviation, represents the 

shift of the mean with respect to the nominal Higgs mass, m'fj, and represents the 

fraction coefficient for the ith Gaussian. The diphoton mass distribution for the production 
process, PxH{fn^^\'^xHim^^)), with the set of parameters, xH{rn^)^ is then the sum of 
right and wrong components, weighted according to the vertex efficiency eR{m'jj) calculated 
from the Monte Garlo events: 


PxH{m^^\'^XH{mH)) = ■ PxHim^'r\^XHimH)) + 

(1 — eji{mjj)) ■ (7.3) 

The complete model for the process, SxHi'n^'y'yl'n^'n)^ is then the product of the expected 
yield and the diphoton mass distribution: 

SxHim^'rl^n) = ^xh{^h) ' PxH{m^'y\'^ XHim^)) (7.4) 


Combined Model 


The model for the SM Higgs boson, S{m^.^\ni^), is constructed as the sum of the models of 
all processes: 

S{m^^\mH) = ^ SxH{'m^'y\mH) ( 7 - 5 ) 


XH&{ggH, VBF, VH,ttH} 


Figure 7-1 shows the Monte Garlo diphoton mass spectrum, along with the fitted distri¬ 
bution (red line), for 77 —)■ 77 at = 125 GeV in 8 TeV untagged 0 class, which has the 
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best resolution among all the event classes. The measures of the resolution, the half of the 
narrowest mass interval containing 68.3% of the area under the distribution (yellow), Ue//, 
and the full width half maximum, FWHM, are 1.04 GeV and 2.16 GeV for this class. The 
corresponding hgures for all the other event classes are provided in Appendix [A] 


Table 7.1 shows the total expected yield, the fraction of contribution from each produc¬ 
tion process (the contribution less than 0.1 % is ignored) and Ue// for if —77 at mu = 
125 GeV for each event class. The expected number of the selected signal events is 475.9 
(96.1) at 8 TeV (7 TeV), corresponding to the acceptance times efficiency 48% (48%). The 
contribution from the corresponding tagged process is dominant in the tagged classes, while 
the contribution from the ggH process is dominant in the untagged class as expected. The 
resolution Ue// increases with the class number of the untagged classes also as expected. 


8 TeV Untagged 0 



Figure 7-1: The 8 TeV untagged 0 class’s diphoton mass spectrum (points) and the fitted 
distribution (red line) of Monte Garlo 77 —)■ 77 events at a Higgs mass of 125 GeV. 

7.1.2 Signal Model as a Function of Higgs Mass 

After building the models at the five Higgs mass hypotheses from Monte Garlo simulations, 
the final signal model for the SM Higgs boson as a function of Higgs mass hypothesis, 
S{m^^\mH)-i is constructed by interpolating between the five masses of the Monte Garlo 
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models. The distribution uses the same functional form as the Monte Carlo models. The pa¬ 
rameters associated with each process, ^= {Nxh xHi'mH)}, are piecewise 

linear functions in rriH'- 


^ ^ ^ , ^XHirni + 5GeV)-V'xHirrii) 

1 XH [^H j — 1 XH\jrii) H-p- ■ — ^i), 


i G {1, 2, 3,4}, ruj G {115,120,125,130}GeV, m* < rriH <mi + 5 GeV. 


(7.6) 


7.1.3 Variations of Signal Model 


The variations of the signal model, with signal strengths or coupling strengths included as 
free parameters, are constructed by modifying the Higgs cross section times its branching 
ratio to two photons (Jxh{^h) ■ These parameters, along with the Higgs mass 

mn, are later measured to quantify the compatibility of the SM Higgs model with respect 


to data. The varied models are summarized as following 45,111 


• S{myy\fiH,^H)- signal model with total signal strength, fin, and Higgs mass, mn, as 
free parameters of interest, in which 
(yxHiran) ■ BH^yyimn) = ■ (^xh (^h) ■ 


• S{myy\figgjjyiHyf^VBF,VH,^H)- signal model with signal strength for ggH and tiH pro¬ 
cesses, fJ^ggH^ttH (sensitive to Higgs coupling strength to fermions), signal strength for 
VBF and VH processes, ^ivbf,vh (sensitive to Higgs coupling strength to bosons), and 
Higgs mass, ttih, as free parameters of interest, in which 

' BH^yyimn) = FggHMH ' (^g^ll{ttH)('^H) ' h), 

VBF(VH){'mH) ■ BH^yyimn) = Fvbf,vh ■ <^v^f(vh)(''^h) ■ Bfl^^^{mH)- 


S{myy\Kv 1 signal model with Higgs coupling strength to bosons, Ky, and 

Higgs coupling strength to fermions, Kf —benchmark parameterization dehned in Ref¬ 
erence 


45 —and Higgs mass, mu-, as free parameters of interest. 


S{myy\Ky, Kg, Tufj)■ sigHul model with effective Higgs coupling strength to photon, By, 
effective Higgs coupling strength to gluon. Kg —benchmark parameterization defined in 


Reference 45 —and Higgs mass, m^, as free parameters of interest. 
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Table 7.1: The expected yield S, the fraction of each production process fggH, fvBF, fvH, 
and the resolution cie// for 77 — )■ 77 at mu = 125 GeV, along with the number of 
background events per GeV at 125 GeV dB/dm^^, S/B and S/\fB for each event class. The 
number of background events under the signal peak B is estimated as dB / dm^^ multiplied 
by 4 aeff. _ 




Expected Higgs Boson at uih = 

125 GeV 

dB / dm^^ 

S/B 

s/Vb 

Event classes 

S 

^ggH 

ivBF 

ivH 


^eff 







(%) 

(%) 

(%) 

(%) 

(GeV) 

(GeV-i ) 




UR Lepton 
+ Multijet 

0.2 

2.9 

1.1 

3.5 

92.5 

1.38 

0.2 

0.18 

0.19 


VH Lepton 
Tight 

0.3 

- 

- 

97.7 

2.3 

1.59 

0.1 

0.47 

0.38 

7 

£ 

VH Lepton 
Loose 

0.2 

3.0 

1.1 

94.9 

1 

1.62 

0.2 

0.15 

0.18 

LO 

VBF Dijet 0 

1.6 

18.1 

81.4 

0.5 

- 

1.43 

0.4 

0.70 

1.06 

> 

VBF Dijet 1 

3.0 

38.1 

59.5 

1.9 

0.5 

1.64 

3.3 

0.14 

0.64 

r- 

VH MET 

0.3 

5.7 

1 

85 

8.3 

1.52 

0.2 

0.25 

0.27 


VH Dijet 

0.4 

28.7 

2.8 

66.4 

2.1 

1.55 

0.5 

0.13 

0.23 


Untagged 0 

5.9 

79.7 

10.0 

9.6 

0.7 

1.12 

11.0 

0.12 

0.84 


Untagged 1 

23.0 

91.9 

4.2 

3.7 

0.2 

1.26 

69.2 

0.07 

1.23 


Untagged 2 

27.2 

91.9 

4.2 

3.8 

0.1 

1.76 

133.5 

0.03 

0.89 


Untagged 3 

34.0 

92.0 

4.1 

3.7 

0.2 

2.32 

311.8 

0.01 

0.63 


UR Lepton 

0.5 

- 

- 

2.8 

97.2 

1.32 

0.1 

0.95 

0.69 


VH Lepton 
Tight 

1.4 

- 

0.2 

96.0 

3.8 

1.60 

0.4 

0.55 

0.88 

1 

VH Lepton 
Loose 

0.9 

- 

1.3 

97.1 

1.6 

1.56 

1.1 

0.13 

0.34 


VBF Dijet 0 

4.4 

17.3 

82.3 

0.3 

0.1 

1.27 

0.7 

1.24 

2.33 

r- 

ai 

VBF Dijet 1 

5.4 

26.0 

73.0 

0.8 

0.2 

1.44 

2.7 

0.35 

1.37 

T—1 

> 

VBF Dijet 2 

13.7 

44.0 

53.1 

2.2 

0.7 

1.56 

21.9 

0.10 

1.17 

VH MET 

1.7 

12.0 

2.3 

74.0 

11.7 

1.58 

1.2 

0.22 

0.62 

00 

UR Multijet 

0.6 

7.5 

1.0 

1.7 

89.8 

1.41 

0.5 

0.21 

0.36 


VH Dijet 

1.6 

31.1 

3.0 

63.3 

2.6 

1.33 

1.0 

0.3 

0.7 


Untagged 0 

5.8 

74.7 

12.3 

11.0 

2.0 

1.04 

4.3 

0.32 

1.37 


Untagged 1 

50.5 

85.1 

7.8 

6.5 

0.6 

1.18 

117.9 

0.09 

2.14 


Untagged 2 

116.5 

91.1 

4.8 

3.9 

0.2 

1.43 

410.9 

0.05 

2.40 


Untagged 3 

151.7 

91.5 

4.4 

3.8 

0.3 

1.99 

856.9 

0.02 

1.84 


Untagged 4 

121.2 

93.2 

3.6 

3.1 

0.1 

2.56 

1395.0 

0.01 

1.01 
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7.2 Treatment of Background for the Signal Extraction 


For each event class, the model of the background diphoton spectrum, the product of the 
expected yield and the diphoton mass distribution, is constructed using parametric functions 
htted to data. The functional forms are chosen to describe the continuously falling character 
of the expected background spectrum. The £t range is 100 GeV < < 180 GeV, such 

that the background under an emerging narrow peak, for any given Higgs mass hypotheses 
within 115 GeV < < 135 GeV, gets constraint from sufficient events in the sidebands 

of the signal region. The differences between the chosen background functions and the 
unknown true function lead to an uncertainty in the extraction of Higgs signal. In our 


previous i/ —)■ 77 analysis 36 , a single background function is chosen for each event class, 
following the criterion that the potential bias of the Higgs results is negligible with respect 
to the statistical uncertainty, at the price of increasing the number of parameters of the 


function and inflating the statistical uncertainty. An updated method 112,113 is used in 
this analysis, which incorporates the uncertainty due to the choice of the functional form 
into the total uncertainty of the Higgs results, and thus avoids the inflation of the statistical 
uncertainty. 


The basic idea of the updated method is: hrst, choose a set of background functions 
describing the data well and that are generic enough to cover the true function. Second, build 
a negative log-likelihood function of Higgs parameter of interest, e.g. signal strength 
for each background function, with a correction term penalizing the increase of the number 
of parameters. Third, construct the envelope negative log-likelihood function by taking the 
minimum value of the individual functions at each from which the best £t fin and the 
associated conhdence interval are obtained. The uncertainty of the background function 
choice is taken into account in the conhdence interval as a result of prohling background 
functions. The implementation of the method is described below, and the performance of 
the method, in terms of the bias of results and the coverage of conhdence interval, is discussed 
afterwards. 
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7.2.1 Selection of the Set of Background Functions 

For each event class, a set of background functions, ■■■, Bri{jn^^\6B^)} with 

Ob^ representing the set of parameters for the ith background function, are chosen from the 
following four function families: 

• Nth order Bernstein polynomial 

NBer(m^^) = ^ rfi = (7.7) 

i=0 ^ ' 

• Nth order exponential sum 

N 

NExp(m^.^) = (7.8) 

i=l 

• Nth order power sum 

N 

NPow(m.,^) = ('^•9) 

i=l 

• Nth order Laurent series 

NLau(m.^.^) = ^ (7-10) 

i=l 


For each function family, starting from the function order iV = 1, except for the Laurent 
series which starts from N = 2 because iV = 1 corresponds to a trivial power law function 
background only hts are performed on data with increasing order N. The goodness of 
the hts is measured using and the so called p-value, the probability of getting a result as 
compatible or less to data than the observed one given that the function under consideration 


is true 114 . If the loose criterion of the ht quality p-value > 0.01 is satished, the function 


is included into the function set. This process keeps going for the {N + l)th order function 
until the higher order function is no longer signihcantly favored by data, quantihed by: 


C 


Pix^ > (-2\n^U) > 0 . 1 . 


C 


N+l 


( 7 . 11 ) 
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In the equation above, is the maximum likelihood for the Nth order function; P{x^ > 
-)obs) is the p-value of the observed —21n ^^^^ for a distribution, with the degree 


of freedom as the difference in the number of parameters between the {N + l)th and Nth order 


function, which is the distribution of —21n ^ ^ in the case that Nfh order function is the 
true function and sufficient number of events is available for htting. The highest {N + l)th 
order function satisfying P{x^ > (—21n ^^^^ )obs) < 0.05 is automatically included into the 
function set without the requirement of goodness of the £t. This function corresponds to the 
true function for the pseudo-experiments used to study the potential biases associated with 


different background functions in the previous if —)■ 77 analysis 36 , and is used as the true 
background function to study the bias and coverage of conhdence interval of the function set 
for the update method. The orders of the hnal input background functions for each event 
class are listed in Table 17.21 


7.2.2 Construction of Envelope Negative Log-Likelihood Function 

After selecting the set of background functions, to extract a Higgs parameter under interest, 
e.g. total signal strength hh at a given Higgs mass hypothesis m^, the so called envelop 
negative log-likelihood function, the envelope function, of the parameter is constructed, with 
signal plus background model on the observed diphoton mass spectrum. The data is binned 
in 320 bins of the diphoton mass with 250 MeV per bin—this choice permits a relatively 
quick extraction process while preserving the precision. 

To construct the envelope function, the likelihood function for individual background 
function, e.g. the likelihood function for the ith background function Ci{^Hi ^Bi), is hrst 
built as a product of Poisson distributions: 

320 

= JJ Poisson I + (7.12) 

i=i 

where rij is the observed number of events in the jth bin of the data, Sj{^HiBrtjj) is the 
expected number of signal events in the jth bin under the m'jh mass hypothesis, which is 
obtained by integrating S'(m.y.y|/r//,m):^) over the bin, and bj^i{6Bi) is the expected number 
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of background events in the jth bin obtained by integrating 

The envelope function —21n££;(ynji^, m^) is then constructed as: 


-21n££;(/i^^,m^) = min {-2\nCi{^H, 


m, 


h 

H-! ^ 


+ ^bJ, 


( 7 . 13 ) 


where 6^ represents the set of values of the background parameters maximizing the 

ith likelihood function at and a given and represents the number of parameters 
of the ith background function, acting as a correction term penalizing the increase of number 
of parameters. For two background functions from the same function family Bn and Bm, 
with Bn having larger number of parameters than Bm, the penalty works in the way that 
the two times negative log-likelihood value after correction for Bn roughly equals to that for 
Bm, if the p-values associated with Bn and Bm are the same. This correction reduces the 
statistical uncertainty, while keeps a small bias and a good coverage of confidence interval 
of fitted signal strength. 

The best £t (jlh is then the fin minimizing —21n££;(/ij:^, m^). The confidence intervals 
are determined from the likelihood ratio —2Aln££;(/i//, m^): 


-2A\nCE{^^H,rnH) 


^'h) 


(7.14) 


For example, the boundary points for the 68.3% confidence interval corre¬ 

spond to: 

-2Aln£i);(/i^-^^“",m^) = -2Aln££;(/i^-^^®+, m^) = 1, (7.15) 


for which the uncertainty of the background function choice is taken into account as a result 
of profiling the background functions. 


7.2.3 Performance 

For each event class, the bias of the best fit fiH, defined as the median difference between 
the measured and true fin relative to the uncertainty, and the coverage of the confidence 
interval are evaluated on toy datasets, which are generated from signal plus background 
model for each background truth function as mentioned above. For untagged classes and 
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production process tagged classes with sufficient large samples, the bias and the deviation of 
the conhdence interval coverage from the nominal value are within 14% and 1% respectively, 
or slightly above, in most cases for the signal region 115 GeV < mn < 135 GeV, which are 
considered as neglegible. For tagged classes with few events, the bias and the deviation of 
the conhdence interval coverage are in general larger, with the maximum value about 30% 
and 10% respectively; and the expected background functions are not so well constrained by 
the data in the sidebands. The inhuence from these classes is negligible since the hnal Higgs 
results are extracted by simultaneous htting over all classes. 
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Table 7.2: The orders of the input background functions for all event classes. 


Event classes 

NBer 

NExp 

NPow 

NLau 

7 

T—1 

> 

r- 

Untagged 0 

12 3 

1 

1 

2 

Untagged 1 

3 

1 2 

1 

2 

Untagged 2 

2 3 

1 2 

1 

2 

Untagged 3 

3 4 5 

12 3 

1 

2 

VBF Dijet 0 

1 

1 

1 

2 

VBF Dijet 1 

12 3 

1 

1 

2 

VH Lepton Tight 

1 

1 

1 

2 

VH Lepton Loose 

1 

1 

1 

2 

VH MET 

1 

1 

1 

2 

VH Dijet 

1 2 

1 

1 

2 

tfH Lepton + Multijet 

1 

1 

1 

2 

8 TeV 19.7 fb-^ 

Untagged 0 

12 3 

1 

1 

2 

Untagged 1 

2 3 

1 

1 

2 

Untagged 2 

3 4 

1 2 

1 

2 

Untagged 3 

4 5 

2 

1 

2 

Untagged 4 

4 5 

2 

1 2 

2 

VBF Dijet 0 

1 

1 

1 

2 

VBF Dijet 1 

1 2 

1 

1 

2 

VBF Dijet 2 

2 3 

1 

1 

2 

VH Lepton Tight 

1 

1 

1 

2 

VH Lepton Loose 

1 2 

1 

1 

2 

VH MET 

1 

1 

1 

2 

VH Dijet 

1 2 

1 

1 

2 

ttH Lepton 

1 2 

1 

1 

2 

ttH Multijet 

1 

1 

1 

2 
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7.3 Systematic Uncertainties Associated with the Sig¬ 


nal Model 


The systematic uncertainties associated with the signal model are considered for the hnal 
Higgs signal extraction. There are two types of uncertainties. One leads to the variations of 
the expected signal yield, and dominates the systematic uncertainty of the signal strength. 
The other leads to the variations of the signal shape, and dominates the systematic un¬ 
certainty of the Higgs mass. These uncertainties are summarized in the following, more 


descriptions are available in Reference 40 . The statistical procedure to incorporate the 


corresponding signal variations into the Higgs signal extraction is described afterwards. 


7.3.1 Systematic Uncertainties Related to the Signal Yield 

There are two kinds of systematic uncertainties influencing the signal yield. The hrst kind 
causes 100% correlated variations of yields of all the event classes under influence. The 
second kind causes migrations of events among classes and so —100% correlated variations 
of yields of the classes the events migrating between. These two kinds of uncertainties are 
introduced below respectively. 


Uncertainties Causing 100% Correlated Variations of Yields 


The systematic uncertainties causing 100% correlated variations of yields of the event classes 


under influence are summarized in Table 7.3 The systematic sources are listed in the first 
column, and their corresponding uncertainties are listed in the second column. Among 
the uncertainties, the cross section uncertainty of each Higgs production process and the 
branching ratio uncertainty of Higgs decaying to two photons are associated with theoretical 
calculations. The former consists of two components: one is from the uncertainty of the 
Parton Distribution Functions (PDF); the other is from the effect of missing higher order 
correction terms, evaluated by varying the factorization scale and the renormalization scale 
(scale). For this analysis, the events from WH and ZH processes are considered together 
as events from VH, and the larger uncertainty of WH and ZH is taken. The rest of the 
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uncertainties are associated with experimental measurements. The theoretical uncertainties, 
especially the cross section uncertainty of the ggH process, drive the uncertainty of the 
expected total signal yield, and thus the uncertainty of the signal strength. 


Table 7.3: The systematic uncertainties causing 100% correlated variations of yields of all 
the event classes under influence. 


Source 

Uncertainty 

Cross Section 
ggH 

VBF 

WH 

ZH 

tiH 

PDF 8 TeV (7 TeV) Scale 8 TeV (7 TeV) 

+7.5%-6.9% (+7.6%-7.1%) +7.2%-7.8% (+7.1%-7.8%) 

+2.6%-2.8% (+2.5%-2.1%) +0.2%-0.2% (+0.3%-0.3%) 

+2.3%-2.3% (+2.6%-2.6%) +1.0%-1.0% (+0.9%-0.9%) 

+2.5%-2.5% (+2.7%-2.7%) +3.1%-3.1% (+2.9%-2.9%) 

+8.1%-8.1% (+8.1%-8.1%) +3.8%-9.3% (+3.2%-9.3%) 

Branching Ratio 

77 —)■ 77 

+5.0%/-4.9% 

Integrated Luminosity 

2.6% (2.2%) 8 TeV (7 TeV) 

Trigger Efficiency 

1 .0% 

Preselection Efficiency 
Per Photon 

1.0% (2.6%) Barrel (Endcap) 


Uncertainties Causing Migration of Events 


The systematic uncertainties causing migration of yields among classes are further divided 
into two groups. One group is related to the DiphotonBDT, and mainly causes the events 
to migrate among the untagged classes, or to migrate into/out of the selection range of the 
analysis which is DiphotonBDT > —0.78 (0.19) for events at 8 TeV (7 TeV). The other group 
is related to the tags of the Higgs production processes and causes the events to migrate 
among the tagged classes, or to migrate between the tagged classes and the untagged classes. 


The systematic uncertainties related to the DiphotonBDT are summarized in Table 7.4 


The uncertainty of each source is propagated to the variation of the DiphotonBDT distri¬ 
bution as already described in Section |4.5.3 The resulting relative yield uncertainty of any 
event class is evaluated as the change of the yield due to the variation, and the maximum 
uncertainty is shown. 

The main systematic uncertainties related to the tags of the Higgs production processes 
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Table 7.4: The systematic uncertainties related to the DiphotonBDT. 


Source 

Yield Uncertainty 

Per Event Class 
(Up To) 

IDBDT Shifting 0.01 

^5% 

(Te/E Scaling 10% 

^16% 

Diphoton Kinematics Varying Higgs px and rapidity 

^20% 


are summarized in Table |7.5[ For each source, the tagged classes under event migration and 
the corresponding migration mode, either among the tagged classes or between the tagged 
and the untagged classes, are shown in the second column. The maximum relative yield 
uncertainty of each relevant Higgs production process for a type of classes, e.g. VBF Dijet 
classes, is shown in the third column. Among all the sources, the uncertainty related to the 
production of additional jets in the events from ggH process has the largest effect on the 
event migration. This contributes to 30% ggH yield uncertainty for all the VBF Dijet classes 
and for the ttH Multijet class, through the ggH event migration between these classes and 
the untagged classes, and up to 14% additional ggH yield uncertainty for the VBF Dijet 
classes, through the event migration among themselves. 


7.3.2 Systematic Uncertainties Related to the Signal Shape 


The systematic uncertainties related to the signal shape include the uncertainties associated 
with photon energy scale and resolution, and the uncertainty of vertex efficiency. The former 
influences the mean and width of both the right vertex and wrong vertex components of the 
signal shape, while the latter influences the relative contributions of these two components. 
The types of systematic sources are listed in the hrst column of Table [7T , and the number 
of corresponding sources, if more than one, are denoted in parenthesis. For each type, the 
largest relative uncertainty of the signal shape parameters for an event class due to a single 
source is shown. A brief description of these systematic uncertainties is provided below. 
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Table 7.5: The systematic uncertainties related to the tags of the Higgs production processes. 


Source 

Class 

(from/to Class) 

Yield Uncertainty 

Per Event Class 
(Up To) 

Production of 

VBF Dijet 

(Untagged) 

30% 

ggH 



Additional Jets 

VBF Dijet 

(Other VBF Dijet) 

14% 

ggH 



in ggH 

tiH Multijet 

(Untagged) 

30% 

ggH 



Jet Energy Scale 

VBF Dijet 

(Untagged) 

10 % 

ggH 

4% 

VBF 

and Resolution 

VBF Dijet 

(Other VBF Dijet) 

6 % 

ggH 

1 % 

VBF 

Muon Selection 

VH Lepton 

(Untagged) 

0.4% 

VH 




ttH Lepton 

(Untagged) 

0 .2% 

tin 



Electron Selection 

VH Lepton 

(Untagged) 

0.4% 

VH 




ttH Lepton 

(Untagged) 

0 .2% 

ttH 



MET Selection 

VH MET 

(Untagged) 

3% 

VH 

4% 

Non VH 

B-jet Selection 

ttH Multijet 

(Untagged) 

1 % 

tin 

2 % 

ggH 


ttH Lepton 

(Untagged) 

1 % 

tin 




Uncertainties Associated With Photon Energy Scale And Resolution 


The systematic uncertainties associated with photon energy scale and resolution originate 
from the imperfect energy correction between data and Monte Carlo simulation using Z —)■ 
e’''e“ events, which are due to three factors. 

The hrst factor is the different effects on the photon and electron energy reconstructions 
from the imperfect Monte Carlo simulations, the photon/electron differences. The main dif¬ 
ferences come from the dehcits of material simulation in the regions before ECAL, effectively 
about 10% dehcit in the region \t]\ < 1 and 20% dehcit in the region I?]! > 1 from esimations, 
and contribute up to 0.2% relative uncertainty of the mean value of the signal shape for 


an event class as shown in Table 7.6 The rest of the differences come from the imperfect 
simulations of the electromanetic shower, and the variation of collection rate of scintilation 
lights with their emission location along the longitudinal direction of the crystal, the light 
collection nonuniformity. 

The second factor is the variation of the energy scale difference between data and the 
Monte Carlo simulation as a function of the particle energy, the energy scale nonlinearity. 
The average electron energy used for the derivation of energy scale correction is lower than 
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the average photon energy from the Higgs decay, since the Z boson mass is 91.2 GeV while 
the Higgs mass in the search range is from 115 GeV to 135 GeV. This energy difference 
contributes up to 0.2% relative uncertainty of the mean value of the signal shape for an 
event class. 

The third factor is the imperfect method for energy scale and resolution correction be¬ 
tween data and Monte Garlo simulation, the energy correction method. This leads to the 
independent energy scale uncertainty and energy resolution uncertainty of each category of 
photons classihed according to the photon location (barrel or endcap) and the shower shape 
(> 0.94 or < 0.94). For photons in the barrel from events at 8 TeV, there are two un¬ 
certainties associated with the resolution, one for the constant smearing term and the other 
for the energy dependent term. All together, there are 10 (8) independent single photon 
energy uncertainties for events at 8 TeV (7 TeV). The largest relative uncertainty of the 
mean value, or of the width, of the signal shape for an event class due to a single photon 
energy uncertainty is 0.04%, or 3%. There is an additional source of uncertainty associated 
with the imperfect simulation of the intrinsic distribution of Z e~^e~, the Z —)■ e’''e“ 
line-shape, which contributes to 0.01% relative uncertainty on the mean value. 

Uncertainty of Vertex Efficiency 

The vertex efficiency is corrected between data and Monte Garlo simulation using Z —)■ 
events. The uncertainty associated with the correction is 1.5% of the right vertex component 
fraction of the signal shape for an event class. 

7.3.3 Correlation of Uncertainties Among Event Classes 

The systematic uncertainties due to different sources are independent. For the systematic 
uncertainties due to the same source, the uncertainties related to the signal yield are 100% or 
—100% correlated among the 7 TeV and 8 TeV event classes under influence. For the signal 
shape, the uncertainties associated with the photon/electron differences, the Z —)■ line- 
shape and the vertex efficiency are 100% correlated among the 7 TeV and 8 TeV event classes. 
The uncertainties assocaited with the energy nonlinerity and the effect of energy correction 
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Table 7.6: The systematic uncertainties related to the signal shape. 


Source 

Shape Uncertainty 

Per Event Class 
(Up To) 

Photon/Electron Differences 

Material Before ECAL (2) 

Light Collection Nonuniformity 

Electromagnetic Shower 

0.2% Mean 

0.02% Mean 

0.05% Mean 

Energy Scale Nonlinearity 

0.2% Mean 

Energy Correction Method 

Single photon energy scale/resolution 
(8 for 7 TeV, 10 for 8 TeV) 

Z —)■ e+e“ line-shape 

0.04% Mean 

3% Width 

0.01% Mean 

Vertex Efficiency 

(w Right Vertex 
° Fraction 


method on single photon energy scale and resolution are 100% correlated within 7 TeV 
classes or 8 TeV classes. These uncertainties are not 100% correlated between the 7 TeV 
and the 8 TeV event classes since they are sensitive to the independent energy calibrations, 
regressions and the differences in the energy correction procedures of 7 TeV and 8 TeV 
events. There are 20% and 50% correlations assigned to the uncertainties associated with 
the energy nonlinearity and the effect of energy correction method on single photon energy 
scale between 7 TeV and 8 TeV classes, respectively, and no correlation assigned to the 
uncertainties associated with the effect of energy correction method on single photon energy 
resolution. 


7.3.4 Procedure to Incorporate Systematic Uncertainties 


The signal model as introduced in Section 7.1 and the corresponding likelihood function for 


each event class as introduced in Section |7.2| are modified to incorporate the signal yield and 
shape uncertainties through nuisance parameters, each associated with a particular source 
of systematic uncertainty. The procedure follows the description in References [111 ,115 and 
is introduced below. 
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Modification of Signal Yield 


The expected yield of a Higgs production process is modified as: 

n(0]v) 

= N^xHi^n) ■ J] (7.16) 

fc=i 

where 6^ represents the set of nuisance parameters associated with the sources of the sig¬ 
nal yield uncertainties, n{dM) represents the number of nuisance parameters, 9% represents 
the nuisance parameter associated with the kth source and 5^ represents the corresponding 
relative yield uncertainty of the process in the event class. 

Modification of Signal Shape 

The mean for the ith Gaussian component of the signal shape is modihed as: 

n{e^) 

9^, 9f,{^/s)) = ^ 

k=l 

'n.{S^(y/s)) _ 

+ (\/l-4(VJ)-«J(8)+Ci.(^).»;;(7)).i;(v^)}, (7.17) 

k=l 

where 9^ represents the set of nuisance parameters associated with the sources of the mean 
uncertainties 100% correlated between 7 TeV and 8 TeV classes, 9^{^/s) represents the set 
of nuisance parameters independent for \/s=7 TeV or y/s=8 TeV, n{9^) and n(6*^(v%)) 
represent the numbers of corresponding nuisance parameters, 9^ represents the nuisance pa¬ 
rameter associated with the kth source of 100% correlated uncertainties, represents the 
corresponding relative uncertainty on the mean, 9’^{7{8)) represents the nuisance parameter 
for 7 TeV (8 TeV) associated with the kth source of partially correlated uncertainties, S^{y/s) 
represents the corresponding relative uncertainty on the mean, Ck{^/s) represents the coef- 
hcient assoicated with the correaltion which is 0.2 (0.5) for uncertainties realted to energy 
nonlinearity (effect of energy correction method on single photon energy scale) aty%=8 TeV 
and is 1(1) at y/s=7 TeV. 

The standard deviation for the ith Gaussian component of the signal shape is modihed 
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as: 






n{8cri\/s)) 


k=l 


(«S(yS) ■ i5(v/i)P}. 


(7.18) 


where 6 * 0 - (a/s) represents the set of nuisance parameters associated with the sources of the 
width uncertainties for y/s=7 TeV or a/s=8 TeV, n{6rj{y/s)) represents the number of nui¬ 
sance parameters, 9'^{y/s) represents the nuisance parameter associated with the kth source 
and d^(\/s) represents the corresponding relative uncertainty on the width. 

The vertex selection efficiency is modihed as: 


enimH, Oy) = 6 ^( 771 ^) ■ min{(l -h 9v ■ Scr), 1}, (7.19) 

where 9v represents the nuisance parameter associated with the vertex selection efficiency 
uncertainty and 6eR represents the corresponding relative uncertainty. 


Modification of Likelihood Function 


A likelihood function, chosen as the standard Gaussian distribution, is assigned to each 
nuisance parameter. The likelihood function for the ith input background function dehned 
in Equation |7. 12 is modihed accordingly as: 


CiifXH, rriH, 9Bi, 9s) = A(/iH, m'^j, 0^/ ■ p{9s), 


(7.20) 


where 9s represents the set of signal nuisance parameters, and p{9s) represents the product 
of the likelihood functions of the nuisance parameters. The envelope function dehned in 
Equation |7. 13 is modihed accordingly as: 


-21n£g(/7/^,m/) = min {-2\nCi{fiH,mH,9 


D ^ ^ O ' 




(7.21) 


where 9^ ■ represents the set of values of the signal nuisance parameters maximizing the 
likelihood function with the ith background function at pn and a given m'jj. The generalized 
envelop function for any signal model S{m.y.y\pH) with Higgs parameters pn is 
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then defined as: 


-2\\iCe{,Vh) 


min {—21n£j {p ^, 6 


Bi ,PH 1 




(7.22) 


where OBi,pH ^s,i,pH represent the set of values of the background parameters and the set 
of values of the signal nuisance parameters maximizing the likelihood function with the ith 
background function at pn- 


7.4 Higgs Signal Extraction Procedure 


The Higgs signal is finally extracted by performing simultaneous profile likelihood fits to the 
observed diphoton mass spectra of all the 25 event classes, including 11 classes for 7 TeV data 
and 14 classes for 8 TeV data. The existence of a signal is demonstrated by a background 
only hypothesis test. The properties of the signal and its compatibility with the SM Higgs 
boson are quantified by measuring various Higgs parameters. The statistical method used is 
described in References 111,115, 1161 and introduced as below. 

For fitting the Higgs parameters, pn, with the associated signal model, S{m^^\pH), the 
parameters pn and the signal systematic nuisance parameters are varied simultaneously 
across all the event classes, while the background nuisance parameters are varied indepen¬ 
dently for each event class. The total envelope function, —2\nCTot{pH), is constructed as: 


25 

-2\nCTot{PH) = (7.23) 

i=l 

where —2lnC*^-{pH) is the envelop function for the ith class. From the total envelope function, 
the best fit pn, the values oi pn minimizing the function, and the associated confidence 
interval or region are extracted. For extracting the confidence interval or region for a subset 
of ph. Phi the remaining parameters of pui P%, are profiled as nuisance parameters, and the 
resulting likelihood ratio function qsipu) i® used: 




-21n- 


^TotiPH,P%pi^) 


C-Totip 


H 


(7.24) 
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where represents the values of maximizing CxotipH) at a given 

For testing the background-only hypothesis against the existence of a signal at Higgs mass 
hypothesis in the presence of an excess of events above the background-only expectation, 
the test statistic qb{rri^) is constructed as: 


(Ibimu) = -21n 




, Ph>0 OT gb(m^) = 0, pH < 0, 


(7.25) 


where pn is the pu maximizing the likelihood at a given m^. The proba¬ 

bility of the test statistic under the background only hypothesis, p{qb{fn^)\pH = 0), is 0.5 
for qb{m'j^) = 0, and follows 0.5 times the distribution with one degree of freedom for 
qbim'jj) > 0 in the limit of large number of events. The probability for observing equal or 
larger excess as the observed one under the background only hypothesis is then quantihed 
by the local p-value, and is translated into the local signihcance aiocai through the standard 
Gaussian distribution g(x): 


local p-value = / p{(lb{'mfj)\pH = 0) dqb{mfj) = / g{x) dx. 


(7.26) 


where ql^^ijn^) is the value of the test statistic observed from the data and aiocai is q'b^^ijn h) 
in the limit of large number of events. 


In the end of this analysis, the Higgs parameters introduced in Section are measured. 
The total signal strength pn is measured using the signal model S{m^^\pH,mH) with uih 
treated as a nuisance parameter. The Higgs mass, mn, the signal strength for ggH and 
tiH processes, PggHpH) ^md the signal strength for VBF and VH processes, Pvbf,vh cire 
measured using the signal model S{m^^\pggHtiH^ pvbf,vhi'^h)- For the measurement of 
each of the parameters, the rest two are treated as nuisance parameters. The Higgs coupling 
strengths to bosons and to fermions, Ky and k/, and the effective Higgs coupling strengths 
to photon and to gluon, and are measured using the signal models S{m^^\Bv, k/, mn) 
and Hg^rriH) respectively, at the measured mn- 
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Chapter 8 


Results of Higgs Search from CMS 
i/ —> 77 Channel 


8.1 Diphoton Mass Spectra and Fits 


The observed diphoton mass spectra are shown in Figure |8-1 Figure 8-2 and Figure 8-3 


for the 7 TeV classes, and in Figure [8^ Figure [8^ and Figure [8^ for the 8 TeV classes. 
A Higgs signal-like excess is observed and quantihed through the simultaneous signal plus 
background £t, using the signal model to the diphoton mass spectra over 

all event classes. The corresponding best-fit values of the signal strength and the Higgs mass 
are jin = 1-12 and = 124.72 GeV. For each event class, the signal plus background model 
at the best-£t (solid red line) is shown. The background component for the £t (dashed red 
line), along with the 68.3% (1 a) conhdence band (yellow) and the 95.4% (2 a) conhdence 
band (cyan) for the expected number of background events from the £t, is shown as well. 

More information for each event class, including the expected number of background 
events per GeV {dB/dm^^) at 125 GeV, and the expected S/B and S/y/B at mn = 125 GeV, 


is presented in Table 7.1, where the number of background events under the signal peak, B, 


is estimated as dB/dm^^ at 125 GeV multiplied by 4 Ceff. The S/B is higher for the tagged 
classes than for the untagged classes in general, and decreases with the increase of the class 
number for the untagged classes, as expected. The S/y/B provides a measure of the signal 
sensitivity of each event class, according to which the 8 TeV untagged 2 class is the most 
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sensitive class though not the one with the highest S/B^ as a result of its relatively large 
signal yield. 

The combined diphoton mass spectrum, with the corresponding signal plus background 


model, of all the 7 TeV and 8 TeV event classes is shown in Figure 8-7 The combined signal 
plus background model is obtained by summing the best-£t signal plus background models 
of all the event classes according to their fractions of the total number of events. The signal 
peak is not obvious because the signals in the high S/B classes are submerged by mixing 
with large number of background events from the low S/B classes. This is the reason that we 
classify events according to S/B and extract the signal by simultaneous £t to the diphoton 
mass spectra over all event classes, instead of fitting an combined diphoton mass spectrum, 
in order to achieve the best signal sensitivity. 

The weighted version of the combined mass spectrum, with the corresponding signal plus 


background model, is shown in Figure 8-8, which provides a better view of the observed 
signal-like excess. The data for each individual class is weighted by the ratio S/{S + B), 
which is evaluated using the values of the signal model and background model at the best- 
£t signal strength and the Higgs mass. A normalization factor is applied such that the 
total number of htted signal events keeps unchanged after the weighting. The signal plus 
background curve shown in the figure for the weighted spectrum is obtained by summing the 
best-fit signal plus background models of all the event classes according to their weighted 
fractions of the total number of events. The weighting is chosen according to the optimal 


signal extraction by htting to the weighted diphoton mass spectrum 117 . This weighting 
procedure estimates and visualizes the contribution of each event class in the simultaneous 
diphoton mass fit, though the fitting to the weighted diphoton mass spectrum is still not as 
optimal as the simultaneous fit used in this analysis for the signal extraction. 
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5.1 fb'' (7 TeV) + 19.7 fb"' (8 TeV) 



5.1 fb'(7TeV) + 19.7fb'(8TeV) ><.|q 3 5.1 fb'(7 TeV) + 19.7 fb'(8 TeV) 



Figure 8-1: The observed diphoton mass spectra of the untagged classes for the 7 TeV dataset 
(points) binned in 1 GeV steps. For each class, the signal plus background model (solid red 
line), at the best-£t fin = 1-12 and mn = 124.72 GeV associated with the signal model 
for the combined 7 TeV and 8 TeV datasets, is shown. The background 
component for the fit (dashed red line), the 68.3% (1 a) confidence band (yellow) and the 
95.4% (2 a) confidence band (cyan) are also shown. 
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Figure 8-2: The observed diphoton mass spectra of the VBF tagged classes for the 7 TeV 
dataset (points) binned in 1 GeV steps. For each class, the signal plus background model 
(solid red line), at the best-fit jlu = 1.12 and uih = 124.72 GeV associated with the signal 
model for the combined 7 TeV and 8 TeV datasets, is shown. The back¬ 

ground component for the £t (dashed red line), the 68.3% (1 a) confidence band (yellow) 
and the 95.4% (2 a) confidence band (cyan) are also shown. 
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Figure 8-3: The observed diphoton mass spectra of the VH and iiH tagged classes for the 
7 TeV dataset (points) binned in 1 GeV steps. For each class, the signal plus background 
model (solid red line), at the best-£t fin = 1T2 and mu = 124.72 GeV associated with the 
signal model 5(m.^.^|/iiy, mj^) for the combined 7 TeV and 8 TeV datasets, is shown. The 
background component for the £t (dashed red line), the 68.3% (1 a) conhdence band (yellow) 
and the 95.4% (2 a) conhdence band (cyan) are also shown. 
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5.1 fb'' (7 TeV) + 19.7 fb"' (8 TeV) 


5.1 fb' (7 TeV)+ 19.7 fb' (8 TeV) 




Figure 8-4: The observed diphoton mass spectra of the untagged classes for the 8 TeV dataset 
(points) binned in 1 GeV steps. For each class, the signal plus background model (solid red 
line), at the best-£t fin = 1-12 and rhu = 124.72 GeV associated with the signal model 
for the combined 7 TeV and 8 TeV datasets, is shown. The background 
component for the fit (dashed red line), the 68.3% (1 a) confidence band (yellow) and the 
95.4% (2 a) confidence band (cyan) are also shown. 
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5.1 fb"' (7 TeV) + 19.7 fb' (8 TeV) 



Figure 8-5: The observed diphoton mass spectra of the VBF tagged classes for the 8 TeV 
dataset (points) binned in 1 GeV steps. For each class, the signal plus background model 
(solid red line), at the best-£t fin = 1-12 and mn = 124.72 GeV associated with the signal 
model for the combined 7 TeV and 8 TeV datasets, is shown. The back¬ 

ground component for the fit (dashed red line), the 68.3% (1 a) confidence band (yellow) 
and the 95.4% (2 a) confidence band (cyan) are also shown. 
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5.1 fb'' (7 TeV) + 19.7 fb"' (8 TeV) 


5.1 fb' (7 TeV)+ 19.7 fb' (8 TeV) 





Figure 8-6: The observed diphoton mass spectra of the VH and iiH tagged classes for the 
8 TeV dataset (points) binned in 1 GeV steps. For each class, the signal plus background 
model (solid red line), at the best-£t fin = 1T2 and mu = 124.72 GeV associated with the 
signal model S'(m.^.^|/iiy, mj^) for the combined 7 TeV and 8 TeV datasets, is shown. The 
background component for the £t (dashed red line), the 68.3% (1 a) conhdence band (yellow) 
and the 95.4% (2 a) conhdence band (cyan) are also shown. 
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xio3 5.1 fb'^ (7 TeV) + 19.7 fb'^ (8 TeV) 



Figure 8-7: The sum of the observed diphoton mass spectra of all the event classes for the 
7 TeV and 8 TeV datasets (points) binned in 1 GeV steps. The corresponding signal plus 
background model (solid red line), obtained by summing the signal plus background models 
of all the event classes according to their fractions of the total number of events, is shown. 
The models correspond to the best-£t fin = 1-12 and tHh = 124.72 GeV associated with the 
signal model S{m^^\iiH,'rnH) for the combined 7 TeV and 8 TeV datasets. The background 
component for the combined model (dashed red line), the 68.3% (1 a) conhdence band 
(yellow) and the 95.4% (2 a) confidence band (cyan) are also shown. 
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Figure 8-8: The S/S + B weighted sum of the observed diphoton mass spectra of all the 
event classes for the 7 TeV and 8 TeV datasets (points) binned in 1 GeV steps. The cor¬ 
responding signal plus background model (solid red line), obtained by summing the signal 
plus background models of all the event classes according to their weighted fractions of the 
total number of events, is shown. The models correspond to the best-£t fin = 1.12 and ifin 
= 124.72 GeV associated with the signal model for the combined 7 TeV and 

8 TeV datasets. The background component for the weighted model (dashed red line), the 
68.3% (1 a) conhdence band (yellow) and the 95.4% (2 a) conhdence band (cyan) are also 
shown. 
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8.2 Local P-Value and Significance 


The local p-value of the background only hypothesis is scanned against the Higgs hypotheses 
in the range 115 GeV < mn < 135 GeV, in steps of 0.1 GeV. The observed local p-value 
and the corresponding signihcance of the excess as a function oi mn for the combined 7 TeV 
and 8 TeV datasets (solid black line), and the ones for the separate 7 TeV (solid blue 


line) and 8 TeV (solid magenta line) datasets are shown in Figure 8-9 The corresponding 
expected local p-value and local signihcance (dashed lines) under the SM Higgs hypotheses 


are also shown. The expected values at each mn are evaluated on an Asimov dataset 116 


a representative dataset following the expected distribution of the corresponding signal plus 
background model with = 1. For the generation of the Asimov dataset, the background 
model at the best-ht fin and ttih are used and the systematic nuisance parameters for the 
signal model are also set to the values at the best-ht. 

The minimum observed local p-value from the combined 7 TeV and 8 TeV datasets is 
7.0 ■ 10“® aX rriH = 124.7 GeV, which corresponds to an excess with a local signihcance of 
5.7 standard deviations. This result, strongly disfavoring the background only hypothesis, 
leads to the observation of a new diphoton resonance—the conventional threshold for an 
observation in particle physics is 5.0 standard deviations. The expected local p-value for the 
Higgs at rriH = 124.7 GeV is 8.5 ■ 10“®, corresponding to an excess with a local signihcance 
of 5.2 standard deviations. 

The observed and expected local signihcance at mn = 124.7 GeV for the 7 TeV, 8 TeV 
and combined 7 TeV and 8 TeV datasets are summarized in Table 18.11 


Table 8.1: The observed and expected local signihcance (Jiocai at = 124.7 GeV. 


^ local 

7 TeV 

8 TeV 

7 TeV + 8 TeV 

Observed 

4.5 a 

4.1 a 

5.7 a 

Expected 

2.1 0- 

4.8 a 

5.2 a 
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Figure 8-9: The local p-value (left axis) of the background only hypothesis and the corre¬ 
sponding signihcance (right axis) of the excess against the Higgs hypotheses in the range 
115 GeV < ruH < 135 GeV. The observed values from the combined 7 TeV and 8 TeV 
datasets (solid black line), and the ones from separate 7 TeV (solid blue line) and 8 TeV 
(solid magenta line) datasets are shown. The corresponding expected values (dashed line) 
are also shown. The excess corresponds to a signihcance of 5.7 standard deviations. 
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8.3 Overall Higgs Signal Strength 


The overall signal strength extracted from the combined 7 TeV and 8 TeV datasets is jlu = 
1.12lg'23 = 124.72 GeV, where the npper and lower nncertainties are the differences 

between the best-£t and the bonndary points of the 68.3% conhdence interval. This obtained 
signal strength is consistent with the SM Higgs expectation within the nncertainty. 

The observed contonr plot of likelihood ratio is shown on the left of Figure 


8-10 The best-£t (red cross), and the 68.3% (solid black line) and 95.4% (dashed black line) 


conhdence contours, correponding to = 2.3 and qs{fJ^H,mH) = 6.17 respectively, 

are also shown. The corresponding likelihood ratio qsi^n) treating mn as a nuisance param¬ 
eter obtained from the combined 7 TeV and 8 TeV (solid black line) datasets, and the ones 
obtained from the separate 7 TeV (solid blue line) and 8 TeV (solid magenta line) datasets 


are shown on the right of Figure 8-10 The boundary points for the 68.3% conhdence interval 
of hh correspond to qsifJ^n) = 1- In order to quantify separately the statistical uncertainty, 
including the uncertainty associated with the background model, and the systematic uncer¬ 
tainty, the qsi^n) with the signal systematic nuisance parameters hxed to the best-ht values 
for the combined 7 TeV and 8 TeV datasets (dashed black line) is obtained, from which 
the statistical upper and lower uncertainties are evaluated as -1-0.21/—0.21. The systematic 
upper and lower uncertainties are computed by subtracting the corresponding statistical 
uncertainties from the overall uncertainties in quadrature, which are -1-0.15/—0.09. 

The observed fiH and the corresponding friH for the 7 TeV, 8 TeV and combined 7 TeV 
and 8 TeV datasets are summarized in Table 18.21 


Table 8.2: The observed signal strength fiH and the corresponding mass mu. 



f^H 

mn 

7 TeV Observed 

2 iq+O-61 

124.03 GeV 

8 TeV Observed 

u.yu_Q 23 

124.93 GeV 

7 TeV -|- 8 TeV Observed 

1.1210:1 = 1.12l0;2}(stat)10;i5(syst) 

124.72 GeV 
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(GeV) 


CMS H^yy (Observed) 5.1 fb"' (7TeV) + 19.7 fb"' (8TeV) 



Figure 8-10: The observed likelihood ratio ^h) and qsil^n)- On the left, the observed 

likelihood ratio qs{fJ^H,fnH) from the combined 7 TeV and 8 TeV datasets is shown as a 
contour plot. The best-£t (red cross) is fin = 1-12 and mn = 124.72 GeV. The 68.3% 
conhdence contour (solid black line) and the 95.4% confidence contour (dashed black line) 
are shown. On the right, the observed likelihood ratio qsif^n) from the combined 7 TeV 
and 8 TeV dataset (solid black line), and the ones from the separate 7 TeV (solid blue line) 
and 8 TeV (solid magenta line) datasets are shown. The qsifJ^n) with the signal systematic 
nuisance parameters hxed to the best-fit values from the combined 7 TeV and 8 TeV dataset 
(dashed black line) is shown as well. The total uncertainty of the extracted signal strength 
from the combined 7 TeV and 8 TeV dataset is -1-0.26/—0.23, which consists of the statistical 
uncertainty -1-0.21/—0.21 and the systematic uncertainty -1-0.15/—0.09. 
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8.4 Mass 


The mass of the observed signal is extracted using the signal model with fi vbf, vh 

treated as nuisance parameters. The measured mass from the combined 7 TeV and 8 TeV 
datasets is mn = 124.72lg;3g GeV. 

The corresponding likelihood ratio gsimn) obtained from the combined 7 TeV and 8 TeV 
datasets (solid black line), and the ones obtained from the separate 7 TeV (solid blue line) 


and 8 TeV (solid magenta line) datasets are shown in Figure 8-11 For the combined 7 TeV 
and 8 TeV datasets, the gsirnu) with the signal systematic nuisance parameters fixed to 
the best-£t values (dashed black line) is also shown, from which the statistical uncertain¬ 
ties are evaluated as -1-0.31/—0.32 GeV. The corresponding systematic uncertainties are 
+0.16/-0.16 GeV. 

The observed mu for the 7 TeV, 8 TeV, and combined 7 TeV and 8 TeV datasets are 
summarized in Table 


Table 8.3: The results for the measurement of mass mn- 



mu 

7 TeV Observed 

124.19lK|^ GeV 

8 TeV Observed 

124.83lKj? GeV 

7 TeV -|- 8 TeV Observed 

124.721°;^ GeV = 124.72l°|^(stat)lK;(^(syst) GeV 
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CMS H^yy (Observed) 


5.1 fb'^ (7TeV) + 19.7 fb'^ (8TeV) 



Figure 8-11: The observed likelihood ratio gsimn) with figgH^tm fivBF,VH treated as 
nuisance parameters. The qs{mH) from the combined 7 TeV and 8 TeV datasets (solid black 
line), and the ones from the separate 7 TeV (solid blue line) and 8 TeV (solid magenta line) 
datasets are shown. The gsimn) with the signal systematic nuisance parameters fixed to the 
best-fit values from the combined 7 TeV and 8 TeV datasets (dashed black line) is shown 
as well. The best-fit is mu = 124.72 GeV. The total uncertainty of the measured mass is 
-1-0.35/—0.36 GeV, which consists of the statistical uncertainty -1-0.31/—0.32 GeV and the 
systematic uncertainty -1-0.16/—0.16 GeV. 
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8.5 Signal Strengths for Separate Higgs Production 


Processes 


The signal strength for ggH and itH processes extracted from the combined 7 TeV and 
8 TeV datasets is figgH,ttH = while the signal strength for VBF and VH processes 

is jlvBF,VH = Both obtained signal strengths are consistent with the SM Higgs 

expectation within the nncertainty. 

The observed likelihood ratio qs{f^ggH,ttHj f^VBF,VH) with uih treated as a nnisance pa¬ 


rameter is shown in Figure 8-12 The best-fit (red cross), the 68.3% (solid black line) and 
95.4% (dashed black line) confidence contours are also shown. The point (magenta triangle) 
corresponding to the SM Higgs expectation figgH,ttH = 1 fivBF,VH = 1 is within the 
68.3% conhdence contour. The corresponding qs{.l^ggH,ttH)i with mn and fivBF,VH treated as 
nuisance parameters, and qs{fJ^vBF,VH), with tjih and figgH,ttH treated as nuisance parameters, 
obtained from the combined 7 TeV and 8 TeV (solid black line) datasets, along with the ones 
obtained from the separate 7 TeV (solid blue line) and 8 TeV (solid magenta line) datasets. 


are shown on the left and right of Figure 8-13, respectively. For the combined 7 TeV and 
8 TeV datasets, the qs{,l^ggH,ttH) and qs{fJ^vBF,VH) with the signal systematic nuisance param¬ 
eters fixed to the best-£t values (dashed black line) are also shown, whose corresponding 
statistical uncertainties dominate the overall uncertainties. 

The observed i^ggH,ttH-i 1 ^vbf,vh and the corresponding mn for the 7 TeV, 8 TeV and 
combined 7 TeV and 8 TeV datasets are summarized in Table 18.41 


Table 8.4: The observed ggH and ttH signal strength figgH^tm and the VBF and VH signal 
strength P^vbf.vh along with the corresponding mass ttih- 








A VBF, VH 

niH 

7 TeV Observed 

1 4S+0-77 

4.18lJ|0 

124.19 GeV 

8 TeV Observed 

1 1 q+0.38 
-'■•-‘-^-0.34 


124.83 GeV 

7 TeV -|- 8 TeV Observed 

1 14+0.36 
-'■•■‘-^-0.31 

i-oslgi 

124.72 GeV 


To have a further look, the separate signal strengths for all the four production processes 
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Figure 8-12: The observed likelihood ratio 1^vbf,vh) with rriH treated as a nui¬ 

sance parameter from the combined 7 TeV and 8 TeV datasets. The best-fit (red cross) is 
figgH,ttH = 1-14 and jlvBF,VH = 1-08. The 68.3% confidence contour (solid black line) and 
the 95.4% confidence contour (dashed black line) are shown. The point (magenta triangle) 
corresponding to the SM Higgs expectation figgn^uH = 1 fivBF,VH = 1 is within the 68.3% 
conhdence contour. 
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CMS H^yy (Observed) 5.1 fb"' (7TeV) + 19.7 fb'' (8TeV) CMS H^yy (Observed) 5.1 fb"' (7TeV) + 19.7 fb"' (8TeV) 



Figure 8-13: The observed likelihood ratio {(ls{.fJ'VBF,VH)) with rriH and ^ivbf,vh 

{jJ^ggHfin) treated as nuisance parameters. On the left (right), the observed likelihood ratio 
Qs{fJ'ggH,ttH) {Qs{fJ^VBF,VH)) hom the combined 7 TeV and 8 TeV datasets (solid black line), 
and the ones from the separate 7 TeV (solid blue line) and 8 TeV (solid magenta line) datasets 
are shown. The (ls{fJ^ggH,ttH) {Qs{fJ'VBF,VH)) with the signal systematic nuisance parameters 
hxed to the best-£t values from the combined 7 TeV and 8 TeV datasets (dashed black 
line) is also shown. The best-£t is f^ggH,ttH = 1-14 {f^vBF,VH = 1-08). The uncertainty of 
the extracted signal strength from the combined 7 TeV and 8 TeV datasets is -1-0.36/—0.31 
(+0.62/-0.56). 
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^‘ggH, fJ'VBF, I^^VH and UtiH are also extracted. For the determination of each signal strength, 
the other three signal strengths and mn are treated as nnisance parameters. Since the domi¬ 
nant event classes for VH and tiH processes have low statistics, the accnracy of the obtained 
fiVH and fitiH along with their nncertainties snffer from the backgronnd estimation as men¬ 


tioned in Section 7.2.3 Instead of providing the most accnrate evalnations of individnal 
signal strength, the resnlts provide an overall estimation of the compatibility with the SM 


Higgs boson. The resnlts are listed in Table 8^ The largest deviation from the expectation 
of the Higgs boson is the signal strength of VH prodnction process, which is still compatible 
with the expectation within 2 standard deviations. 


Table 8.5: The observed signal strengths for all the fonr prodnction processes P - vbf -, 

jlvH and fittH along with the corresponding mass mn- 



kggH 

kvBF 

kvH 

kttH 

mu 

7 TeV -|- 8 TeV Observed 

1 1 i:-l-0.37 
-'■•-‘-^-0.32 

-'-•^■'--0.68 

— f) 

2.561?;™ 

124.60 GeV 


8.6 Higgs Coupling Strengths 


The likelihood ratios qs{K,v,F,f) and qs{B^,Bg) scanned at uih = 124.72 GeV from the com¬ 
bined 7 TeV and 8 TeV datasets, along with the corresponding best-hts (red cross), the 
68.3% (solid black line) and 95.4% (dashed black line) conhdence contonrs, are shown on the 


left and right of Fignre [8-14[ respectively. 

For the likelihood scan of Ky and Kf, it assnmes Ky > 0 as only the relative sign between 
Kv and K/ is measnrable. The best-£t is ky = 1.05 and kj = 1.03, which snpports the same 
sign scenario and is consistent with the SM Higgs expectation ky = 1 and kf = 1 (magenta 
triangle) at 68.3% conhdence level. The opposite sign scenario is not exclnded thongh, and 
the local mimimnm in the region Kf < 0 is within the 68.3% contonr. The 68.3% conhdence 
interval (CL) for Ky is [0.61, 0.77] U [0.90, 1.24], and that for Kf is [—0.95, —0.50] U [0.69, 
1.75]. 


166 












For and Kg, the extracted values are k^y = 1.10^q; 23 and kg = O.Qd^Q;^®, consistent with 
the SM Higgs expectation. 



Figure 8-14: The observed likelihood ratios qs{f^v,i^f) and qs{K^,Kg) at tiih = 124.72 GeV 
from the combined 7 TeV and 8 TeV datasets shown on the left and right. The best-£t (red 
cross) for the Higgs coupling strengths to bosons and to fermions is ky = 1.05 and kf = 1.03. 
The best-£t (red cross) for the effective Higgs coupling strengths to photon and to gluon is 
k^ = 1.10 and kg = 0.94. The associated 68.3% (solid black line) conhdence contours and 
95.4% (dashed black line) conhdence contours are also shown. The points (magenta triangle) 
corresponding to the SM Higgs expectation Ky = 1 and Kf = 1, and k^ = 1 and Kg = 1 are 
within the 68.3% conhdence contours of the best-hts. 
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Chapter 9 


Other CMS and ATLAS Higgs 
Results 


To provide an overall picture of the Higgs searches at LHC, other main results from CMS 
and ATLAS experiments using LHC Run I data are briefly summarized below. 


9.1 Signal Significance, Mass and Compatibility with 
SM Higgs in Terms of Signal and Coupling Strengths 

9.1.1 CMS Results 

For the Higgs searches at CMS experiment through main decay channels, the H ^ ZZ ^ 4i 


channel 118 reports the observation of a narrow resonance with a local signihcance of 


6.8 standard deviations. The measured mass is mn = 125.6 ± 0.4(stat) ± 0.2(syst) GeV— 
compatible with the measured mass from the —)■ 77 channel within 2 standard deviations, 
and the best-£t overall Higgs signal strength is /Ih = 0-93lo;23(stat)lQ;o9(syst), consistent 
with the SM Higgs expectation. The H —>■ H^+H^” —)■ 2i2iy channel |119| reports an excess 
of events above background with a local signihcance of 4.3 standard deviations at the Higgs 
mass of 125.6 GeV measured from the H —ZZ —>■ 4i channel, and the corresponding 
best-£t signal strength (Ih = 0.72tg;^g, consistent with the SM Higgs expectation as well. 
Besides the bosonic decay channels, the two fermionic decay channels H —>■ 120] and 
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H ^bb 


121 


report an excess with a local significance of 3.2 and 2.1 standard deviations for 


a Higgs mass of 125 GeV, and the corresponding best-fit signal strength (Ih = 0.78 ± 0.27 
and jin = 1-0 ± 0.5, respectively. The combination of these two channels 122] leads to 
the strong evidence for the 125 GeV Higgs decaying into down-type fermions with a local 
significance of 3.8 standard deviations, for which the corresponding best-fit signal strength 
is ph = 0.83 ± 0.24. To test the direct Higgs coupling to up-type top quark, a search for 


itH production 123 is performed by analyzing the events from the above decay channels 
and the two photon decay channel tagged according to the tiH signature, assuming a Higgs 
mass of 125.6 GeV. An excess with a local significance of 3.4 standard deviations is observed, 
and the best-fit signal strength is fitm = 2 . 8 lo! 9 , which is compatible with the SM Higgs 
expectation at 2 standard deviations level. 

In addition, searches are performed through H —)■ and H —)■ (analysis only 

performed on events at 8 TeV for channels (HI as well, despite their very small 

branching ratios and low sensitivity. The observed (expected) 95% GL upper limits on their 
branching ratio for a Higgs mass of 125 GeV—assuming the SM cross section—are 0.0016 and 
0.0019, corresponding to 7 . 4 ( 6 . 5 l^; 9 ) and 3.7 x 10^ times the SM value, respectively. Since 
the result from H —)■ is consistent with the SM Higgs expectation with a branching 

ratio 0.0632 ± 0.0036 larger than the limits for and e^e“, the leptonic couplings of 

the Higgs are shown as not flavour-universal as expected by the SM. Furthermore, a search 
is performed for the Higgs decaying into particles not interacting with the detector—the 
invisible decays {H —)• invisible) |125 , targeting the non-SM decay particles such as dark 
matter candidates. The observed data is consistent with the SM background expectation, 
and the observed (expected) 95% GL upper limit on the invisible branching ratio for the 
125 GeV Higgs is 0.58(0.44). 


For the combined GMS results 126 , the Higgs mass measured through both the H —)■ 
77 and H —)■ ZZ —)■ Ai channels is itih = 125.02](Q;27(stat)]'^g;}5(syst) GeV. The over¬ 
all Higgs signal strength—the relative Higgs production cross section with respect to the 
SM expectation—as well as the signal strengths for different production processes are ex¬ 
tracted at this mass combining the main decay channels, H —)■ 77 , H ^ ZZ ^ 4i, 
H —)■ W^W~ —)■ 2^2//, H ^ bb and H —)■ , with multiple Higgs production tags ex- 
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plored. The best-fit overall signal strength is fin = 1-00 ± 0.09(stat)lQ;Q7(theo) ± 0.07(syst), 
where the systematic uncertainty is further decomposed into the theoretical related compo¬ 
nent (theo) and the rest (syst). The best-£t signal strengths for the individual production 
processes are figgH = 0.85lo;i6, fivBF = flyn = 0.92+Q;3g and jxtm = 2-90lJ;94. Both 

the overall signal strength and the individual production signal strengths are compatible 
with the expectations of the SM Higgs—for the iiH signal strength, agreeing with the result 
from the dedicated iiH search as mentioned above, the compatibility is at about 2 standard 
deviations level. Furthermore, various Higgs coupling strengths are probed under different 
physics scenarios, using the inputs from the main decay channels as well as the H —)■ 

assuming no in- 


and H —)■ invisible channels. For the benchmark scenarios of Reference 45 


visible or undetectable Higgs decays, the best-hts for the Higgs coupling strengths to bosons 
and fermions are ky = 1.01 ±0.07 and kf = 0.87^'^g Jg, which supports the same sign scenario 
between By and k/ as expected by the SM. The data excludes the opposite sign scenario 
at the 95% CL while not at the 99.7% CL. For the effective Higgs coupling strengths to 
photon and to gluon, the best-hts are k^ = l.ld^g;}^ and kg = 0 . 891 q; 4 q. The above results 
from combination are consistent with the results from H ^ 'j'j channel alone, and with 
smaller uncertainties due to the extra constraints from the other decay channels. The full 


combined results are provided in Reference 126 , which are all compatible with the SM Higgs 
expectation. 


9.1.2 ATLAS Results 


For the corresponding results from ATLAS experiment 127,128 , all the main bosonic decay 


channels report observation of excess with signihcance beyond 5 standard deviations—5.2 
standard deviations from H —)■ 77 . Strong evidence for the Higgs coupling to down-type 
fermions is obtained with a significance of 4.5 standard deviations. The measured Higgs mass 
from the 77 —)■ 77 and 77 —)■ ZZ —)■ 47 channels is = 125.36 ± 0.37(stat) ±0.18(syst) GeV. 
At this mass, the best-£t overall signal strength from the 77 —)■ 77 channel is fin = 
1.17 ± 0.27, which agrees with the result from the CMS 77 —77 channel. Combining 
all the main decay channels, together with the 77 —)■ Z'y and 77 —)■ /i/x channels, the best-£t 
overall signal strength is jin = 1.18 ± 0.10(stat)l'iQ Q 7 (theo) ± 0.07(syst). The best-£t signal 
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strengths for the individual production processes are figgH = 1-23 -a 20 ) P'VBF — 1-23 ih 0.32, 
jj-VH = 0.80 ± 0.36 and fitm = 1-81 ± 0.80. Various Higgs coupling strengths are probed as 
well. In particular, the best-hts for the Higgs coupling strengths to bosons and fermions are 
kv = 1.09 ± 0.07 and kj = The best-hts for the effective Higgs coupling strengths 

to photon and to gluon are k^ = 1.00 ± 0.12 and kg = 1.12 ± 0.12 with the effective Higgs 
coupling strength to Z 7 , Kz-y, prohled. These combined results are summarized in the right 


column of Table 9T, and are compared with the CMS combined results summarized in the 
left column. 


ATLAS and CMS, with different detector design, independent analysis methods and 
similar luminosities for the analyzed data, obtain results compatible with each other, which 
lead to the observation of a new particle with the signal and coupling strengths consistent 
with the Standard Model Higgs boson. 


Table 9.1: The comparison between combined CMS results (left) and ATLAS results (right). 



CMS 

ATLAS 

ifiH (GeV) 

125.02;°;2®(stat)l[(;J^(syst) 

125.36 ± 0.37(stat) ± 0.18(syst) 

P-H 

1.00 ± 0.09(stat)tQQ7(theo) ± 0.07(syst) 

1.18 ± 0.10(stat)+o;o7(theo) ± 0.07(syst) 

PggH 

0.851°:^® 

1 (^q+0.23 

J--^o_0.20 

PVBF 


1.23 ±0.32 

PVH 

n no+O-SS 

'^•^^-0.36 

0.80 ±0.36 

PttH 

^•^'^-0.94 

1.81 ±0.80 

kv 

1.01 ±0.07 

1.09 ±0.07 


0.87l°;i^ 


k-y 


1.00 ±0.12 

kg 

omtlil 

1.12 ±0.12 
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9.2 Spin and Parity 


The new particle is identified as a boson since it is observed through the hT —?• 77 and 
H ^ ZZ ^ 4£ channels. Its observation through the TT —)■ 77 channel further indicates 


that its spin is not equal to 1 129,130 and its charge conjugation is positive. All observations 
are in favor of the SM Higgs spin-parity hypothesis with spin-0 and even parity, while disfavor 
opposite parity under spin -0 hypothesis, spin -1 hypothesis and several models under spin -2 


hypothesis tested so far 40,131-135 


173 










174 



Chapter 10 


Conclusion 


Passed by mornings and nights, bright and dark, we are now at the end of this 
Odyssey, searching for the Higgs boson through its decay into two photons at the 
CMS experiment at CERN’s Large Hadron Collider. This thesis concludes here with 
our hnal results concerning the observation of a new particle and the measurements of its 
properties from the refined and extended analysis, using the advanced multivariate analysis 
techniques, that we have developed since 2011, on the full LHC “Run I” data collected by 
the CMS detector during 2011 and 2012, consisting of proton-proton collision events at ^/s 
= 7 TeV with L = 5.1 fb“^ and at ^/s = 8 TeV with L = 19.7 fb~^, with the final calibration. 


An excess of events above the background expectation is observed, with a local signifi¬ 
cance of 5.7 standard deviations at a mass of 124.7 GeV. This result confirms our observation 
of an excess of events, with a local significance of 4.1 standard deviations near 125 GeV in 
2012, which provided the strongest evidence among all the Higgs search channels for the ob¬ 


servation of a new particle from the GMS experiment 36,37 . This result further constitutes 
the standalone observation of the new particle through the two photon decay channel. 


A further measurement provides the precise mass of this new particle as 
rriH = 124.72;°J^ GeV = 124.72;°J^(stat)l[J;}^(syst) GeV, 
with a relative total uncertainty less than 0.3% dominated by the statistical uncertainty. 
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The production cross section times the two photon decay branching ratio of this new 
particle relative to that of the Standard Model Higgs boson, the signal strength, for all the 
Higgs production processes combined, is extracted as 

= 1.12+°|}(stat)l[^;J^(syst). 

The relative uncertainty is about 20% dominated by the statistical uncertainty. This result 
is compatible with the Standard Model Higgs boson expectation within the uncertainty. 


The separate signal strengths for VBF and VH production processes, sensitive to Higgs 
couplings to bosons, and for the ggH and tiH prodnction processes, sensitive to Higgs cou¬ 
plings to fermions, are further extracted as 

^^VBF,VH — J-'Uo.o.se! 
n - — 114 + 0-36 

which have large uncertainties and are consistent with the Standard Model Higgs boson. 


The signal strengths for individnal prodnction processes are also extracted as 

r, — 1 1 c;+0.37 
h^ggH — -‘-•-‘-+-0.32) 

flVBF = 1 - 51 ^ o ; 68 , 




2.56 


+2.50 

-1.79- 


These results, especially for VH and tiH, are limited by the large statistical uncertainties. 
The largest deviation from the expectation of the Higgs boson is the signal strength of 
VH production process, which is still compatible with the expectation within 2 standard 
deviations. 
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The couplings of this new particle to bosons and to fermions relative to the key predictions 
from the Standard Model about those of the Higgs boson, proportional to boson mass squared 
and to fermion mass, respectively, assuming the existence of the Yukawa interactions between 
the Higgs boson and fermions, are further extracted as 

kv = 1-05 with 68.3% conhdence interval [0.61, 0.77] U [0.90, 1.24], 

kf = 1.03 with 68.3% conhdence interval [—0.95, —0.50] U [0.69, 1.75]. 

The extracted Ky shows that the coupling of the new particle to bosons is compatible with 
the Standard Model prediction at 68.3% conhdence level. The extracted Kf supports the 
existence of the interaction between the new particle and fermions, and further shows that 
the coupling of the new particle to fermions is compatible with the Standard Model prediction 
at 68.3% conhdence level. 

The ehective couplings of the particle to photon and to gluon relative to the Standard 
Model Higgs boson are extracted as 

— 1 1 n+0.21 
Ky — 1.1U_Q23, 

_ Q 04+0-38 
Kg — u.au_Q 23. 

These results are also compatible with the Standard Model Higgs boson expectation and 
provide no evidence for the existence of new heavy particles in the loops given the current 
precision. 


177 



The observation of a new particle from the if —77 channel is snpported by the hnal 
search results from the other main Higgs decay channels at the CMS experiment on the 
LHC “Run I” data, including the standalone observation of a new particle from the H —)■ 


ZZ —>■ channel 118 , and strong evidences from the H —)■ H^+H^ —)■ 2i2u channel 


119 


and also from the combination of fermionic decay channels H —)■ and H ^ bb 


122 . These results conhrm the observation of a new particle from the CMS experiment 


in 2012 37 . Combining the if —)■ 77 and if ZZ —)• 4f channels, the mass of the 
new particle is measured precisely as 125.02^Q;27(stat)lQ;}5(syst) GeV. Combining all main 
decay channels, the total production cross section of this new particle relative to that of the 
Standard Model Higgs boson is extracted as 1.00 ± 0.09(stat)lQ;o7(theo) ± 0.07(syst), with 
the relative uncertainty reduced to 10 % with respect to the result from the if —)■ 77 channel 
alone, and is compatible with the Standard Model Higgs boson expectation. All the other 
CMS combined results on the relative cross sections for separate Higgs prodcution processes 
and couplings are compatible with the Standard Model Higgs boson expectations as well. 
In particular, the coupling to bosons and that to fermions relative to those of the Higgs 
boson are extracted as 1.01 ± 0.07, with an uncertainty within 10%, and 0.87lo;J3, with an 


uncertainty of about 15%, respectively 126 


The above observation and measurements of a new particle from the CMS experiment 
are conhrmed by the results from the ATLAS experiment also at LHC—with different design 
of detector, independent analysis methods, and similar luminosity of analyzed data—in the 


ff —)■ 77 channel and all the Higgs decay channels combined 127,128 


The new particle is identihed as a boson since it is observed through the f7 —)■ 77 and 
ff —)■ ZZ —)■ 4f channels. Its observation through ff —>■ 77 channel further indicates that its 


spin is not equal to 1 129 130 and its charge conjugation is positive. All the studies regarding 


the spin and parity of this new particle are in favor of the SM Higgs boson hypothesis with 


spin-0 and even parity 40,131-135 . 
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Final Remarks 


Standing at the end of this Odyssey of searching for the Higgs boson at the Large Hadron 
Collider, we have a list of results on our hands, which points to a new particle looking very 
similar to the Higgs boson in terms of the production rate, couplings, and spin and parity. 

What is the ultimate reality behind it? 

Is this particle the quantum of the scalar field, slowing down the particles with masses 
such that they could get together to form the structures in the universe including ourselves, 
as described by the Standard Model of particle physics? Does it relate to the phenomena 
beyond the description of the Standard Model such as dark matter? Further measurements 
of this particle from LHC “Run H”, with the center-of-mass energy increasing to 13 TeV and 
total luminosity about 100 fb“^, and from other experiments in the future would provide 
more information to tell. 

What is sure for the moment- 

Searching for the Higgs boson does bring us together to experience a series of events in 
that space and time, becoming a fundamental part of our existence and a monumental part 
of human history. All the information we having obtained from the proton-proton collisions 
at the Large Hadron Collider, along with the epic efforts of generations of physicists and 
engineers from all over the world, all the sleepless nights, all the memorable moments, all 
the collisions among ourselves, all the beautiful minds and hearts, all the emotions, and all 
the stories, are folded into the results sent towards the future, passing through layers of 
time, as pairs of photons passing through layers of nights, as what we have received from 
our predecessors. 
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Appendix A 


Figures of Signal Model 


7 TeV Untagged 0 7 TeV Untagged 1 




7 TeV Untagged 2 7 TeV Untagged 3 




Figure A-1: The 7 TeV untagged classes’s diphoton mass spectra (points) and the htted 
distributions (red lines) of Monte Carlo —)■ 77 events at a Higgs mass of 125 GeV. 
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7 TeV Dijet 0 


7 TeV Dijet 1 



Figure A-2: The 7 TeV VBF tagged classes’s diphoton mass spectra (points) and the htted 
distributions (red lines) of Monte Carlo —)■ 77 events at a Higgs mass of 125 GeV. 
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7 TeV VH Lepton Tight 



(GeV) 


7 TeV VH MET 



7 TeV VH Lepton Loose 



7 TeV VH Dijet 



Figure A-3: The 7 TeV 
distributions (red lines) 


VH tagged classes’s diphoton mass 
of Monte Carlo 77 —)■ 77 events at a 


spectra (points) and the 
Higgs mass of 125 GeV. 
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7 TeV ttH Lepton + Multijet 



Figure A-4: The 7 TeV ttH tagged class’s diphoton mass spectrum (points) and the fitted 
distribution (red line) of Monte Carlo FT —)■ 77 events at a Higgs mass of 125 GeV. 
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Figure A-5: 
distributions 


The 8 TeV untagged classes’s diphoton mass spectra (points) and the htted 
(red lines) of Monte Carlo —)■ 77 events at a Higgs mass of 125 GeV. 
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Figure A-6: The 8 TeV VBF tagged classes’s 
distributions (red lines) of Monte Carlo H —)■ 


diphoton mass spectra (points) and the 
77 events at a Higgs mass of 125 GeV. 
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8 TeV VH Lepton Tight 8 TeV VH Lepton Loose 




8 TeV VH MET 



8 TeV VH Dijet 



Figure A-7: The 8 TeV VH tagged classes’s diphoton mass spectra (points) and the htted 
distributions (red lines) of Monte Carlo —)■ 77 events at a Higgs mass of 125 GeV. 
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8 TeV ttH Lepton 8 TeV ttH Multijet 



Figure A-8: The 8 TeV ttH tagged classes’s diphoton mass spectra (points) and the htted 
distributions (red lines) of Monte Carlo H ^ 'y'j events at a Higgs mass of 125 GeV. 
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Appendix B 


Variables for Higgs Production 
Tagging 

B.l Variables Related to Jets 

• '■ ratio between the scalar sum of tracks in the jet which 
match any of the pileup vertices and the scalar pt sum of all tracks in the jet. 

• average square of Ai? between the particle-flow can¬ 
didate momentum within the jet and the jet momentum weighted by the p^ of the 
particle-flow candidate. This measures the width of the jet. 

• the transverse momentum of the jet (leading, sub-leading). 

• r/t(i-2). pseudorapidity of the jet (leading, sub-leading). 

• nijj: the dijet mass. 

• Nj-. the number of jets. 

• Nb-j'- the number of b-jets 
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B.2 Variables Related to Electrons 


• d’^y-. the absolute impact parameter of the electron track with respect to its closest 
vertex in the transverse plane. 

• dl: the absolute impact parameter of the electron track with respect to its closest 
vertex in 2 ;. 


• PconvVtx'- the p-value for the vertex £t of the conversion matching the electron. 

• NMiss'- the number of missing hits before the first hit of the track. 


EleMVA: the identihcation score evaluated by a Multivariate Technique estimating the 
likelihood of being a prompt electron over the likelihood of being an electron from a 
jet 
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• V^OReipuCorrPFCombineO'i'- fh® pileup Corrected pt sum of particle-flow charged hadrons, 
neutral hadrons and photons within Ai? < 0.3 of the electron divided by the electron 
Pt- The pileup contamination is estimated and subtracted by pevent times an effective 
area. 


• 77 ®: the pseudorapidity of electron. 

• p^\ the transverse momentum of electron. 

• rriee'- the dielectron mass. 


B.3 Variables Related to Muons 

• NPixel- the number of hits in pixel detector. 

• NtrkL ayer'- the number of tracker layers with hits. 

• NMuonChamber'- the number of hits in muon chamber. 

• NMatching- the number of muon stations with muon segments matching the tracker 
track. 
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• the absolute impact parameter of the muon track with respect to its closest vertex 
in the transverse plane. 

• the absolute impact parameter of the muon track with respect to its closest vertex 
in 

• /^DF\ x^ divided by number of degrees of freedom for the global muon track £t. 

• ReiBetaPuCorrPFComhineOi- the pilcup Corrected pt sum of particle-flow charged hadrons, 
neutral hadrons and photons within Ai? < 0.4 of the muon divided by the muon px. 
The pileup contamination is estimated and subtracted by 0.5 times the px sum of the 
charged particle-flow particles within the cone associated with pileup vertices. 

• the pseudorapidity of muon. 

• pi^\ the transverse momentum of muon. 

• 'm/iti'- fhs dimuon mass. 


B.4 Variables Related to Transverse Missing Energy 


MET: the magnitude of transverse missing energy. 


B.5 Variables Related to Photons 


1^77 ~ 2 ^^ I- separation between the diphoton pseudorapidity and the average 


pseudorapidity of the dijet 106 


the separation in the azimuthal angle between dijet and diphoton. 


AR^^e'- AR between photon and electron. 


AR^^etrk- AR between photon and electron track. 


the photon-electron mass. 
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• AR^ n'. AR between photon and muon. 


• cos{6*): cosine of the angle between the diphoton momentum in the center of mass 
frame of diphoton-dijet and the total momentum of diphoton-dijet in the lab frame. 

• \A(j)^^ji\: the separation in the azimuthal angle between the diphoton and the leading 
jet. 

• \A(j)^^^MET\'- the separation in the azimuthal angle between the diphoton and MET. 


192 



Bibliography 


[1] S. L. Glashow, “Partial-symmetries of weak interactions,” Nucl. Phys., vol. 22, pp. 579- 
588, 1961. 

[2] S. Weinberg, “A model of leptons,” Phys. Rev. Lett., vol. 19, pp. 1264-1266, 1967. 

[3] A. Salam, “Weak and electromagnetic interactions,” in Elementary particle physics: 
relativistic groups and analyticity (N. Svartholm, ed.), p. 367, Almquvist & Wiskell, 
1968. Proceedings of the eighth Nobel symposium. 

[4] H. Politzer, “Reliable perturbative results for strong interactions?,” Phys.Rev.Lett., 
vol. 30, pp. 1346-1349, 1973. 

[5] D. Gross and F. Wilczek, “Ultraviolet behavior of non-Abelian gauge theories,” 
Phys.Rev.Lett., vol. 30, pp. 1343-1346, 1973. 

[6] F. Englert and R. Brout, “Broken symmetry and the mass of gauge vector mesons,” 
Phys. Rev. Lett., vol. 13, pp. 321-323, 1964. 

[7] P. W. Higgs, “Broken symmetries, massless particles and gauge helds,” Phys. Rev. 
Lett., vol. 12, pp. 132-133, 1964. 

[8] P. W. Higgs, “Broken symmetries and the masses of gauge bosons,” Phys. Rev. Lett., 
vol. 13, pp. 508-509, 1964. 

[9] G. S. Guralnik, C. R. Hagen, and T. W. B. Kibble, “Global conservation laws and 
massless particles,” Phys. Rev. Lett., vol. 13, pp. 585-587, 1964. 

[10] P. W. Higgs, “Spontaneous symmetry breakdown without massless bosons,” Phys. 
Rev., vol. 145, pp. 1156-1163, 1966. 

[11] T. W. B. Kibble, “Symmetry breaking in non-Abelian gauge theories,” Phys. Rev., 
vol. 155, pp. 1554-1561, 1967. 

[12] F. J. Hasert et ai, “Search for elastic muon neutrino electron scattering,” Phys. Lett. 
B, vol. 46, pp. 121-124, 1973. 

[13] F. J. Hasert et ai, “Observation of neutrino-like interactions without muon or electron 
in the Gargamelle neutrino experiment,” Phys. Lett. B, vol. 46, pp. 138-140, 1973. 


193 



[14] G. Arnison et al., “Experimental observation of isolated large transverse energy elec¬ 
trons with associated missing energy at ^/s = 540 GeV,” Phys. Lett. B, vol. 122, 
pp. 103-116, 1983. 

[15] M. Banner et al., “Observation of single isolated electrons of high transverse momentum 
in events with missing transverse energy at the GERN anti-p p collider,” Phys. Lett. 
B, vol. 122, pp. 476-485, 1983. 

[16] G. Rubbia, “Experimental observation of the intermediate vector bosons , W~, 
and Rev. Mod. Phys., vol. 57, pp. 699-722, 1985. 

[17] J. M. Gornwall, D. N. Levin, and G. Tiktopoulos, “Uniqueness of spontaneously broken 
gauge theories,” Phys. Rev. Lett., vol. 30, p. 1268, 1973. 

[18] J. M. Gornwall, D. N. Levin, and G. Tiktopoulos, “Derivation of gauge invariance from 
high-energy unitarity bounds on the s matrix,” Phys. Rev. D, vol. 10, p. 1145, 1974. 

[19] G. H. Llewellyn Smith, “High-energy behavior and gauge symmetry,” Phys. Lett. B, 
vol. 46, pp. 233-236, 1973. 

[20] B. W. Lee, G. Quigg, and H. B. Thacker, “Weak interactions at very high energies: 
the role of the Higgs-boson mass,” Phys. Rev. D, vol. 16, p. 1519, 1977. 

[21] ALEPH Gollaboration, GDF Gollaboration, DO Gollaboration, DELPHI Gollaboration, 
L3 Gollaboration, OPAL Gollaboration, SLD Gollaboration, LEP Electroweak Working 
Group, Tevatron Electroweak Working Group, SLD electroweak heavy flavour groups, 
“Precision electroweak measurements and constraints on the standard model,” CERN- 
PH-EP/2010- 095;arXiv:1012.2367 . 

[22] ALEPH, DELPHI, L3, OPAL Gollaborations, and the LEP Working Group for Higgs 
Boson Searches, “Search for the standard model Higgs boson at LEP,” Phys. Lett. B, 
vol. 565, pp. 61-75, 2003. 

[23] GDF and DO Gollaborations, “Gombination of tevatron searches for the standard model 
Higgs boson in the W~^W~ decay mode,” Phys. Rev. Lett., vol. 104, p. 061802, 2010. 

[24] Fermilab, “Run H handbook.” http://www-ad.fnal.gov/runH/index.html. 

[25] LHG Higgs Gross Section Working Group, https://twiki.cern.ch/twiki/bin/ 
view/LHCPhysics/CrossSections. 

[26] LHG Machine Outreach. http://lhc-machine-outreach.web.cern.ch/ 
Ihc-machine-outreach/collisions.htm, 

[27] L. Evans and P. Bryant, “LHG machine,” JINST, vol. 3, p. S08001, 2008. 

[28] GMS Gollaboration, “The GMS experiment at the GERN LHG,” JINST, vol. 3, 
p. S08004, 2008. 


194 




[29] ATLAS Collaboration, “The ATLAS experiment at the CERN large hadron collider,” 
JINST, vol. 3, p. S08003, 2008. 

[30] CMS Collaboration, “Search for the standard model Higgs boson decaying into two 
photons in pp collisions at y/s = 7 TeV,” Physics Letters H, vol. 710, pp. 403-425, 
2012 . 

[31] CMS Collaboration, “A search using multivariate techniques for a standard model 
Higgs boson decaying into two photons,” CMS-PAS-HIG-12-001, 2012. 

[32] CMS Collaboration, “Search for the standard model Higgs boson in the decay channel 
H —>■ ZZ —)■ Ai in pp collisions at ^/s = 7 TeV,” Phys. Rev. Lett., vol. 108, p. 111804, 
2012 . 

[33] CMS Collaboration, “Combined results of searches for the standard model Higgs boson 
in pp collisions at a/s = 7 TeV,” Phys. Lett. B, vol. 710, pp. 26-48, 2012. 

[34] ATLAS Collaboration, “Combined search for the standard model Higgs boson using 
up to 4.9 fb“^ of pp collision data at a/s = 7 TeV with the ATLAS detector at the 
LHC,” Phys. Lett. B, vol. 710, pp. 49-66, 2012. 

[35] CMS Experiment, “Life@CMS: The Higgs story - unblinding, top-up and seminar.” 
Available at http://www.youtube. coni/watch?v=gmpqakF7_ME, 

[36] CMS Collaboration, “Observation of a new boson with mass near 125 GeV in pp 
collisions at a/s = 7 and 8 TeV,” JHEP, vol. 06, p. 081, 2013. 

[37] CMS Collaboration, “Observation of a new boson at a mass of 125 GeV with the GMS 
experiment at the LHG,” Phys. Lett. B, vol. 716, p. 30, 2012. 

[38] ATLAS Gollaboration, “Observation of a new particle in the search for the standard 
model Higgs boson with the ATLAS detector at the LHG,” Phys. Lett. B, vol. 716, 

p. 1, 2012. 

[39] GMS Gollaboration, “Measurements of the new Higgs-like boson at 125 GeV in the 
two photon decay channel,” CMS Analysis Note, CMS-AN-2013-253, 2014. 

[40] GMS Gollaboration, “Observation of the diphoton decay of the Higgs boson and mea¬ 
surement of its properties,” Eur. Phys. J. C, vol. 74, p. 3076, 2014. 

[41] G. Hooft, “Renormalizable lagrangians for massive Yang-Mills helds,” Nuclear Physics 
B, vol. 35, no. 1, pp. 167 - 188, 1971. 

[42] F. Halzen and A. D.Martin, Quarks and Leptons. New York, USA: John Wiley & Sons, 
1984. 

[43] A. Bettini, Introduction to Elementary Particle Physics. Gambridge, UK: Gambridge 
University Press, 2008. 


195 


[44] H. V. Klapdor-Kleingrothaus and A. Staudt, Non-accelerator Particle Physics. Bristol 
; Philadelphia: Institute of Physics Publishing, 1995. 

[45] LHC Higgs Cross Section Working Group, “Handbook of LHC Higgs cross sections: 3. 
Higgs properties,” CERN Report CERN-2013-004, 2013. 

[46] The TEVNPH Working Group for CDF and DO Collaborations, “Updated combina¬ 
tion of CDF and DO searches for standard model Higgs boson production with up 
to 10.0 fb“^ of data,” arXiv:1207.0449 [hep-ex], FERMILAB-CONF-12-318-E, CDF- 
NOTE-10884, DO-NOTE-6348, 2012. 

[47] ALICE Collaboration, “The ALICE experiment at the CERN LHC,” JINST, vol. 3, 
p. S08002, 2008. 

[48] LHCb Collaboration, “The LHCb detector at the LHC,” JINST, vol. 3, p. S08005, 
2008. 

[49] J.-L. Caron, “Accelerator complex of CERN: an overview of all accelerators of CERN.” 
Available at https://cds.cern.ch/record/42384, 

[50] J. Friedman, T. Hastie, and R. Tibshirani, “Additive logistic regression: a statistical 
view of boosting,” Annals of Statistics, vol. 28, pp. 337-407, 2000. 

[51] J. H. Friedman, “Greedy function approximation: a gradient boosting machine.,” An¬ 
nals of Statistics, vol. 29, pp. 1189-1232, 2001. 

[52] A. Hoecker et al., “TMVA: Toolkit for multivariate data analysis,” 
arXiv:physics/0703039 [physics, data-an]. 

[53] 1. Antcheva et al., “ROOT - a G-I--I- framework for petabyte data storage, statistical 
analysis and visualization,” Computer Physics Communications, vol. 180, pp. 2499- 
2512, 2009. 

[54] GMS Gollaboration, “GMS physics : Technical design report volume 1: Detector per¬ 
formance and software,” CERN-LHCC-2006-001; CMS-TDR-8-1. 

[55] J. Anderson et al., “Snowmass energy frontier simulations,” arXiv:1309.1057. 

[56] GMS Gollaboration, “The GMS tracker system project : Technical design report,” 

CERN-LHCC-98-006 ; CMS-TDR-5. 

[57] GMS Gollaboration, “The GMS tracker : addendum to the technical design report,” 

CERN-LHCC-2000-016 ; CMS-TDR-5-add-l. 

[58] GMS Gollaboration, “Description and performance of track and primary-vertex recon¬ 
struction with the GMS tracker,” arXiv:1405.6569. 

[59] GDF Gollaboration, “The GDF-H detector: Technical design report,” FERMILAB- 
DESIGN-1996-01; PERMITAB-PUB-96-390-E. 


196 


[60] CMS Collaboration, “The CMS electromagnetic calorimeter project : Technical design 
report,” CERN-LHCC-97-033; CMS-TDR-l 

[61] CMS Collaboration, “Addendum to the CMS ecal technical design report: Changes to 
the cms ecal electronics,” CERN-LHCC-2002-027. 

[62] CMS Collaboration, “Energy calibration and resolution of the CMS electromagnetic 
calorimeter in pp collisions at ^/s = 7 TeV,” JINST, vol. 8, p. P09009, 2013. 

[63] CMS Collaboration, “The CMS hadron calorimeter project: Technical design report,” 

CERN-LHCC-97-031; CMS-TDR-2. 

[64] CMS Collaboration, “The CMS muon project : Technical design report,” CERN- 
LHCC-97-032; CMS-TDR-3. 

[65] CMS Collaboration, “CMS TriDAS project : Technical design report, volume 1: The 
trigger systems,” CERN-LElCC-2000-038; CMS-TDR-6-1. 

[66] CMS Collaboration, “CMS the TriDAS project : Technical design report, volume 2: 
Data acquisition and high-level trigger,” CERN-LHCC-2002-026; CMS-TDR-6. 

[67] R. Friihwirth, “Application of Kalman hltering to track and vertex htting,” Nucl. 
Instrum. Meth. A, vol. 262, pp. 444-450, 1987. 

[68] R. E. Kalman, “A new approach to linear hltering and prediction problems,” J. Eluids 
Eng., vol. 82(1), pp. 35-45, 1960. 

[69] R. E. Kalman and R. S. Bucy, “New results in linear hltering and prediction theory,” 
J. Eluids Eng., vol. 83(1), pp. 95-108, 1961. 

[70] K. Rose, “Deterministic annealing for clustering, compression, classihcation, regression 
and related optimisation problems,” Proceedings of the IEEE, vol. 86, pp. 2210-2239, 
1998. 

[71] R. Fruhwirth, W. Waltenberger, and P. Vanlaer, “Adaptive vertex htting,” CMS- 
NOTE-2007-008. 

[72] CMS Collaboration, “Description and performance of track and primary-vertex recon¬ 
struction with the CMS tracker,” arXiv:1405.6569. 

[73] M. Anderson et ai, “Review of clustering algorithms and energy corrections in ECAL,” 
CMS-IN-2010-008. 

[74] E. Meschi et ai, “Electron reconstruction in the CMS electromagnetic calorimeter,” 
CMS-NOTE-2001-034, 2001. 

[75] H. Liu, G.Hanson, and N. Marinelli, “Conversion reconstruction with tracker-only 
seeded tracks in CMS 900 GeV data,” CMS-AN-2010-039, 2010. 


197 



[76] CMS Collaboration, “Electron reconstruction and identification at ^/s = 7 TeV,” CMS- 
PAS-EGM-IO-OO 4 , 2010. 

[77] W.Adam, R.Friihwirth, A.Strandlie, and T.Todorov, “Reconstruction of electrons with 
the Gaussian-sum filter in the CMS tracker at LHC,” J. Phys. G: Nucl. Part. Phys., 
vol. 31, p. N9, 2005. 

[78] CMS Collaboration, “Performance of CMS muon reconstruction in pp collision events 
at = 7 TeV,” JINST, vol. 7, p. P10002, 2012. 

[79] CMS Collaboration, “Particle-flow event reconstruction in CMS and performance for 

jets, tans, and GMS-PAS-PFT-09-001, 2009. 

[80] CMS Collaboration, “Commissioning of the particle-flow reconstruction in minimum- 
bias and jet events from pp collisions at 7 TeV,” GMS-PAS-PFT-10-002, 2010. 

[81] M. Cacciari and G. P. Salam and G. Soyez, “The anti-fe let clustering algorithm,” 
JHEP, vol. 04, p. 063, 2008. 

[82] CMS Collaboration, “Identification of b-quark jets with the CMS experiment,” JINST, 
vol. 8, p. P04013, 2013. 

[83] CMS Collaboration, “Measurement of the inclusive W and Z production cross sections 
in pp collisions at ^/s = 7 TeV with the CMS experiment,” JHEP, vol. 10, p. 132, 
2011 . 

[84] P. Nason, “A new method for combining NLO QCD with shower Monte Carlo algo¬ 
rithms,” JHEP, vol. 11, p. 040, 2004. 

[85] S. Frixione, P. Nason, and C. Oleari, “Matching NLO QCD computations with parton 
shower simulations: the POWHEG method,” JHEP, vol. 11, p. 070, 2007. 

[86] S. Alioli, P. Nason, C. Oleari, and E. Re, “A general framework for implementing NLO 
calculations in shower Monte Carlo programs: the POWHEG BOX,” JHEP, vol. 06, 
p. 043, 2010. 

[87] S. Alioli, P. Nason, C. Oleari, and E. Re, “NLO Higgs boson production via gluon 
fusion matched with shower in POWHEG,” JHEP, vol. 04, p. 002, 2009. 

[88] P. Nason and C. Oleari, “NLO Higgs boson production via vector-boson fusion matched 
with shower in POWHEG,” JHEP, vol. 02, p. 037, 2010. 

[89] T.Sjostrand, S.Mrenna, and P. Z. Skands, “PYTHIA 6.4 physics and manual,” JHEP, 
vol. 05, p. 026, 2006. 

[90] G. Bozzi, S. Catani, D. de Florian, and M. Grazzini, “The q-p spectrum of the Higgs 
boson at the LHC in QCD perturbation theory,” Phys. Lett. B, vol. 564, p. 65, 2003. 


198 



[91] G. Bozzi, S. Catani, D. de Florian, and M. Grazzini, “Transverse-momentum resum¬ 
mation and the spectrum of the Higgs boson at the LHG,” Nucl. Phys. B, vol. 737, 
p. 73, 2006. 

[92] D. de Florian, G. Ferrera, M. Grazzini, and D. Tommasini, “Transverse-momentum 
resummation; Higgs boson production at the Tevatron and the LHG,” JHEP, vol. 11, 
p. 064, 2011. 

[93] LHG Higgs Gross Section Working Group, “Handbook of LHG Higgs cross sections: 2. 
Differential distributions,” GERN Report GERN-2012-002, 2012. 

[94] L. J. Dixon and M. S. Sin, “Resonance-continuum interference in the di-photon Higgs 
signal at the LHG,” Phys. Rev. Lett., vol. 90, p. 252001, 2003. 

[95] J. Alwall et al, “MadGraph 5 : going beyond,” JHEP, vol. 06, p. 128, 2011. 

[96] T. Gleisberg et al, “Event generation with SHERPA 1.1,” JHEP, vol. 02, p. 007, 2009. 

[97] GMS Gollaboration, “Measurement of the production cross section for pairs of isolated 
photons in pp collisions at y/s = 7 TeV,” JHEP, vol. 01, p. 133, 2012. 

[98] GMS Gollaboration, “Measurement of the differential dijet production cross section in 
proton-proton collisions at y/s = 7 TeV,” Phys. Lett. B, vol. 700, p. 187, 2011. 

[99] S. Agostinelli et al, “GEANT4 - a simulation toolkit,” Nuel. Instrum. Meth. A, 
vol. 506, p. 250, 2003. 

[100] M. Gacciari and G. P. Salam, “Pileup subtraction using let areas,” Phys. Lett. B, 
vol. 659, p. 119, 2008. 

[101] M. Oreglia, A study of the reaetions tjj' —)■ 77 '^. PhD thesis, Stanford University, 1980. 
SLAG Report SLAG-R-236. 

[102] J. Beringer et al. (Particle Data Group), “Review of particle physics,” Phys. Rev. D, 
vol. 86 , p. 010001 , 2012 . 

[103] GMS Gollaboration, “Determination of jet energy calibration and transverse momen¬ 
tum resolution in GMS,” JINST, vol. 6 , p. P11002, 2011. 

[104] M. Gacciari, G. P. Salam, and G. Soyez, “The catchment area of jets,” JHEP, vol. 04, 
p. 005, 2008. 

[105] M. Gacciari, G. P. Salam, and G. Soyez, “FastJet user manual,” CERN-PH-TH-2011- 
297, 2011. 

[106] D. L. Rainwater, R. Szalapski, and D. Zeppenfeld, “Probing color-singlet exchange in 
Z J- 2-jet events at the LHG,” Phys. Rev. D, vol. 54, p. 6680, 1996. 

[107] 1. W. Stewart and F. J. Tackmann, “Theory uncertainties for Higgs mass and other 
searches using jet bins,” Phys. Rev. D, vol. 85, p. 034011, 2012. 


199 



[108] K. S. Cranmer, “Kernel estimation in high-energy physics,” Comput.Phys.Commun., 
vol. 136, pp. 198-207, 2001. 

[109] CMS Collaboration, “Absolute calibration of the luminosity measurement at CMS: 
Winter 2012 update,” CMS-PAS-SMP-12-008, 2012. 

[110] CMS Collaboration, “CMS luminosity based on pixel cluster counting - summer 2013 
update,” CMS-PAS-LUM-13-001, 2013. 

[111] CMS Collaboration, “Combined results of searches for the standard model Higgs boson 
in pp collisions at ^/s = 7 TeV,” Phys. Lett. 5, vol. 710, p. 26, 2012. 

[112] P. D. Dauncey, G. J. Davies, M. Kenzie, N. Wardle, “Handling background shape 
function uncertainty as a nuisance parameter, with reference to Higgs to two photons,” 
CMS-AN-2013-162, 2013. 

[113] P. D. Dauncey, G. J. Davies, M. Kenzie, N. Wardle, “An application of the treatment 
of parametric model choice as a discrete nuisance parameter to the GMS 77 —)■ 77 
analysis,” CMS-AN-2013-230, 2013. 

[114] G. Gowan, Statistical Data Analysis. Oxford University Press, 1998. 

[115] ATLAS and GMS Gollaborations, LHG Higgs Gombination Group, “Procedure for the 
LHG Higgs boson search combination in summer 2011,” Tech. Rep. ATL-PHYS-PUB- 
2011-011, GMS-NOTE-2011-005, GERN, 2011. 

[116] G. Gowan, K. Granmer, E. Gross, and O. Vitells, “Asymptotic formulae for likelihood- 
based tests of new physics,” Eur.Phys.J.C, vol. 71, p. 1554, 2011. 

[117] R. Barlow, “Event classihcation using weighting methods,” J. Comp. Phys., vol. 72, 
p. 202, 1987. 

[118] GMS Gollaboration, “Measurement of the properties of a Higgs boson in the four-lepton 
hnal state,” Phys. Rev. D, vol. 89, p. 092007, 2014. 

[119] GMS Gollaboration, “Measurement of Higgs boson production and properties in the 
WW decay channel with leptonic hnal states,” JHEP, vol. 01, p. 096, 2014. 

[120] GMS Gollaboration, “Evidence for the 125 GeV Higgs boson decaying to a pair of tan 
leptons,” JHEP, vol. 05, p. 104, 2014. 

[121] GMS Gollaboration, “Search for the standard model Higgs boson produced in associ¬ 
ation with a W or a Z boson and decaying to bottom quarks,” Phys. Rev. D, vol. 89, 
p. 012003, 2014. 

[122] GMS Gollaboration, “Evidence for the direct decay of the 125 GeV Higgs boson to 
fermions,” Nature Physics, vol. 10, pp. 557-560. 

[123] GMS Gollaboration, “Search for the associated production of the Higgs boson with a 
top-quark pair,” JHEP, vol. 09, p. 087, 2014. 


200 



[124] CMS Collaboration, “Search for a standard model-like Higgs boson in the ^ and 

decay channels at the LHC,” arXiv:1410.6679;CERN-PH-EP-2014-243, 2014. 

[125] CMS Collaboration, “Search for invisible decays of Higgs bosons in the vector boson 
fusion and associated ZH production modes,” Eur. Phys. J. C, voh 74, p. 2980, 2014. 

[126] CMS Collaboration, “Precise determination of the mass of the Higgs boson and tests 
of compatibility of its couplings with the standard model predictions using proton 
collisions at 7 and 8 TeV,” arXiv:1412.8662; CERN-PH-EP-2014-288, 2014. 

[127] ATLAS Collaboration, “Measurements of the Higgs boson production and decay rates 
and couplings using pp collision data y/s = 7 and 8 TeV in the ATLAS experiment ,” 
ATLAS-CONP-2015-007, 2015. 

[128] ATLAS Collaboration, “Measurement of the Higgs boson mass from the 77 —)■ 77 and 
H —)■ ZZ* —)■ 4/ channels in pp collisions at center-of-mass energies of 7 and 8 TeV with 
the ATLAS detector,” Phys. Rev. D, voh 90, 052004, 2014. 

[129] L. D. Landau, “On the angular momentum of a two-photon system,” Dokl. Akad. 
Nauk, voh 60, p. 207, 1948. 

[130] C. N. Yang, “Selection rules for the dematerialization of a particle into two photons,” 
Phys. Rev., voh 77, p. 242, 1950. 

[131] CMS Collaboration, “Measurement of the properties of a Higgs boson in the four-lepton 
hnal state,” Phys. Rev. D, vol. 89, p. 092007, 2014. 

[132] CMS Collaboration, “Measurement of Higgs boson production and properties in the 
WW decay channel with leptonic final states,” J. High Energy Phys., vol. 01, p. 096, 
2014. 

[133] CMS Collaboration, “Study of the mass and spin-parity of the Higgs boson candidate 
via its decays to Z boson pairs,” Phys. Rev. Lett., vol. 110, p. 081803, 2013. 

[134] ATLAS Collaboration, “Evidence for the spin-0 nature of the Higgs boson using AT¬ 
LAS data,” Phys. Lett. B, vol. 726, p. 120, 2013. 

[135] CMS Collaboration, “Constraints on the spin-parity and anomalous HVV couplings 
of the Higgs boson in proton collisions at 7 and 8 TeV,” arXiv:1411.344k, CERN-PH- 
EP-2014-265, 2014. 


201 



